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Abstract 



of "Recent approaches to high-dimensional American and Bermudan option pricing", 
by Frcdcrik S Herzberg, Morton College, University of Oxford 

A number of Bermudan option pricing methods that are applicable to options on 
multiple assets are studied in this thesis, one of the dominating questions being the 
natural scaling needed to extrapolate from Bermudan to American (both approximate 
and "exact") option prices. Among the Bermudan option pricing techniques discussed 
in more detail will be (1) the use of cubature formulae for symmetric measures to price 
non-perpetual Bermudan options, and (2) rcduitc-bascd approximation of subharmonic 
functions (corresponding to piecewise harmonic interpolation in the one-dimensional 
setting). These (iterative) algorithms shall be proven to be sound in a sense that is yet 
to be made precise. Also, for each iteration sequence resulting from these algorithms, 
monotone convergence to the least fixed point of the iteration procedure will be shown. 
(As part of an appendix, the fixed points of a Bermudan option pricing algorithm based 
on polynomial interpolation shall be characterised as well.) 

If the iteration procedure is based on cubature, a linear convergence rate of the 
iteration sequence can be derived. Moreover, for a couple of practically relevant set- 
tings one can find estimates for a non-perpetual American option pricing algorithm 
based on cubature. 

At the outset of this dissertation, the existence of exercise regions for multi- 
dimensional Bermudan options is established; afterwards one can proceed to prove 
bounds on the American-Bermudan barrier put option price difference ("continuity 
correction") when the argument of this function - as a function of the logarithmic 
start price - approaches the exercise boundary. In particular, results of Feller's shall 
be generalised to show that an extrapolation from the exact Bermudan prices to the 
American price cannot be polynomial in the exercise mesh size in the setting of many 
common market models, and more specific bounds on the natural scaling exponent of 
the non-polynomial extrapolation for a number of (both one- and multi-dimensional) 
market models will be deduced. 

Finally, three approximate A-hedging algorithms for high-dimensional derivative 
securities are proposed and implemented, alongside with a measure of comparing their 
effectiveness. 
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Chapter 1 



Definitions and basic facts on 
Bermudan and American 
options 



In order to clarify terminology, we start by introducing the mathematical notions 
corresponding to the financial concepts that we shall allude to. 
Our first definition is a notational convention. 

Definition 1.1. Let d € N. By exp : R'' ^ R>o'* and In : R>o'' we denote com- 

ponentwise exponentiation and taking natural logarithms componentwise, respectively. 

Remark 1.1. For any d G N, R'' is a Lie group with respect to componentwise multi- 
plication ■ : {x,y) 1-^ {xiyi)i£{i^,,,^d} ■ Its Lie algebra is the vector space M.'^ with its usual 
(componentwise) addition. The exponential map from the Lie algebra (R'', +) into the 
Lie group (R'',-) is componentwise exponentiation exp : x i— > (e^')^^^-^ Therefore 
the abbreviation introduced in Deftnition \l.l\ is consistent with standard notation. 

Definition 1.2. Let T be a positive real number. Consider a real-valued stochastic 
process X :— {Xt)t£[o,T], adapted to a filtered probability space (fl, {!Ft)te[o,T], P) ■ We 
will call X a logarithmic price process for a non-dividend paying asset (for short, 
a logarithmic price process or simply log-price process^, if and only if there exists 
a probability measure Q equivalent to P on Tt and a constant r > such that the 
stochastic process exp [Xt — rt^^^y^ is a martingale with respect to the filtration T := 
i^t)te[o,T] and the probability measure Q. Ln this case, such a Q is called a martingale 
measure and r a market price of risk or a dicsount rate for the stochastic process X 
and the probability measure P. 

Definition 1.3. Let d G N. A c?-dimensional basket is a d-tuple of logarithmic price 
processes such that there exists a probability measure Q and a market price of risk r > 
such that Q is a martingale measure and r a market price of risk for all components 
of the d-tuple. 

For the rest of this Chapter, we will adopt the terminology and the notation for 
Markov processes of Revuz and Yor pB^ . 




In particular, for all probability measures v on B {R'^) , is the probability mea- 
sure induced by the transition function ((5s)s>q via the lonescu-Tulcea-Kolmogorov 
projective limit construction, cf Revuz and Yor |2(il Theorem 1.5]). 

For any d € N, we will denote the cr-algebra of Borel subsets of R'' by B (W^) . 

Definition 1.4. Let again d N. A family Y :— (i^^)^gRd of M.'^ -valued homogeneous 
Markov processes adapted to a filtered probability space (n,J-,Q] with respect to 




5 



T , with transition function (-Ps)s>o and initial measure Sx, is called a d-dimensional 
Markov basket if and only if there is a homogeneous transition function (Qs)s>o on 
the measurable space (W^,B(M.'^)^ and a constant r > such that the following three 
assertions hold: 

1. The process is a Markov process with transition function {Qs)s>o '^iih respect 
to T for all x e . 

2. The process exp {Y^t — ''^)tg[o t] martingale with respect to T and Qs^. 

3. The measures Ps^ and :— Qs^ are equivalent for all x ^W'- . 

In this case, P is called a family of martingale (or: risk-neutral) measures associated 
with Y , and r is called the discount rate for Y . 

The expectation operator for the probability measure will be denoted by for 
all X £ R''. 

If the transition function P is a Feller semigroup, then we shall refer to Y as a 
Feller basket. 

If P is a translation-invariant Feller semigroup, we shall call Y a Levy basket. 

Remark 1.2. A priori, it is not clear if there are logical connections between the 
three assertions in the previous Definition \1.4\ in particular the author does not know 
whether the third assertion implies the first one. 

Notational convention 1.1. If no ambiguity can arise, we will drop the superscript 
of a Markov basket. Thus, in the notation of Definition \1.4[ we set 

E^[f{Y,„..., Y,J\T.] [/ (y;^, . . . , Yl) \ Ts] 

for all s > 0, n £fi and n-tuples of stopping times t = (ti, . . . , t„) whenever f : R" 
M is nonnegative or / {Y^^, . . . ,Y^^^ e (P^). Here we are using the term "stopping 
time" as a synonym for M.^-valued stopping time, that is a stopping time with values 
in [0, +oo]. 

Also, since we are explicitly allowing stopping times (with respect to the filtration 
generated by a process X) to attain the value +oo, we stipulate that the random variable 
f (Xr) (for any Lebesgue-Borel measurable function f) should be understood to be 
multiplied by the characteristic function of the event {r < +oo}. Formally, this can be 
done by introducing a constant A ^ , called cemetery, and stipulating that Xt = A 
on {t = +oo\ and /(A) = for all measurable functions f (cf eg Revuz and Yor ]2f)\ 
pp 84,102]). 

We will not formally define what we mean by an option itself, but we will rather 
define what expected payoffs and prices of some classes of financial derivatives are. 

Definition 1.5. Consider a d-dimensional Markov basket Y with an associated family 
P' of martingale measures and discount rate r > 0. 

The expected payoff of a Bermudan option with (log-price) payoff function g -.W^ 
M>o on the underlying Markov basket Y with exercise times in J C [0, +oo), log start- 
price X and maturity T e [0, +oo] is defined to be 

U'(T){x)= sup E4e-'-(^^^)ff(X,^T) 

r stopping time., r(0)C JU{-t-oo} 

The expected payoff of a perpetual Bermudan option is the expected payoff of a 
Bermudan option of maturity +00. 

The expected payoff of a Bermudan option with exercise mesh size h > is the 
expected payoff of a Bermudan option with exercise times in h -Nq . 
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The expected payoff of an American option is the expected payoff of a Bermudan 
option with exercise times in [0, +00). 

We shall call the expected payoff of a Bermudan option (or an American option) 
a Bermudan option price ( or an American option price ) if and only if the martingale 
measures associated with the underlying basket are unique (that is, if the market model 
described by P, T and X is complete^. 

In recent years, there has been increasing interest in incomplete market models 
that are governed by general Levy processes as log-price processes, as is not only wit- 
nessed by a tendency in research papers to focus on Levy process settings (for instance 
Boyarchenko and Levendorskii 0; Asmussen, Avram and Pistorius 0; 0ksendal and 
Proske pS], to take a random sample). Even textbooks, such as Karatzas' Wf\ and 
Mel'nikov's |22i] introductory works, are putting considerable emphasis on incomplete 
markets. Finally, "Levy finance" has already been treated in survey articles intended 
for a general mathematical audience, e g Applebaum's article [2- We will try not to 
deviate too much from this consensus that tries to accomplish as much mathematical 
generality as possible, while stopping short of studying Markov process models in their 
full generality. Instead we note that a substantial proportion of our results is concerned 
with perpetual Bermudan and American options, and it is precisely the medium and 
long-term risk theory where Levy finance seems to be applied most frequently. As a 
last remark on this issue, we consider it as beyond the scope of this thesis to ques- 
tion whether it is reasonable from an economist's point of view to study incomplete 
markets. 

Whilst there are some points to be made about market failures on stock markets 
that might entail arbitrage opportunities (for example, when assets are traded simulta- 
neously on several stock exchanges, or in the event of insider trading), the transaction 
costs to exploit these arbitrage opportunities usually tend to be close to the actual 
gain that can be achieved through taking advantage of the arbitrage. Therefore we 
shall, for the sake of mathematical simplicity, merely refer to the works of Corcos et 
al PI as well as Imkeller et al |14[ IT^ , and impose a strict no-arbitrage assumption 
(which under certain regularity conditions on the basket is equivalent to the existence 
of a martingale measure, cf Karatzas J3 Theorem 0.2.4]). 

Example 1.1 (A few common examples). 1. The price of a European call op- 
tion on a single asset with maturity T and strike price K is the price of a Bermu- 
dan option with the set of exercise times being the singleton {T} and the (log- 
price) payoff function (exp(-) — K) V 0. 

2. The price of a perpetual American put of exercise mesh size h > on the 
arithmetic average of two assets in an underlying basket with strike price K 
is the price of a Bermudan put option with the set of exercise times being the 
whole of the half-line [0,+oo), the maturity being T and the payoff function 

( X — °^p((')i)+°^p((')^) ^ V 



3. Consider a perpetual Bermudan call option on a single asset that continuously 
pays dividends at a rate S and whose logarithm follows a Markov process Z 
adapted to some probability space {fl, {J-'t)t>o, P)- Then, in order to exclude ar- 
bitrage, we will have to require the existence of a family of measures P' such that 
each is equivalent to P^ (in particular, P^^ = 6^) and such that (e~'"*+'^*+^*)^^ 

is a ¥^ -martingale for all a; e M''. The expected payoff' of the option will then be 
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sup 

r stopping time, r(^^)ChNoU{+oc-} 




As an auxiliary result, let us remark 
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hence 



Lemma 1.1 (Lower semi-continuity of sup). If I is a set and («/c,^)^g/ ^gj^Q is a 
family of real numbers, then 

sup lim inf afe^£ < liminf supa^,^. 

Proof. We have 

sup ak,e > afe/o 
i 

for all fc e No and £o G I, therefore for all n e N and £o G I, 

inf supofe,^ > inf ak,eo, 

k>n I k>n 

thus 

inf supafe,£ > sup inf ak,io, 

k>n i k>n 

sup inf sup ak^i > sup sup inf ak,eo = sup sup inf ak,eo ■ 

n k>n £ ' n eo ' to " 

This is the assertion. □ 

This estimate enables us to prove the following Lomma that is asserting the approx- 
imability of expected payoffs or prices of American options by sequences of expected 
payoffs or prices of Bermudan options, respectively. 

Lemma 1.2. Let d G f^, T > 0, x G W^, consider a bounded continuous function 

g > (the payoff function^, and a d-dimensional basket X having a modification with 
continuous paths. If the expected payoff of an American option of maturity T , log 
start-price x and payoff function g on this basket X is less than infinity, then the limit 

limf/™°(r)(:c) = SM^U^''^''{T){x) 
HO ken 

exists and equals the American expected payoff. 

Proof. Consider a sequence {hk)keHo € (lR>o)^° such that /ifc J, as A; — > oo. Choose a 
sequence of stopping times (t^)^^^^ such that for all x gW^, 

sup [e-'-f"^^)^ (^rAT)] = supE- [e-'-(^^^^)ff (X,,at)^ 

r stopping time ^ 

and define 

Ti,k := inf {t e /ifcNo : t > Tf } . 

Then, due to the continuity conditions we have imposed on g and on the paths of (a 
modification of) the basket X, we get 

supe-'-("^^^)5(^T,AT)=sup lim e-'^^^^.'^^^^^ (^r, ,at) 

f f fe^oo 

and hence by the lower semi-continuity of sup, one obtains 

supe-''(^^^^)5(X,,^T) < liminf supe-''(^^.'=^^)ff (X,, ,at) • 

Now we can use the Montone Convergence Theorem and Lebesgue's Dominated Con- 
vergence Theorem (this is apphcable because of the boundedness of g) to swap hm- 
its/suprema with the expectation operator. Combining this with the specific choice of 
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the sequence (r^j^gpj , this yields for all x € W^, 



sup 

T stopping time 



-r(TAT) 



g {Xtat) 



= supE^ 



supe 



-r{TeAT) 



9 {XreAT) 



liminf supe-'-^"^-"''^).? (X^,.,at) 



k^oo 



< liminfE^ 

k — '■oo 



liminf supE" le-"-^^'-"^^^ g (X,, ,at) 



< lim inf 



sup 

T(n)C?n,NoU{+oo} 



< lim sup sup E' 

fe^oo T(a)C?ifcNoU{+oo} 



< 



sup E^ 

T stopping time 



This finally gives 



sup E^ 

T Stopping time 



-r(TAT) 



hm stip E^ 



9 (Xtat) 



e-(-^^)g{XrAT)\=U'''''°{T){x). 



Since the left hand side does not depend on {hk)k, we conclude that lim/jjo U'^^^°{T){x) 
exists and is equal to sup^ stopping time [e-^^'^^'^^g (X^at)] • □ 



1.1 Some classes of Bermudan option pricing algo- 
rithms 

Let C° {R^, [0, +oo)) and L° {R^, [0, +oo)) , as usual, denote the spaces of nonnegative 
continuous functions defined on ffi'', and of nonnegative measurable functions defined 
on Mf^, respectively. 

The purpose of the following definitions is merely to introduce a fagon de parler 
which will allow us to quickly describe desirable properties of approximative Bermudan 
pricing algorithms in the later parts of this thesis. 

Definition 1.6. A map D : L° {R'^, [0, +oo)) L° (R'', [0, +00)) is said to be a sound 
iterative Bermudan option pricing algorithm (for short, a sound algorithm ) for a payoff 
function g € C° (M'', [0, +00)) if and only if Df > g for all f £ L° {R'^, [0, +00)) and 
the map D is pointwise monotone, that is 

V/o,/iGL° (R^[0,+cx))) 

(V.T e m." fo{x) < fi V.T e Dfo{x) < Df,{x)) . 

A sound iterative Bermudan option pricing algorithm D is said to have a perpetual 
limit if and only if 

sup D°"5 e i° (R'^, [0, +00)) 

nSNo 

(rather than this supremum being allowed to equal +00 on a subset of positive measure 
of its range). In that very case, the function in the last line is simply referred to as 
the perpetual limit of the algorithm. Finally, D is said to converge linearly in to 
the perpetual limit if and only if there exists a c € (0, 1) such that 



< c - II 



n-l\ 
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Z?" being shorthand for D°" for all n e Nq. 

Remark 1.3. The elements 0/ L° (M'', [0, +00)) should be conceived of assigning the 
value - that is, the expected payoff - of an option to the vector of logarithmic start 
prices of the components of the basket ( at least on the complement of a Lebesgue null 
set). 

Remark 1.4. The monotonicity condition imposed on sound iterative Bermudan pric- 
ing algorithms entail that the sequence of functions (£'"5)„gNo ~ i-^°"9)neNo (^^w^V^ 
pointwise increasing. Thus, this sequence has a limit: 

g < snp D''g = lim D'^g. 

The infimum of all D-fixed points is always an upper bound for the perpetual limit: 

Lemma 1.3. Let D be a sound iterative Bermudan pricing algorithm for g with a 
perpetual limit u. Then the function u is smaller than any fixed point of D; moreover, 
Du > u. 

Proof. Any fixed point /i of D is in the image of D and therefore, due to our assump- 
tions on sound algorithms, pointwise greater or equal g. Now, as D (and thus Z?") is 
pointwise monotone, 

Vn € No h = D''h > D''g, 

therefore 

h = sup D'^h > sup D"^g, 

m m 

where the right hand side is just the perpetual Hmit. Hence, any fixed point of D is 
greater or equal the perpetual limit. Furthermore, observe that due to the pointwise 
monotonicity of D, 

D (^snp D^^g^>D"^+'g, 

therefore for all m G No, 

D (supL'"g ) > snp D'^+^g = supL'"'g. 

\ n / m m 

□ 

Later on, it will turn out that if D is based on either cubature or piecewise harmonic 
interpolation or the reduite, the perpetual Hmit is, in fact the minimal fixed point (cf 
Theorem 15.11 and Lemma 15.61 for piecewise harmonic interpolation. Theorem 15.21 for 
reduite-based approximation, and Theorem l(i.ll for a result on a map D which is based 
on cubature). 

Moreover, we shall show that the algorithm based on cubature converges linearly 
in the sense of the definition above. 



1.2 Outline of this thesis 

We will postpone giving a more informal account of our motivation to use cubature 
formulae for Bermudan option pricing until we study the numerical implementation 
of some Bermudan pricing algorithms. In this thesis, we will first of all prepare the 
derivation of bounds on the natural scaling of the difference between an American and a 
Bermudan perpetual barrier option price (conceived of as a function of the Bermudan's 
exercise mesh size). Later on, we will show that this is sufficient to obtain bounds on 
the natural scaling for the difference of certain non-perpetual American and Bermudan 



10 



barrier options. For the one-dimensional setting, analogous results have been obtained 
by Broadie, Glasserman and Kou 0. 

Furthermore, we will prove soundness and existence of a perpetual limit for a 
number of Bermudan pricing algorithms, including pricing based on cubature. In 
particular, as was previously remarked, we will obtain a linear convergence rate for 
the latter class of algorithms. 

Later on, we shall prove convergence bounds for a non-perpetual American pric- 
ing algorithm in which one is computing non-perpetual Bermudan prices of a certain 
exercise mesh size via cubature and successively halves their exercise mesh size while 
leaving the number of paths at which the option is evaluated constant. 

The natural scaling for which bounds are derived in the first chapters of this report, 
can be used to consistently extrapolate from a finite number of Bermudan barrier prices 
to an approximation for the American price. En passant, we will sketchily explain 
how some of the features of object-oriented (C++) programming can be exploited to 
enhance the efficiency of a Bermudan pricing algorithm based on cubature (that is to 
say, how to circumvent exponential complexity by achieving recombination through 
the map class template). 

1.3 Notation 

We are following largely standard probabilistic notation, as can be found for instance 
in the works by Ito and McKean jr. QS] or Revuz and Yor |2()) . 

Both A C B and A C B for sets A and B will mean that ^ is a subset of B 
(possibly A — B). 
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Chapter 2 

Exercise regions 



In this Chapter, we will give a rigorous proof for the fact that an American/Bermudan 
option price coincides with the payoS that is expected if one exercises at the first 
possible entry of the log-price process into the immediate exercise region (G C M.'^). 
Characterisations of such regions for special cases have been proven in recent years 
[2] 



Definition 2.1. Given a countable subset 2 C [0, +00) and a Lebesgue-Borel measur- 
able set G C M'', often referred to as exercise region, we define the stopping time 

:= min{i el : Xt e G} , 

(the superscript will be dropped when no ambiguity can arise) which is just the first 
(nonnegative) entry time in X into G. If Q is a subset of space-time, that is G C 
K'' X [0, +00) rather than space (ie W^) itself, we use the space-time process rather 
than just the process itself to give an analogous definition: 

r| := min {teta+T : {Xt,t)eg}, 

where t^ is the time- coordinate at which the space-time process was started. Also, for 
h > we set 

h _ hfi h _ hN 

) 

to denote the first positive entry time in /iNq into G or Q , respectively, whilst finally 
fg :— T™° and Tq :— Tq^" will denotes the first nonnegative entry time into Q and G, 
respectively. 

For convenience, we will also adopt the following convention for this Chapter: 

Definition 2.2. Let I C [0,+oo). A stopping time t is called 2 -valued if the range 
of T, denoted by ran t, is a subset ofJU {+00}. 

Lemma 2.1. Consider a countable subset 2 C [0,+oo). Let X be a d-dimensional 
basket with an associated risk-neutral measure P and discount rate r > 0. Suppose 
g = {K — /) V 0, s > and e~^' f{X.) is a P-submartingale. For all 2-valued stopping 
times T there is a space-time region B — UmgxI"} ^ such that 



E 



for all T £ [0, +00) where 

T ■.= tI = inf {uel : Xue Bu} . 

7/ sup {e~'"^g (Xt-) : t stopping time} is ¥-integrable, then the latter inequality will 
also hold for T ~ +00. 

The Lemma holds in particular for X — sNq for arbitrary s > 0. 
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Proof. Firstly, we will treat the case of T < +00. Define 

\/teI Bt := Xr ({r = t}) C . 
Let us first of all assume that 

e X g{Xt) > a.s. on {Xt € Bt} 
and let us also for the moment suppose 

yteJ {f = t}n{T>T}^$. 



(2.1) 



(2.2) 



Both of these assumptions will be dropped at the end of the proof for the case T < +00 
in order to show the Lemma in its full strength. Now, from equations Ij2.2ll and (I2.1|l 
one may derive 



E 



e-^(^^^)g{X,^T),{r^t} 

E[e-^^g{Xr),{r = t}] 
E[e— (X-/(X,)),{f-t}] 

,-r(r^T) (^K-f{Xr^T)),{r^t} 



E 



(2.3) 



for all t e [0, T] n I. 

Furthermore, observe that r > f a.s. Using Doob's Optional Stopping Theorem 
(see eg Varadhan |^ Theorem 5.11]), we infer from our assumption of e~^'f{X.) 
being a P-submartingale with respect to the canonical filtration the assertion that 
(^g-r-vj {Xjj))^f^^f/\T tat} a P-submartingale with respect to the filtration {J-f/\T, ^tat}- 
Hence, if we combine this with equation H2.3II and note that {f = = {f A T = i} e 
Ttat for all t g [0, T) n 2, we obtain for every t e [0, T) n X, 



E 

= E 

< K -E 

< K -E 



e-^(^^^)g{XrAT),{r^t} 



-r(f AT) 



{r = t} 
{r = t} 



,-r(f AT) 



/ (X^^t)) ,{r = t} 



E 



"■(^'^^)/(X,^i.)),{f = 



= E 



'e-^^^^^^giXf^T)),{r^t} 



On the other hand, since r < r, if f > T, then also t > T, entailing 

fAr = T = TATon {f>r}. 
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Summarising these last two remarks, one concludes 

tein[o,T) 

tein[T,+oo) 
teJn[o,T) 

+ E IE 

tein[T,+oo) 



-r(fAT) 



g {Xf;,T) ,{r = t} 



E 



In order to complete the proof for the case of T < +00, let us show that the 
assumptions 112. l|l and Ij2.2|l are dispensable. 

If the assertion Ij2.1|l failed to hold, we would simply define the stopping time 

•= Xc(U,gj{r=t}n{g(XO>0}) • °° + E X{T=t}n{giXt)>0} ■ t 

and based on this definition, we would set 

ytel B[ ■= Xr' ({r' = t}) = Xr ({r' = t}) C M'*. 

Then we would have on the one hand H2.1II for B' instead of B which, according to 
what we have been able to show under the assumption of H2.1II . yields 

provided the condition H2.2II is satisfied. However, in any case 



E 



<^'^^)g{Xr,^T)] = E le-^'^^^^^giXr^T) 



Thus, with B' we have found a set that can play the role of B in the Lemma's statement, 
under the assumption Ij2.2|l . 

Now suppose the condition Ij2.2|l was not satisfied, and (12.111 may or may not hold 
(in the former case, one may even replace r' by r in what follows). In this situation 
we consider the stopping time 

^" = ^{r'<T}u({r'>T}n{r^,>T}) ' + ^C({r'<T}U ({r'>T}n{T5, >T}) ) ' °°- 

If one now defines 

ytel B't' := Xr" ({r" = t}) 

then 

ytein [0, T] B[' = Xr" ({r' = t" = t}) ^ Xr' ({r' = t}) = B[ 

and 

e X n (T, +00) B'l = Xr" ({r" = t}} = Xr" ({r" ^ t' ^ t} D {r|, > T}) 
C X.,({T' = t}) = B;, 
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hence 112. l|l holds for B" instead of B (since 112. l|l holds for B' instead of B and as we 
have just seen B'l C B[ for all t £2). Furthermore, 

'^B" = ^{T'<T}u{{r'>T}n{rl,>T}) ' '^B' + ^C({r '<T}U ({r'>T}n{r5, >T}) ) ' °° 
= X{r"<T} ■ T%, + X{+oo>T">T}n{rl,>T} ' '^B' 
+XC({r'<T}u({r'>T}n{r5,>T})) ' °° 

(the first line because of Xr" (C ({r' < T} U ({r' > T} n {t|, > T}))) = 0). There- 
fore 

tI„ >T on {t" > T} , 

thus 112.211 holds for t^„ instead of f and r" instead of r. But we have already seen 
that 112. l|l holds for B" instead of B. Therefore, using what we have proven under the 
assumption of both Ij2.1|l and Ij2.2ll , we get 

J L s . 

On the other hand, however. 



E 



(as t" = r' on {r < T}, as well as t" > T on {t > T}, thus r" A T = T = r A T on 
{r > T}) and we have already seen that 



E 



Finally, 



E 



'■(-'^^)5(X.,^^)1 = E \e-^^^^^^giXr^T) 



whence with i?" we have found a set that can play the role of B in the Lemma's 
statement. 

Finally, we need to consider the case where T — +co. The random variable 
snp{e~^'^g {Xr) ■ r stopping time} is an upper bound on e^^^ ^'^g {Xr'/\T) for all 
T > 0. Therefore, as soon as sup^ stopping time ^"'^^f (^t) is P-integrable, we are al- 
lowed to apply Lebesgue's Dominated Convergence Theorem and the assertion for 
T = +00 follows by letting T tend to infinity in 



E 



-r(TAT) / V \ 

e ^ 'g(XrAT) 



< E 



(which we have already proven for all T e (0, -|-oo)). 



□ 



Corollary 2.1 (Formula for an option price using hitting times). Let X be 

a d-dimensional basket with an associated risk-neutral measure P and discount rate 
r > 0. Consider a countable subset! C [0, -l-oo). Suppose g — {K — f)\J Q, s > Q and 
assume that the process e'^^' f{X.) is a V-submartingale. Then one has 

sup E \e-''^s g (x^, 

ecR''x[0,T] measurable '- ^ ' 

sup E[e-'^"5(X0] 

r Xr\[{).T\-valued stopping time 



for all T < +ao . If the random variable sup^ stopping time ' 
then the equation 



sup E 

6CR''x[0,+oo) measurable 

sup E [e" 

T X-valued stopping time 



^g (Xr) is V-integrable, 



X^ 



holds. 
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Definition 2.3. Let I C [0, +00) be countable, G C M.'^ x [0, +00) and G C R'^ 
measurable, and X a d- dimensional Markov basket with an associated family of risk- 
neutral measures ¥' and discount rate r > 0. We define 



^rtvrf(x,t) 



well 



G,X '■ ^ ' 



g{x), {x,t)eg. 



, x(^G, 
xeG. 



Instead ofVg^{x,0), we shall often simply write Vg^{x). Also, the subscript X will 
be dropped when no ambiguity can arise. Also, Vg and Vq will be shorthand for y^^" 
and Vq^° , respectively. 

As another notational convention, let us from now on use supg^j{dx[o,T] supg^^d 
to denote sup£;cR<ix[o,T] measurable and supccK^i measurable' respectively. 

Theorem 2.1 (Optimality of the immediate exercise region). Let X be a d- 

dimensional Feller basket with P' being an associated family of risk-neutral measures 
and r > being the discount rate belonging to P' . Suppose g — {K — /) V 0, X C 
[0, +00) is countable, and T G [0, +00]. Assume, moreover, that e^^'f{X.) is a P^- 
submartingale for all x . Define 



F- 



.I,T 



(x, t) eM"^ X [0, T] : sup V^{x, t) < g{x) 

ecR''x[0,T] 



if T < +CX) (we may drop the superscript T wherever this is unambiguous) and else 

pi, +00 _ 



{x,t) eR'^ X [0,00) : sup V^{x,t)<g{x) 

ecR''x[0,+oo) 



The 



V^,^rix,0) = U'''°{T){x) 



if T < +00, and Vp:r,+oa{x,0) = 11"^° {+oo){x) for all x e R'^ such that the random 
variable sup^ stopping Ume^^^'"' 9 i^r) is F^'-integrable. 

Proof. Let T < +00. Using Corollary |0] and recalling the definition of U^^° , all we 
have to show is 



VxeR'' V^^.t{x,0)= sup E(^'°) 

ecR'^x [0,T] 



9 if',) 



(where we recall that fg < Tg denotes the first nonnegative entry time into G). How- 
ever, after exploiting the special particular shape of F^ , we can - due to the bound- 
edness oi g < K which yields Vg < K for all Q which allows us to apply Lebesgue's 
Dominated Convergence Theorem - swap sup and E to get for all x €R'^, 



^/x(x,0) 



E(-,o) 



K^x) 



sup (x^T ij,^') 



sup e'"^FxEV 

ecR'^xfCT] 



sup E^^'*') 

ecR'^x [0,T] 



x„ >^pX 
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(where for notational convenience Tg should denote the first nonnegative entry time 
into Q). Now, let us use the strong Markov property of the Feller process X , and for this 



purpose, let 9 denote the shift operator on the space-time path space D ([0, +00) 



(which is the set of all cadlag functions from [0, +00) into W'' x [0, +00) - recall that 
all Feller processes have a cadlag modification). We obtain 



[0,+oo)) 



V^x(a;,0) = sup E^^-") 

acR'*x[0,T] 

sup E(^'°) 

0CK<*x[O,T] 

sup E(^'°) 

£(CK<*x[0,T] 



E 

]E(-^o,0) 



-rfSoS_x 



9 (^f J ° 



e 5 P^f 



But e '^'g{X.) is a P^-supcrmartingale for all x G W^, therefore by Doob's Optional 
Stopping Theorem, {e~^'"g {X^)) i -i.j,, ^x^^x ^j^a must also be a P^-submartingale 

for all X e M'^ and T' e (0, 00) (note that fg < fgof^j: a.s. because of the fact that 9 is 
the shift operator for the space-time process {t, Xt)t>o, rather than simply for X). Let- 
ting T' tend to infinity, we can employ Lebesgue's Dominated Convergence Theorem 
{us g < K yields e-™c/ {X^) < K e ^^(P^) for u e jfl A T',f^ o ^-x^ A T'j for all 

-rfgoS.x ( 

T e (0, -l-oo) and x e M ) in order to get that the expected value of e g X^x^a 

\ ^ 

is always greater or equal than the expectation of e^'"^i g (X^x ) . Hence 



-•rrgoti-x 

e ( 



V^x(a;,0) = sup E("'°) 

acK<*x[0,T] 

> sup E(^'°) 

5CK''x[0,T] 

sup E(^'°)l/g^(a;,0). 

ecK''x[0,T] 



pi 



The case T = +00 can be dealt with analogously. 



□ 



Lemma 2.2 (Time-stationarity of immediate exercise regions for perpetual 
Bermudans). Let X he a Levy basket with P' being an associated family of probability 
measures and discount rate r > 0. Then for all s > we have 

U^^°{+00){x) = ^ U-^o(+oo){x)<g{x)}^^^ 

for all X G Mf^ satisfying the condition that the random variable 
sup {e~'"^£f (X-r) : T stopping time} be ¥^ -integrable. 



Proof. Consider an integern e No, and an a; S M*^ such that sup {e ^'^g^Xr) : r stopping time} 
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is P^-integrable. Then we shift the time scale by ns to get 

rnsTffix.ns) 



sup I 

ecR''x[0,+oo) 



sup e 

ecK''x[ns,+oo) 



sup 

e'CR''x[0,+oo) 



]E(^'0) 



sup e '"^S' g (Xr = , ) 



where 6 denotes the shift operator on the space (as opposed to space-time) path space 



D{[0, 



Because of the boundedness oi g < K which entitles us to apply 



Lebegue's Dominated Convergence Theorem, we may swap sup and E to obtain 



sup ( 

ecR''x[0,+oo) 



sup e"''^E(^^"") 

gcR''x[0,+oo) 



sup yg 

ecR'*x[0,+oo) 



Vf°{x,ns) 



for all 71 G Nq. Thus we conclude 



sup yg 

acK''x[0,+oo) 



sup Vg 

ecK<*x[0,+oo) 



for all t e sNq. If we insert this equality fact into the definition of F'^^" , we see that 
the condition determining whether a pair (x, t) belongs to F^^° does not depend on t. 
On the other hand, by Corollary 12 .11 

sup Vf°{x,0) = W'''>{+oo)ix), 

ecR''x[0,+oo) 

and the left hand side equals - by our previous observations in this proof - the term 
featuring in the definition of F^^^" . 

□ 

Summarising the two previous Lemmas and applying them to a more concrete 
setting, we deduce that the expected payoff of a perpetual Bermudan option of mesh 
size h > equals 



gix), x€G ' 



C/™«(+oo)(.t) = 
where G = {t/™°(+oo)(-) < g(-)}. 

Lemma 2.3. Let us fix a Levy basket with an associated family of risk-neutral prob- 
ability measures P' and discount rate r > 0, as well as a region G C M.'^ and a real 
number h > 0. Then we have 

Vx^G V^{x) = e--''Px„-x, * V^{x). 

In particular, using Lemma \2.2l one has the following equation for the expected per- 
petual Bermudan option payoff: 



-rhv 



■ Xq — X 



,,U'^''«{+^)ix). 
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Proof. Using the Markov property of X, denoting by the shift operator on the path 
space £> ([0, +00), M'') of a Levy process X, and taking into account the fact that 
Tq > (i.e. Tq > h) in case x ^ G, we obtain: 



^rg ° f'lj • X{r^>h} 



= e-^^Fxo-x,*V^{x). 





= e~ 


-rh^x^x 




r{r^-h 




= e~ 


rh^x^x 


e' 






= e~ 


rh^x-^x 


er 






= e~ 




e~ 





□ 
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Part II 

Bounds on the 
American-Bermudan barrier 
option price difference 
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Chapter 3 

Scaling the difference between 
perpetual American and 
Bermudan barrier options 



Embracing the terminology of Broadie, Glasserman and Kou 0, we shall refer to the 
difference between an American and the corresponding Bermudan options (on the same 
basket and with the same payoff function) as "continuity correction". 

3.1 The exercise boundary and its relevance for con- 
tinuity corrections 

Lemma 3.1. Consider a d-dimensional Feller basket X with an associated family 
of risk-neutral probability measures P' and discount rate r > 0. Furthermore, let 
f -.W^ ^ R4. be a nonnegative continuous function such that e^^' f{X.) is a martingale, 
K a nonnegative real number, and define g := K — f . Finally, let G be a measurable 
set such that g > on G. Then we have for all x e W'' and s > 0, 



V^{x)=K-Y^. 



— r7lSTn)X 



fl {x,s ^ G} n {x^s e G} 



Similarly, if g :— f — K instead and one assumes that this g is nonnegative on G, the 
identity 



V^{x)^f{x)~K-Y^e 



n=l 



fl {x,s ^ G} n {x^s e G} 



holds for all x E M.'^ and s > 0. 



Proof. Let a; e M'' and s > 0, and let us first set consider the case oi g = K — f . Since 
g > on G by assumption, one has the identity 



= W 



■■{K-f) (X^s 



Moreover, 



lim e-'-(^SA-)(^ _ f) (X,.^„,) = e-^-o^K - f) (X,. ) 



<K 
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i"^-almost surely, therefore by Lebesgue's Dominated Convergence Theorem, 



VSix) = Hm 



{K - f) 



(3.1) 



But since e^*"/ (X) = (e"''*/ (Xt))j>Q is a P'^-martingale, we may apply Doob's Op- 
tional Stopping Theorem to get that (^e"''^'^^^"'*) J (X^g/\„s))j^gj^ is a -martingale, 
too, whence 



Vn £ No 



-r{TQAns) 



fix). 



This finally yields, because of equation H3.1II and the monotonicity of the sequence 



which completes the proof as 



K ■ lim E^ 

n — >oo 



-r(r^ Ans) 



fix) = K^^ 



fix) 



E^ 



-rnsmx 



fl {X,, ^ G} n {X„, G G} 



for all n e No- 

The case oi g — f ~ K can be treated analogously. 



□ 



An important feature of the immediate exercise region for a perpetual Bermudan 
option with payoff function of the form (i^T— /)V0 where / is monotonely increasing in 
each component, and exercise mesh size h, is that - owing to the fact that J7''(-|-oo)(-), 
the option price as function of the logarithmic start price, is monotonely decreasing in 
each component - it is south-west connected in the following sense: 

Definition 3.1. A set E C Mf^ is called north-east connected if and only if for all 
a; e K'* such that x > componentwise, E + x Q E. Likewise, any set F CM.'^ is called 
south-west connected if and only if for all cc S M'' such that x > componentwise, 
F + x^ F. 



'xeF ■^i)i=l 



Remark 3.1. IfF is a south-west connected subset ofMf^, thensnpF := (sup, 
is an element of the boundary of F. Analogously, if E is a north-east connected subset 
o/R'', then inf i? := {infxeE Xi)f^-^^ is an element of the boundary of E. 

Lemma 3.2 (Characterisation of the American-Bermudan barrier difference 
for perpetual puts). Let X be the logarithmic price process of the multidimensional 
Black-Scholes model with constant volatility and interest rate, that is 



yt>0 Xt=l (Xo), + C7, ■ {Bt)^ 



2 



(where B is the d-dimensional Wiener process) for some r > and a G M>o • Let 
g = K — f , wherein K > be a real number and f > be a continuous function 
that is monotonely increasing in each component and such that (e~''*/(-'^t))t>o ^ 
martingale. Finally, consider a measurable set of the shape G = j — H for some 
7 G M'* and some convex north-east connected set H C ]R_|_'^ (making G convex and 
south-west connected) such that g = K — f is nonnegative on G. Then we have for all 
s > 0, 



G 



CG. 
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Proof. Fix an s > 0. Introduce stopping times Tq ^ for n € N, t > through 

Vn e N > r^_„ := A ns 

and define 



Vn e N Vf > := E' 



on 



Now consider an arbitrary n e N and t > 0. Due to Doob's Optional Stopping 
Theorem, applied to the two-component sequences of stopping times ^0, Tq „^ and 

(o,Tq „^ (which both are bounded by ns), combined with the fact that e~^' f{X.) is 
a martingale, 



thus asg = ii' — / >OonG and Tq^ < Tq, 



= E 

+E 
+E' 
= E' 
+ E' 



e— In - /) (X 3/.) - e— {K - /) J , {r^ < ns} 

e—l.. (X - /) (x^,/, ) V - e-'-^S.n (jc - /) (x,.. J V , [rf ,t^g > ns} 
e— l.n - f) (X^.n^ V - e— [K - /) (x,. J V , [rf <ns< r^} 

e---* {K - f) {X,) V 0|;:;J;; , {r| „, > ns} 



=0 



+E 
E 
+ E' 



e-"-* {K - /) (X,) V 0|;:;f" , {tJ,„ < n^ < raj 
e-'-* (i^ - /) , {r* „ < ra < ns} 

e-*(if-/)(XO|;:;f;;;:::,{r|„,ra>ns 



+E' 
-E 
E 

-E 



=0 



e-'-* (i^ - /) {Xt)\lZll2 , < ns < r^} 

e-™« (if - /) (X„,) A , {r|_„ < ns < r^}] 



(X-/) (X„,)AO,{r*„<ns<TS} 



Note, however that the subtractor converges to zero exponentially. For, exploiting 
not only the Markov property, but also / > and that e~'^'^^-^'^ is a martingale, we 
have 
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> £■ 

= £■ 

+] 
= E' 



)s < ns < 



-K 



^ {K-f) (X„,)AO,{rJ„<ns<Ta} 
^K-f) (X„,) A 0, {r| „ < (n - 1)^ 

(K - f) (X„,) A 0, |t|^„ 

[e-™^ {K - /) (X„,) A 0| {X(„_i), ^ G}] F {t| „ < (n - 1) s < < r^} 



— (i^-/)(X„,)AO|{X(„_.),eG} 
> E- [e-™^ (if - /) (X„,) A 0| {X(„_i), ^ G}] F 



,^ = I n — - I s < ns < T, 



{X(„_i), ^ G} n fl {X,, ^ G} 

n-1 



-E' 



(if - /) (X„,) A o| {X(„_ .),,^}] F {X(„_ e g} n fl {X,, ^ G} 



> 



IE- [e-™^ (-/) (X„,)|{X(„_i), ^G}] 



{X(„_i), ^ G} n fl {X,, ^ G} 



1=1 

ri-1 



-E' 



(-/) (x„,)| {X(„_.)^,^}] F G g} n fl {X,, ^ G} 



i=l 



> E- [e-"- (-/) (X„,) , {X(„_i), ^ G}] 



(-/) (Xn.),{X(„_.),eG} 



> E- [e-™^(-/)(X„,)]F 
+E- [e-™^ (-/) (X„,)] 

= E-[(-/)(Xo)]P 
+E-[(-/) (Xo)] 



i=l 
n-1 



1=1 

n-l 

'i=i 

{X(n-l)s ^ G} 



{-'^(n-l)s ^ G} 

{^(„-.). G g} 



{^(„-.). G g} 



=1 

n-l 



{-''^(n-l)s ^ G} 



-/(.) F 



{-'^(n-l)s ^ G} 



{^(„-.). ^ g} 



on CG - and both of the conditioned probabiHties in the last Hne converge to nought 
n tends to +oo. 

Therefore, due to equation H3.2II . 



lim E' 

71 — >OQ 



lim K ■ E' 



(3.2) 
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on CG. 

However, we also know that 



■(K-f) {X,)\ 



e-^' (K-f) {Xt)\ 



as 71 ^ oo 



almost surely. Combining this convergence assertion with Lebesgue's Dominated Con- 
vergence Theorem yields 



vj-vs 



lim 



G,n \ I 



and recalling equation Ij3.2ll . we finally arrive at 



CG. 



(3.3) 
□ 



Next we shall prove that the expression on the right hand side of the last equation 
is monotonely increasing as • J, 7 componentwise. 

Lemma 3.3. Under the assumptions of Lemmc ^S^ one has for all s > and x,y E B."^ 
that satisfy the relation 

X > y > J componentwise 

the following lower bound: 



00 f 



4=0 

1 



(3+«o)' 



In the proof of Lemma l331 we will apply the following 

Auxiliary Lemma 3.1. Consider a measurable set A CM. of positive Lebesgue mea- 
sure and two continuous positive integrable functions p,q E {A) on A, as well as 
a nonnegative monotonely decreasing function f £ {p dX) C\ {q dX), and assume 
that not only is the function | monotonely decreasing, but also jj^p dX = J^q dX. 
Then 

' fpdX> f fq dX. 

J A 

Proof. If the continuous monotonely decreasing function p — q was either > or < 
on all of A, then one would get J^p dX > J^q dX or J^p dX < J^q dX, respectively. 
Hence there exists a real number b E A such that p{a) = q{a) as well as p > g on 
A C (—00, a) and p < q on A C (a, +00). But since / is monotonely decreasing and 
nonnegative, these inequalities yield 



An(— 00,0) 



f-{p~q)dX + 



yln(a,+oo) 



f-ip-q)dX 



> 



An{ — oo,a) 



f{a){p-q)dX 



An(a,+oo) 



f{a){p^q)dX>0, 



thus 



f fpdX~ f fqdX= f f-{p- q)dX > 0. 
JA J A J A 



□ 
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Proof of Lemma \'-i.'A Let x > y > which in particular entails x > y > supG. Let 
us define, for all z GMf^ and z G No, a measure Q(^'*) on G by 



Q 



(z,i) _ I'+^j" 



Then, due to our choice of X as being a Brownian motion with drift, Q^^^*^ will have 

a positive continuous Lebesgue density, denoted by ^^^ja— • 
Now whenever zq,zi £ G such that 



3q > 1 (sup G ~ zo) — a ■ (sup G — zi) > 
(in particular zq < zi), note the identity 



(3.4) 



%diag (^i^...,<T,2)(^+i)s (• - ^ - + 5) 



(where ga,c shall for every symmetric positive semidefinite G G R''^'' and a G M'' 
denote the density of h'a.c, the Gaussian measure with covariance matrix G and mean 
a) for all z eW^. For, we shall then obtain 



|-o---(2+iH! 



^ Nol^-|^lP-2-(i+^)s-'M(^0-»l) 



-|;==ll^-2'((.+ l)s^. + ;,)(^0-^l) 



and therefore 



This implies 



> 



dA 



— (2:0) 



dQO^/ X - dQ(^/ X- 

dA'i y^^) dX'i v^ij 



K^o) < 37v;7-t(^i) 



(3.5) 



whence we have shown that ^{y,i) is monotonely increasing on the ray-segment 
R{z2) := (^supG - Z2 ■ IR+''^ n G. But in addition. 



(i+«o) = 





(i+«o)= 









will hold (since G was assumed to be convex and south-west connected). Thus 

(i+No)= 

E' e^'"^G decreases monotonely on any ray (or ray-segment) R{z2) = (sup G — Z2 

G (for Z2 > componentwise) if we look at this ray (-segment) as a linearly ordered 
set with respect to the componentwise real order relation <. From these two sets of 
monotonicity assertions (for all Z2 > componentwise), together with the fact that 
Q(^'^)(_G) = Q^y'''^\G) = 1, one deduces by means of the above Auxiliary Lemma, l3?Tl 





(i+«o)/ 







)(^'^)(dz) < 





(i+No) = 




'r{z2) 
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therefore 



< 



r 


r 


V 2 ^ / 




/ 




/aBi(o)nE+^ 


'r{z2) 








(i+i>io) = 


/ 


[ E'- 




/OSi(o)nE+^ 







{dz)- 



dZ2 

9Bi(0)nM+'^ 
(^2:2 



A 



asi(o) n: 



(where i?i(0) denotes the d-dimensional unit ball) which via the Fubini-ToneUi Theo- 
rem amounts to 



(i+No) = 



(3.6) 



On the other hand, the definition of Q^"^'^^ gives 

Q^'"^''\dz) = 



(i+No)< 



{^(^+^)^ ^ ^} 



for arbitrary w €E M'^. Hence inequality (j3.(i|l becomes 



E^ 



< E^ 



(3.7) 



The Markov property of X, moreover, entitles us to state 



E" 



(l+»o) = 



for all i and w € M''. Thus we can use inequality 113. 7|l to derive the following estimate: 



E^ 



1=0 



(3.8) 



1=0 



1 



• 1 - E' 



(i+N„)s 



(3.9) 



^ El 

i=0 



z + - I .s > e-K'+i)'^ 



1 -E"= 



(5+«o) = 



which via Lemma l3?2l gives the desired result. 



□ 



Corollary 3.1. Suppose the assumptions of the preceding Lemma VTl^ hold. Consider 
X > y > J and assume, in addition, there is a positive lower bound for the sequence 

V«{r|=(^+i).| 



Then one has 



{v3-VS)iy)>{M 



H = (»+i)4 ' 
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Proof. Equation and Lemma lil 21 yield 



oo 



. 2 

G ~ 



« + - I s > e 



r(»+i)s 



• 1 - 



(3.10) 



for every z ^ G. Inserting x and y for z in this estimate, the Corollary can be deduced 
via the estimate of the preceding Lemma 

□ 

We can explicitly state a partial differential equation that the said difference Vq — 
Vq obeys: 

Lemma 3.4. Consider a measurable set G C . Let again X he the logarithmic price 
process of the multidimensional Black-Scholes model, that is, 



( 



yt>o X, = 



\ 



V 



and let us assume that g and f are as in Lemma XS^ Then for allt > s > 0, the partial 
equation 



11=0 



holds on £iG, wherein for all u, v > 
and (always following the notation introduced in Chapter^ 

Proof. Fix s > 0. We shall prove the Lemma by studying the space-time Markov 
process {Xt, s — i)tg[o,s] and the functions 



: CG 



(x, t) h— > 



where 



yt e [0, s) T^{t) := T^-* = inf {m e sN - i ; X„ e G} . 

It is clear that 

yU>te[Q,s) Tq{u) O 9^_t ^ TQ{t) - {U- t), O e'„_f = X^a(t), 

hence if we employ the Markov property in two directions we can for all u > t get the 
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following P' -almost sure identities: 

(X.. („)) 

= E-^'E-^" 
= e'^("-*)E-^' 



g-rrJ(«)oe„_t^ 



u- 



— r(u—t) —rT%. (t) 



Thus, the process {f^ {Xt, s — i))jg[o s] ^ martingale. The infinitesimal generator of 
the space-time Markov process {Xt, s — t)tg[o,s] is 



d_ 

dt ' 2 



^ a 1 " 

- '^^^dA - VV = -dd+1 + '^^^dA - 



due to the well-known result on the infinitesimal generator of Brownian motion with 
drift (see eg Revuz and Yor ^26| p. 352]). Therefore we have proven 



vs>o o = Lr = -^r 



Now one observes that 



which yields 



Vs > 



Vs>0 /"(^O) = V(5 on CG 



du 



and thus brings the proof of the Lemma to a close. 



□ 



Remark 3.2. In the remainder of this Chapter, we will continue to largely focus 
on put options, thus always setting g = K — f for some componentwise monotonely 
increasing f and assuming the G occurring in the defition of V^, s > 0, to be south- 
west connected. However, one can easily derive analogous results for call options on 
dividend-paying assets, by simply cutting the interest rate to discount the dividends 
and by replacing g by —g and G by some north-east connected measurable subset ofM.'^ 
that is assumed to satisfy the condition f — K > on G. 



3.2 Continuity corrections in a one-dimensional set- 
ting 

Definition 3.2. For measurable G CM.'^ and s > define 
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Remark 3.3. Note that in general, VA ^ g ^ Vq on G, hut 



= on CG. 



Theorem 3.1. Suppose d ~ 1, let G — {~-co,j] and g ^ K — exp, and assume 
{Xt)t>[) = {^Xq + a ■ Bt + {r — ^)) ' words: X is the logarithmic price process 
of the one- dimensional Black-Scholes model with constant volatility a and discount 

2 

rate r > Q. Set fi :^ r ~ ^ . Then one has for all s > the relations 



Hmt/^(7)-F^(7) 



K ■ 



exp 



E 



fl {x,t > 7} n {Xnt < 7} 
° {Xnt < 0} 



.i=l 



(3.11) 



Furthermore, > 0, there exist constants co,Cq > such that for all sufficiently 
small s > 0, 

cqs^ < limV^h) - VAh) < Cqs^ 
tio 

2 

// both /i < and r > there exist constants ci, Ci > such that for all sufficiently 
small s > 0, 

cisTJ < 111111/(^(7) - Vg(7) < Cis*. 

Remark 3.4. Although computing the constants cq, Cq, ci, Ci explicitly is possible, we 
refrain from it for the moment, as it is not required to find the right scaling for an 
extrapolation for Vq from s > to s — and it would not provide any additional useful 
information for our extrapolation purposes. The same remark applies to all examples 
and generalisations that are studied subsequently. 

Proof. The existence of linitjo ^g(7) ^ consequence of Lemma [r!2l The first identity 
in the statement of the Theorem is a consequence of the previously estabhshed Lemma 
13.11 whereas the second equation in the statement of the Theorem follows from a 
result by Feher jlll p. 606, Lemma 3] on processes with stationary and independent 
increments. For, if we define 

00 fn— 1 

Vs > e [0, 1) ^q, s) ■■= E n > 0} ^ ^ 0} 

n— 1 Li— 1 

then Feller's identity ^l] p. 606, Lemma 3] reads 

Vs>0 yqe[0,l) - In (1 - s)) = E — ^ 0} (3.12) 



n=l 



and holds whenever X has stationary and independent increments, in particular for all 
Levy processes (note that our definition of a Levy process requires them to be Feller 
processes in addition). This entails 



- E 

n=l 

= 1 — exp 



rnsmO 



fl {x,s > 0} n {Xns < 0} 

i = l 

^ p—rns \ 



(3.13) 
(3.14) 



Tl = l 
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which is enough to prove the second identity Ij^.llll in the Theorem. This ushers in 
the derivation of the estimates on exp ^— '}2^=i — {^ns < 0}^ which are needed 
in order to prove the inequahties of the second half of the Theorem. We shall show 
that if /i > 0, there exist constants co, Cq > such that for all sufficiently small s > 0, 

(oo — rns \ 
2 

and if both /i < and r > there exist constants ci,Ci > such that for all 
sufficiently small s > 0, 

Now, the scaling invariance of Brownian motion yields for all n G N and s > 0: 

PnX„. <0} = P"{b„, <-^ns}=P"{i?i<-^Mi/2} 

, /-^("^)'^' /-x^\ 
= (27r)-i/2 / exp — dx. (3.15) 



We divide the remainder of the proof, which will essentially consist in finding estimates 
for the right hand side of the last equation, into two parts according to the sign of fi. 
Case I: n >0. In this case we use the estimates 

Va;<OVy<0 -y -x < - ^ Y' 

thus 

Vy < 

e-y' f e~^'dx < f exp f-^^) dx < e^^' f e'^'/'dx, 
hence by transformation for all ?; < 

/tT _,2 fy ( x'^\ /TT 



^-e-y<J expl^--jdx<^-e--. (3.16) 

Due to equation II3.15|I . this entails for all n G N, s > 0, /i > (if we insert — f (fis)^/^ 
for y) 

e V CT ; g 2cr2 

^ < P" {Xns < 0} < . 

Therefore for arbitrary r, s > 0, 

^ Ti— 1 n— 1 n— 1 

The sums in equation (j3.17ll have got the shape of X] 9"/"- for g < 1. Now one performs 
a standard elementary computation on this power series: 

EtT - E / T.'''dr^ dr^-Hl-ql (3.18) 
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which immediately gives 



1/2 

1 - e"H'-+f 



< exp - ^ {Xns < 0} 



< 1 - e 



n=l 



2,/2 



when applied to equation lj^H7|l . Due to de I'Hospital's rule, the differences in the 
brackets on the left and right hand sides of the last estimate behave like s when s | 0. 
This is sufHcient to prove the estimate in the Theorem for the case of ^ > 0. 

2 

Case II: ^ < and r > ^ . In that case we employ the estimates 

Vx<o Vy<o _-l_vl^J^!i^^_-l^yl 

^ - 2 2- 2 - A 2 



and proceed analogously to Case I, to obtain 

< J exp ^-y ) da; < ^Tre^. (3.19) 

In the special case of y :— ^(ns)^/^ < 0, this leads to the esimate in the statement of 
the Theorem via 

oo —ns(r+-^^^ oo —m^i -t oo —ns(r — 

n— 1 n—1 n—1 



Therefore in case fj, = the scaling exponent is exactly 5. □ 

Corollary 3.2. Assume d ^ 1 and let, as in the previous Theorem JL {Xt)t>o = 
r — ; in words: X be the logarithmic price process of the one- 

dimensional Black- Scholes model with constant volatility a and discount rate r > 0. 
Furthermore, suppose g = K — exp and let G" := (—00,7''] denote the optimal exercise 

region for a (one- dimensional) perpetual Bermudan put option of exercise mesh size 

2 

s and strike price K on the (one- dimensional) basket X . Define fi := r — ^ , jq 
suPs>o7'' '^'^'^ •= (^o°j7°]- Then we have for all s > 

Hm V^o (7°) - VS^ irl = 0(7° ~Y)+K- 

Moreover, there are constants co, Co, ci, Ci, such that «/ /i > 0, 
Co-s^2> Hm 17,^0 (7°) - VS.iY) >co-sK 
2 

and if both fJ, < and r > , 

Ci-s^> limV'*o(7*') - V'g47') > ci • s^f 
for all sufficiently small s > 0. 



\ n=l / . t\n 
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Proof. The first asymptotic identity in the statement of the Corollary follows from 
equation Ij3.11ll in Theorem 13.11 as soon as we have remarked that 



and 



limV-^„(70)=liml-^o(7°)=5(7°) 



limF^.(7^)=.limy^.(7^)=.g(7^), 



for these equations yield 

limy^o(7°) - n^iY) = g(7°) - gii^+limVUin " V^^ii^ 
which by the differentiability of g means 

We can now use results on the exercise boundary for perpetual Bermudan options 
obtained by Boyarchenko and Levendorskii equation (5.3)] who showed 

7^ - 7" ^ si 

for sufficiently small s, and the estimates in the Corollary follow directly from the 
estimates of Theorem 13. II □ 

Remark 3.5. Up to this point, we have derived estimates for the American- Bermundan 
option price difference at the boundary 7 of the exercise region G = (—00, 7] (in case of 
a put) or G = [7, -\-oo) (in case of a call with dividends) . We can extend these hounds 
of the American- Bermudan difference from the exercise boundary to the complement 
of the exercise region: By continuity, we can even extend the lower bounds or upper 
bounds, respectively, to a neighbourhood of the exercise boundary: For, if we consider 
a put for the moment, we get from Lemma \3.4\ that {x, s) 1-^ Vq{x) is continuous for 
all G, implying that if x ^ G = (—00,7] 

^.^ ln(Vg(7)~-Fj(7)) ^ ^.^ J. \n{VS{x)-VS{x)) ^ 
sio In s xIj sio In s ' 

where Vq is shorthand for IhntioVQ. Thus, if limsjo '"(^'^'■'j'^ a £ (0,1) - 

where we have, thanks to Theorem AS . 1\ estim.ates for the limit lirngjo (^^(t'^^^gC-)')) _ 
we will for any e > get a S > such that 

Vx.[7,7 + ^) u^1^(^c(^)-^c(-)), + 
sio Ins 

Analogously, we can proceed to derive bounds for the American- Bermundan call option 
price difference (for an option on a dividend-paying asset) in a neighbourhood of the 
exercise boundary. 

3.3 One-dimensional continuity corrections outside the 
Black-Scholes model 

The identity 113.1111 of Theorem l3.1l ca,n be used to derive estimates in the spirit of the 
second half of Theorem l3.1l in more general situations. We will illustrate this by means 
of the following example: 
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Example 3.1 (Merton's jump-diffusion model with positive jumps and "mod- 
erate" volatility). Suppose the logarithmic price process X is governed by an equation 
of the form 

Vt > Xt^ Xa+at + pZt + aBt 

where a G M, /?, cr > 0, Z is the Poisson process (thus, in this setting, only positive 
jumps are allowed for simplicity) and B a normalised one- dimensional Brownian mo- 
tion, and the stochastic processes B and Z are assumed to he independent. Let P' he 
an associated family of risk-neutral measures and r > the discount rate. In order to 
employ we shall compute the sum Yl^=o — ■'^^ i-^ns < 0} for all s > 0. Since 

{Xns < 0} = P° 1^^^ — ^} arhitrary n, s we may without loss of generality take 
P = 1. Let us also assume a > 0; note that since exp (Xt ~ i"t)t>o martingale - as 
X is a logarithmic price process - , a must he such that a — r + ^ + [e^^] = (if 
r > and a > are given), hence a > implies a < y^2{r + EP [e^i]). Now, hy def- 
inition of the Poisson distrihution together with the symmetry and scaling invariance 
of Brownian motion 

OO — VTLS 

V - — P" {x^s < 0} 

^-^ n 

OO OO / \ ^ , TTiS 

= y.y.^ V°{aB^s<-ans-k} 

^-^ fc! n 

n=0 fe=0 

n=0 fc=0 ■ ^ 



(with the convention that 0° = 1/ Now let us first of all try and find estimates for the 
prohahility in the last line. By equation f.'^.lfip applied to y :— —%{ns)^ — ^{ns)^^ < 0, 



< 



< 



n ( „ OL , .1 k , ,l1 

^ {Bi< (7is)2 - -{ns)-^ \ 



£(ns)2+i(„s) 2 



e 



2 

which yields, using the ahhreviation a' ^, 



e 



< P°<Bi< (ns)2 (ns) 

a a 



< 



1-2 

e 
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' (3.22) 



2V2 

< P"<^Sl< (71S)2 - -(ns)-2 I 



< , 



Thus, we can perform the following estimates to derive an upper hound of the sum in 



00 00 



kl n [ ~ CF a 

n=0 k=0 ^ 



^ — Efcii^^---' 

00 

E 

00 



^ 00 ns(^l+r+Vj 

2 ^ n 

n=0 



00 — ns| l+r+^=:, e 

1 \ ^ e 



n=0 



where the last line uses that ^ = ^ > a • (r + E° [e^^]) > and we need to impose 

the condition that e~°''('^~^^ < 1 + r + ^3 (which, given r > and a, will be 

satisfied if a > is sufficiently small) to employ the identity 

°° „n 1 

Vq < 1 V — = In . (3.23) 

11=0 ^ 

The lower bound follows simply from 

°° 1 ' k ^ 

VneNoVs>0 '^-^(ns-e~^'^~^^ >1 

fe=0 

(forn — recall that 0^ — 1 in this paragraph by our earlier convention) as this entails 
(when exploiting the estimate \S.2^) and finally \S. 2Sp ): 

1 

2Vl: 



n=0 k=0 

1 - g-n.(H-r+a'^) - ^ 

^ — E — ^ — Efcii---^- 



> 



n=0 fc=0 



2\/2 



Lin (1-6-^1+'^+"") 

2V2 V 



35 



As a consequence of these estimates and using the Taylor expansion of exp around 0, 
we now get the existence of two constants C3 > and C3 > (which can be computed 
explicitly) such that for all sufficiently small s, 

- P° < 0} < C3 • S^. 

n=0 J 

Finally, we may apply identity f,9.j j|) from Theorem - as this is an immediate 
consequence of Feller's identity '11.^ p. 606, Lemma 3] and our Lemma \S. 1\ - and 
conclude that if g = K — f and G = (—00,7], then 

Kc3 ■ s5 < UmV^clT) - VSi-y) < KC3 ■ s"^ 
for all sufficiently small s > 0. 



3.4 Continuity corrections in higher dimensions 



The proof of Theorem l3.1l relies heavily on the use of Feller's result Jll p. 606, Lemma 
3] which in turn is proven by means of elementary Fourier analysis and a so-called "basic 
identity" [TT| P- 600, equation (1.9)]. 

Hence, if one aims at generalising Theorem 13.11 to higher dimensions, one should 
first of all find a multi-dimensional analogue of the said basic identity. 

Indeed, we shall see that this is feasible. Let us for the following fix a stochastic 
process X = {Xt)t>o on R"* with stationary and independent increments. 

Lemma 3.5. Suppose H is a measurable subset ofM."^, and s > 0. Define for all n E N 



yKeB{R'^) Rn{K):= 



fl {x,s e Ci/} n {x^s eKHH} 



l<i<n 



as well as 



fl {X,, e Ch} n {Xn, eCnnK} 



l<i<n 



(in particular Rq = 5q [■ O H] = and Qa = (5o [• H Cff] = Sq ). Then for all ri e No, 

Q„+i + i?„+i = g„*p%,. 

Proof. Consider a measurable if C M''. Clearly, 



{Qn+l+Rn+l){K) 



fl {X,, e Zh] n {X(„+i), e k] 



(3.24) 



On the other hand, since X is a Markov process, we have 

Qn{K)^{Ps{x^H-)rXK{Q) 

(where {Pt)t>n {^Xt * ')t>o translation-invariant Markov semigroup of tran- 

sition functions for the process X whose increments are stationary and independent), 
thus 

' f{y)Qn{dy)^iPs{xcH-)rfio) 
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for all nonnegative measurable functions /. But this implies 

{Qn*KJiK) = f f XK{z + y)P''xSdz)Qn{dy) 



= (Ps ixCH-)r XK-{z)P'xM^)j (0) 

= {Ps{XCH-)roP^XK{0), 

and the right hand side of this equation coincides with the one of identity lj^i.24|l . □ 

Applying Fourier transforms we obtain 

Corollary 3.3. Let us adopt the notation of the preceding Lemma and define the 
Fourier transform of a countable sequence {^'n)^ of finite measures on , denoted by 
(;Z0; = Ai: (0,l)xM'^^C, by 

oo „ cso 

yq e (0, 1) VC e C) = E / = V ^"^^(C)- 

n=0 -^R" n=0 

Then for all q E (0, 1), and ( E R'^ the equation 

1 - (RXii^ = c) (i - (o) 

holds. 

Proof. The result of the previous Lemma reads 

when we apply the Fourier transform. After multiplication with g""*"^ and summing 
up over rt € No, one arrives at 

Vge(0,l) VCeR^ C) - XChW +%, C) - XhW = 9^(9, C)l«^(^^ 

=Qo(C) =So(C) 

hence 

yq e (0, 1) R{q, •) - 1 = qQ{q, O^VH " Qil, O- 
This is our claim. □ 

Definition 3.3. A subset A C M."^ is called +-closed if and only if A is measurable 
and A + A C A, that is sums of elements of A are again elements of A. 

Lemma 3.6 (a la Feller, Wiener, Hopf). Suppose H is a +-closed set and its 
complement Ci? is a +- closed set as well. Assume furthermore ^ H 
Ro = 0), and let In the main branch of the logarithm on C. Then 



n=l ''^ 



(dx) 



for all {q, Q E (0, 1) x 'R'^ such that the left-hand side is well-defined. In general, for 
all q E (0, 1), one has at least 



1 - R{q, 0) = exp 
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Proof. Let q £ (0, 1). According to the previous Corollarv KlHI we have 

VCeC/ In =hi i^-hiQ(9,C) (3.25) 

wherever this is defined. Due to the identities IT ^ ^-"^ r e i?i(0) C C 

(cf equation Ij3.18|l in the proof of Theorem 13. II above') and P^x^ = P°Xs*" this can 
also be written as 

/ e* *C.(pO^J-(rf^) (3.26) 



oo - oo / ^ \n, 

E^(%.c))%E^(q('^'C)-i 



n 

n—l n—1 



However, at least for C = and arbitrary choice of q, one may still state identity II3.25|I 
as this follows from Corolla,rv 13.31 more or less directly: First we note that 



ln((l-gP^(g,0)) •(aX('Z,0)) 
In (oXiq, 0) - In (l - gP^ (9, 0) 



(as in these statements the arguments of In are positive, hence surely in the domain 
of In) and written in series notation 

E-(l-W«)n('?'0)) +E^P"x.(g,0)" 



n 

n—l 



E ^ ((i - '^1*^(9' 0)) • (oX (9,0) - 1 



n=l 



But Corollarv 13.31 implies 

VneN l-<7P0x.(g,0) •(Q„)„(9,0)-l = -(i?„) Jg, 0). 

Combining these two equations yields H3.25|l . Next, note that 



^J■R,q 



is still a finite measure - concentrated on H - and thus possesses a Fourier transform. 
Analogously, the measure /j,Q,g :— X]^o I^Qn is concentrated on Cff and also has a 
Fourier transform as it is finite. Now, for arbitrary n € N, the properties of the Fourier 
transform imply 

{Q{q,-)-lj = yf^Q.q-6oj = {{flQ.q- So)*"). 

But since H and CH are +-closed sets, i.e. H + H C H and CH + CH C CH , the 
measures on the right hand sides of these two equations, ^iR,q*^ and (/J.Q,g — 5o)*"; 
have to be (signed) measures on H and CiJ, respectively. Let us now split the sum in 
l|3.26|l and insert the terms we have previously identified: 

VC e i7 

V^/ e'-'C.(pO^J-(rf^) + V^ / e^-'':-(P%J*"(dx) 

n=l •'H „=i " 

OO 



E -(a^«.«*")(0 + E ^(MQ.. - ^o)*"(C). (3.27) 



n=l 
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It is the injectivity of the Fourier transform that yields from this 



oo 



-((r°x. 

n \ 



i-nH) 



oo 



^((P%J™(.nCi7)) 



(-1)' 



{PQ,q - So)*" 



Either side of this equation equals the sum of two (signed measures), and we recall 
that the first measure on the left hand side and first measure on the right hand side 
are both concentrated on H, whilst the second measure on the left hand side as well 
as the second measure on the right hand side are both concentrated on CH. The only 
way for this to be true is that the two measures that are concentrated on each of H 
or CH are equal: 



^ n 

n=l 



oo 

((pV)*"(-nif))=5:i(M«/"), 



and also 



(-ir 



the former identity being exactly what the statement of the Lemma expresses in the 
language of Fourier transforms. □ 

Based on this result, we may partially generaHse Theorem l.S.ll to higher dimensions 
when we require G (the set that we refer to the exercise region) to be +-closed set. 

Theorem 3.2. Let us make the assumptions nf Lem.m.a \,'^.Sl viz: Let X be the logarith- 
mic price process of the multidimensional Black-Scholes model for independent assets 
with constant volatility and interest rate, that is 



> Xt^[ (Xo), + a, ■ {Bt)^ + ( r - ]^a,^ 



(where B is the d-dimensional Wiener process) for some r > and a G ]R>o . Let 
g = K — f , wherein K > be a real number and f > be a continuous function 
that is monotonely increasing in each component and such that (e~''*/(-'^t))t>o 
martingale. Finally, consider a measurable set of the shape G = j — H' for some 



7 e 



and some convex north-east connected set H' C 



(making G convex and 



south-west connected) such that g = K — f is nonnegative on G and ^ H' . Suppose 
furthermore that both H := —H' and ZH = — CiT' are +-closed. Then for all s > 0, 



OO 

= • E e"™*P'' 

n=l 



oo 

K ■ e-"'*P" 

n=l 



^exp -E 



n-1 



fl {Xu e Cg} n {Xnt e G} 

Li=l 

n-1 

fl {x,t e Ch} n {Xnt e H} 

t=s 



i=l 



tio 



t=s 

tio 



tlO 
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Proof. The existence of limt^o ^g(7) is a consequence of Lemma IT!^ The subsequent 
identity follows directly from Lemma 13.61 (just as in the proof of the corresponding 
equation in Theorem 13. II except that en lieu of Lemma the proof of Theorem 13. II 
makes use of Feller's original result ^1 p. 606, Lemma 3]): For, the second equation 
in Lemma [3.61 may be read 



for all g € (0, 1), in particular for q = e . 

□ 

Again, an analogous result can be accomplished when the function g is replaced by 
f - K axiA the set H by ~H: 

Remark 3.6. In special cases, one can find estimates for f''^ {Xg e H] = P° {Xg G G — 7} 
that are strong enough to establish multi- dimensional generalisations of the estimates 
in the second half of Theorem, \S . 1\ We shall give a few examples. 

In general our results can he used for the extrapolation from (multi- dimensional) 
Bermudan to American harrier knock-in option prices when the harrier regions and 
their complements are, up to a constant factor, closed with respect to multiplication, 
and when, in addition, the harrier region is convex as well as south-west connected (in 
the case of put options) or north-east connected (in the case of call options), and is 
contained in the immediate exercise region of the corresponding American option. 

Example 3.2. Let d he any natural numher. Consider the convex, north-east con- 
nected, -\--closed set H := M<o'' (whose complement is also +-closed) and set G = 
"f -\- H . Let us impose the same assumptions on X , f , and g as in the statement of 
the previous Theorem AS.Sl Then one has, due to the independence of the components 
of X , the following hounds for all s > 0: 



where c,C,a,f3 depend on fi. a — 2 and (3 = 2 "^^"^ if fj. > 0, a = 2 ^''/^ and 
(3 = 2-"^ iffi< 0. Hence also. 




that is 



ri=0°= l<i<n 





K ■ exp 





hm 1/^(7) - ^(3(7) > ^^V^{^) - V6{x) 



for arhitrary s > and a; > 7 componentwise. 
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Example 3.3. Suppose d — 2 and let again X, f, K and g be as in Theorem I.V.ijl 
though we will later on have to impose the condition of ^ > (componentwise). Fur- 
thermore consider the convex south-west connected +-closed set 

H := {{x,y) e : x + cy < 0} 

whose complement is also +-closed for c > and set G = j + H . Note that in this 
situation exp G (which one might refer to as the non-logarithmic exercise region) equals 

expG = • e : u,v > 0, v < u^''} 

fexp denoting componentwise exponentiation as before). Then we get for allt > and 
arbitrary d G N, 

P'' {Xt EG} = P° {Xt e H} (3.28) 



i-i/^B, e i-i/2i7 _ ti/2 (t^'^ 



1=1 ^ 



« / 1=1 



v,/l^{x,y)eW'-.{x-^^^t^'\y-^^t^'^'^e-H 
v^/{{x,y)^R^ ■.x + cy> ^t^'^ + c^t^'^ 

L CTi (72 



(3.29) 

(where vq,! is the normal Gaussian measure on R and i'o,i'' = i'o,i^''y'- 

Next observe that for any real number a > 0, from rotating the set {{x, y) G : x + cy > a] 
by ^, n and we obtain, via exploiting the translation-invariance of the two-dimensional 
normal Gaussian measure I'o.i^, the relation 

4!^o,i^ {[x^v) e ■.x + cy>a]> 2i/o,i^ {{x,y) eB? : x + cy > a] . 

Using the trivial estimate vo,i^ {{x, y) £ : x + cy > a] > vq,!^ {{x, y) € : x + cy > a}, 
we arrive at 

1^0,1^ {ix,y) : x + cy>a} < 1^0,1'^ {{x,y) € R"^ : x + cy > a} 

< 2 • 1/0,1^ {(x,?/) e : x + cy > a} . 

But of course, by a change of coordinates, viz z := ^, one has 

t'o.i^ {{x, y) £ : X + cy > a} — ct'o.i^ {ix,z) e : x + z > a} 

and 

^'o.i^ { [x, y) G K.+^ : X + cy> a] — cv^ i^ {(x, z) G : x + z > a] . 

Applying this to the equation \S. 29\) for {Xt G G} and using the assumption 
H < (componentwise) yields 

ayo,i^{{x,z)eR+^ ■.x + z>a} < P'' {Xt G G} 

< 2cz^o,i^ {(a;,^) G : a; + z > a} (3.30) 
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for 

a 

We can find the following hounds for the measure in the previous estimate: 
Lemma 3.7. For all a > and t > 0, 

1 . ™o , ... 1 



<vo,i{{x,z)eR^ : x + z>a, a:,z>0}<^. 



Proof. The elementary proof has two parts. Firstly, we observe that for all a, z > 

and X € [0, a], 

|2 \ „2 , / \2 



which implies 



\z + a - > + {a - xY 
< / e-=^'/2 / e-'"^^dz dx 

Jo Ja-x 



\ 2 2 / 



e 

Jo 

-a/2 
oo 

< / e"-^dx.e-"V4 = J_,-«V4. (3.31) 



ITT TT 

2 " 75 

Secondly, we have for all a > and a; > 0, 

2x^ + 2o? > |x + af >x^ + a^ 

Thus, 



/TT 

-e = / e c(x • e 







2 

< / e-l^+"l'/2 dx 



which via 



and l|OTll gi-^ 



poo pOQ I pa 



2V2 

< /"e-^'/2 / e-~'-^/^dzdx+ I e-"'/2 / e^^'/^dz dx 



00 



< V"^/2 + JLe-V4<^V2e-"V^ 
~ 2 V2 ~ 



But 



/ e-^'/2^z dx + / e-^'/2 / e~'"^^dz 



e 



e-- /2 / /^dz dx 

J (a-x)va 

p2 



a— a: 



'2iT ■ 1^0,1 {{x,z) eR^ : x + z>a, x,z>0} 
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from which the Lemma follows. 



□ 



This Lemma 's inequalities admit by means of identity kfi. the following conlu- 
sion: 

4v 2 V 2 



By Theorem. \,'^.S\ and the formula ^„ q'^/n — In for all q G (0, 1) 
analogously to the deliberations in the proof of Theorem that 

cV2 



dude. 



I 



K 



1 - e 



) 



< UmF^(7)-y^(7) 



< K ll 



4\/2 



After applying de I'Hospitals rule to the bases of the powers on each side of this esti- 
mate, we get constants Co, Ci > such that for all sufficiently small s, 

Co ■ s'^ < VSii) - VSii) <Ci-s^. 

Example 3.4 (a special Extended Black-Scholes Model). In this example we do 
not assume a multi- dimensional Black-Scholes model, but we presume the discounted 
price process vector S to satisfy the stochastic differential equation 

dSt ^ C ■ Stdt + D ■ StdBt, 

where C, D : fl ^ M.'^'^'^ are mutually commuting symmetric random matrices and B 
is a one- dimensional Brownian motion, subject to the initial condition 

So = 

Then, due to Albeverio and Steblovskaya J2l Proposition 4], we have got an explicit 
solution of that stochastic differential equation, given by 



VtG[0,T] St 



exp ( i • ( C - iz? 



D ■ B, 



exp ( i ( C - -D"^ ] ] cxp {tBt ■D)-e' 



(3.32) 



Then In S and hence the logarithmic non- discounted process In iS+r- are Levy processes, 
and thus Theorem I.V.jjl applies. In this setting we can compute the expression in the 
last line kS. S3p by applying the Spectral Theorem to the symmetric matrices C — ^nD^ 
and tBt ■ D- 
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Chapter 4 

From perpetual to 
non-perpetual Bermudan barrier 
options 



Recall how the function U\T) := U^{T) : ^ R, the expected payoff of a non- 
perpetual Bermudan option on a Feller basket with validity T, log-price payoff function 
g and exercise mesh h as a. function of the logarithmic start price vector, given that 
the option is exercised on the first entry into G C M'*, was defined: 



AT 



behaviour of the difference Uq {Nt) (y) 
remains constant must be the same as the one of the difference Vq (y) 



The purpose of the following Lemma lTTl is to see see that for all y ^ G, the limiting 

UQ{Ns){y) as t tends to zero whilst Nt = Ns 

VQ{y). In words: 

In the (sub-optimal case) of a non-stationary exercise policy for a non-perpetual option, 
the American-Bermudan barrier option price difference has the same limiting behaviour 
as the American-Bermudan difference for the corresponding perpetual barrier options. 

Lemma 4.1. Suppose X is a d- dimensional Feller basket with P' and r > being an 
associated family of risk-neutral probability measures and discount rate, respectively. 
If we define 

fc-i 

V/i > OVfc e N Hk,h := {t^ = kh} = {Xkk. G G} n f| {X, ^ G} , 
then one has for all h > 0, N G N, measurable G C R'' and x G CG, 



V^{x) - U^{Nh){x) 

N 



rNhv 



— e 



3=0 



-rNhm,x 



N 



k=0 



Cg 



[g iXNh)\{XNheCG}] 



Proof For all TV e N, a; e M'^, /i > we can use the Markov property of the Feller pro- 
cess X and the definition of the sequence of events (iJfc_/i)^gpj^ to obtain the following 
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expressions for Ug and Vq: 



JV 



k=0 





" N 






N 














. fc=o 






fe=0 



AT 



= ^ e-'^'=''P-Ffe,^ • [g {Xkh)\ {Xkh e G}] 



k=0 



+e 



-rNhinjx 



N 



L fe=0 



E-[g{XNh)\{XNh&CG}], 



as well as 



fe=0 



J2 e-'''''P" [iffe,/.] • [5 {Xkh)\ {Xkh G G}] . 



k=0 



Then immediately for all h, N, x, 



N 



L fe=0 



■E^ [g{XNh)\{XNh&^G}] 



■rkh^x 



[HkA-^''[9{Xkh)\{Xkh&G}] 



Regarding the difference between Vq and Uq, observe that again for all ft. > 0, iV e N 
and X e M'', 



^ e-'-'^'^P- [Hk,h] ■ [5 (^fcft)l {Xkh e G}] 

/c=JV+l 

CXD 

/c=JV+l 

DC 



= e 



-rA'"?i ^ ^ e~^^^ ■ E^ 



j=0 



JV 



9{Xih)-XH,,,,f]{Xjh^G} 
3=0 



^-rNh^x 



-rNhwax 



N 



{x^h) , n i^jh G} 
3=0 



N 



j=0 



E-[V'^(X;Vfc)|{X;Vfc0G}], 
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where again we have exploited several times the Markov property of X, and - in order 
to interchange E and J2 ^ the assumption that g > on G. 

□ 

This proves 



Corollary 4.1. Under the assumptions of Lemma \4-l\ as well as g = K — f for 
nonnegative / : R"* ^ R and g > on G , we have the identity 



limUUT){x)-U^- 



„-rTTT 



2" 



fl {X,.,-.T ^ G} 

J=0 



ljinV^iy)-V^ 
Cg V*io 



+e" 



. k=0 



G} 



m—n 

■m — >-oo 



limF^(2/)P^^(dy) 
Cg 



[giXr) \{Xt G CG}] . 



<K 



Remark 4.1. Put informally, this Corollary |^.i| means that as soon as one has es- 
tablished order estimates (in n) on the difference \m\tiQVQ{x) — Vq'^ {x) (for in- 
stance the ones from Theorem one only needs to find estimates on the proba- 



bilities 



ana 



)k=0 ^k,2- 



to obtain 



order estimates on the difference limt|o Uq{T){x) — Uq {T){x 
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Part III 



Convergence of some 
approximate pricing algorithms 
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Chapter 5 



Bermudan option pricing based 
on piecewise harmonic 
interpolation and the reduite 

5.1 Introduction 

We intend to approximate the function that assigns the value of a Bermudan option 
with payoff function g and no dividends to the logaritiimic start prices of the underlying 
assets by piecewise harmonic functions. In the first step, we will compute a piecewise 
harmonic approximation to the function that assigns the European option price asso- 
ciated with g and the Bcrmudan's maturity T > to the logarithmic asset prices at 
the penultimate time T — t where exercise is possible. Then we iteratively compute the 
expectation of this function after time t, discount, take the maximum with the payoff 
function g, and perform a reduite-based interpolation (in the one-dimensional setting: 
a piecewise harmonic interpolation). 

Now we would like to answer the following questions: Given the stationarity of 
perpetual Bermudan option prices, can we prove that there exists a minimal fixed 
point of the iteration step described above (which would then be an approximation to 
the perpetual Bermudan price)? If so, can we characterise it explicitly? Is the iteration 
step monotone? 

First, we will discuss these questions in the one-dimensional setting - very little 
knowledge of potential theory has to be assumed for the proofs in that section. Second, 
we shall generalise that approach to higher dimensions; this will entail a few technical 
subtleties. 

5.2 Piecewise harmonic Bermudan option pricing for 
options on one asset 

Consider {ao,...,ar7i} C M, the set of (mutually distinct) support abscissas, and let 
L : C°°(M, R) C°°(M, R) be the infinitesimal generator of a Markov semigroup 
of operators on Lcbcsguc measurable functions from M to M. We call a function 
/ € C°°(M, R) P -harmonic (or shorter: harmonic, if no ambiguity can arise) if and 
only if Lf = 0. Let {Pt)t>o denote the semigroup generated by L. 

For the following, assume L to be a second-order differential operator, that is, there 
are constants /3i,/32 G R such that 

i:/^/3i/' + /?2/". 
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A function g : R ^ M is said to be subharmonic {superharmonic) if and only if g is 
right- and left-differentiable (thus, letting Lg become well-defined as a function from 
K to RU {±00}) and Lg > {Lg < 0, respectively). 

In particular, the supremum (infimum) of countably many harmonic functions is 
subharmonic (superharmonic). 

Lemma 5.1. Given two support abscissas and ordinates, there is a unique harmonic 
interpolation, provided L is a second-order differential operator with a non-trivial 
second-order part (i e ^ 0) or a non-zero first-order part (i e P2 7^ 0/ 

Proof sketch. The uniqueness is a consequence of the maximum principle for harmonic 
functions. The existence follows (in our one-dimensional setting) by distinguishing the 
cases delineated in the statement of the Lemma. If L is a second-order operator and 
it has only a non-zero term of second order, then the space of solutions are ah affine- 
linear functions from R to R. This space is two-dimensional. If there are terms of 
different order, the space of solutions will have basis elements of the form exp(Q!-) and 
we have to solve a linear or quadratic equation to find the a (or a's) satisfying this 
linear or quadratic equation. Since L is sub-Markovian, there will be at least one real 
solution to this equation for a. □ 

The Lemma implies 

Corollary 5.1. There cannot be more than two linearly independent harmonic func- 
tions: There is a canonical monomorphism from the space of functions to the - two- 
dimensional - space of pairs of subordinates. 

Lemma 5.2. A subharmonic function from R to R is constantly zero if it has three 
zeros. 

Proof. The left- and right-differentiablility of subharmonic functions entail that for all 
subharmonic g, Lg will be defined as a function from R to R U {±00}. □ 

If there is only a first order non-zero term, the space of harmonic functions will 
just coincide with the space of constant functions. 

Lemma 5.3. 1. Piecewise harmonic interpolation with respect to the support ab- 
scissas {ao,...,am} preserves subharmonicity on [ao,am]; The interpolating 
function dominates the interpolated function on A :— [00,0™], and if the in- 
terpolating function f equals the harmonic function fi on [a;, ai+i] for all i < m, 
then we have f = sup{/o, . . . , fm-i}- 

2. The interpolating function f is strictly dominated by the interpolated function 
T{f) on the intervals (— oo,ao) and {am, +00). 

Proof sketch. 1. The domination part follows from the maximum principle for har- 
monic functions. From the maximum principle, we also get for alH < m — 1 that 
if fi 7^ fi+i, then 

{ft = fi+i} = {at+i} ■ 

Now there are two possibilities: either fi < fi+i on (—00, ai^i) and fi > fi+i on 
(oi+ij-l-oo) or the other way round fi > fi+i on (—00,0^+1) and fi < fi+i on 
(oi+i, -foo). However, in the former case, the interpolating function would equal 
fi A fi+i on [0^,0^+2], which is superharmonic, and it would also dominate the 
subharmonic interpolated function 77 on [ai,ai+2]- Then, rj — {fi A fi+i) would 
be nonpositive and subharmonic on [0^,0^+2] and it would have three zeroes, in 
Ui, Qi+i and a,i+2. By Lemma this can only be true if 77 — {fi A fi+i) = on 
[ai, ai+2]- Thus, rj — fi A fi+i on [oi, 0^+2]- Since 77 is subharmonic on [a^, 0^+2], 
so must be fi A fi+i then, and therefore, fi A fi+i is harmonic on [0^,0^+2]. 
This means fi = /i+i (as both fi and /i+i are harmonic) which contradicts our 
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assumption that fi ^ fi+i. Therefore, > fi+i on (-oo,ai+i) and fi < 
on (fli+i, +00) for all i < m. 

Inductively, this yields f > fi on [ag, a^] for all i < m, hence / sup{/o, . . . , fm-i} 
on [ao,am]- 

2. The function / — is subharmonic on (—00, ai] and it has two zeroes in oq 
and oi. Moreover, it is nonpositive on (ao,ai). Because of Lemma [5. 21 then 
/ — I{f) has to be positive or negative on (—00,00). In the former case, we 
are done. In the latter case, due to the maximum principle, / — must be 
decreasing and therefore in oq we would have L{f — X{f)) (oq) < 0, which is 
absurd. A symmetric argument works for the proof of the domination of 
by / on the interval (a„i, +00). 

□ 

Lemma 5.4. Piecewise harmonic interpolation to a set of support absicssas {ag, . . . , 
is monotone on [ao,am] in the sense that if f < g on [ao,am], then the piecewise har- 
monic interpolation of f will be dominated by the piecewise harmonic interpolation of 
g on [ao,am]. 

Proof. Use the maximum principle on each of the intervals [a^, Oj+i] for i < m. □ 

Lemma 5.5. Let X : rI^o^"™! denote the operator of piecewise harmonic inter- 

polation with respect to the set of support abscissas {cq, . . . ,am}. Let / : M M fee 
subharmonic on M. Consider a harmonic function h : R ^ M., assumed to dominate 
f: f <h on R. Then I{f) <h onM.. 

Proof. From the previous Lemma [5.41 we already know that 2{f){x) < I{h){x) holds 
for all X g [ao, However, I{h) = h, hence < h on [oq, a,„] and from Lemma 
15.31 we conclude that h> f > T{f) on the intervals (—00, ao) and (a™, +00). □ 

Theorem 5.1. LetT : M}°-o,am] ^ j^R g^ajjj denote the operator of piecewise harmonic 
interpolation with respect to the set of support abscissas {ao,...,am}. Let g be a 
subharmonic function, let c be nonnegative and subharmonic, and let h be harmonic. 
Let c be, moreover, harmonic on each of the intervals [aijfli+i] for i < m. Suppose 
c < g on [ao, and c, g < h on M., r > and let t > 0. Now define 

K-.f^I {e-'-'Pt (!(/) V c) V g) r [ao, a^] 

as well as 

/ : R ^ M subharmonic, f > c on [ao, a^], 1 
Vi e {1, . . . , m — 2} f harmonic on [ai,ai+i], >. 
/ harmonic on (— 00, oi), (a„i-i, +00), f < h J 

Then K maps the convex and bounded subset Q of C°[ao,am] continuously to it- 
self. Moreover, due to Lemma \ 5.1\ Q is a subset of a finite- dimensional subspace 
o/C°[ao,am] (this subspace being the space of all functions from [ao,a„i] that are har- 
monic on each of the intervals [ai, a^+i] for i < m. By Brouwer's Fixed Point Theorem, 
K, has got a fixed point in Q. Finally, K is a composition of monotone functions on 
[aojOm] and therefore monotone as well. 

Proof sketch. We can divide the proof for 1C{Q) Q Q into three parts: 

1. The cone of subharmonic functions is closed under V, under Pt, under multiplica- 
tion by constants and under piecewise harmonic interpolation J (cf Lemma lOll . 
therefore the image of Q under /C can only consist of subharmonic functions. 



Q \ f \ [ao,am] 
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2. The upper bound on the elements of the image IC{Q) follows from the mono- 
tonicity of Pt and I fLemma I5.4|l . combined with the equations Pfh — h and 
I{h) — h as well as the Lemma[01 First, we may state e~^^Pt {I{f) V c)V g < h 
for all / < /i, which by Lemma allows us to deduce 

I {e-^'Pt {I{f) Vc)Vg)<h 

for all f e Q. 

3. The lower bound follows again from the monotonicity of I, but this time only by 
exploiting c < g on [gq, a„i] and employing the fact that the space of those func- 
tions that are harmonic on each of the intervals [0^,0^+1] for i < m is invariant 
under the composition of X with the restriction to [ao, am] (yielding X(c) = c on 
[ao, flm]). 

Since c is nonnegative, we get that Q is bounded by supj^^j „^ j /i > as a subset of 
C°[ao,am], and because Q is finite-dimensional, we may apply Schauder's Theorem, 
provided we are given the continuity of /C. However, this last assertion follows from 
the maximum principle. □ 

The existence of a minimal fixed point for K, can be proven constructively as well: 

Corollary 5.2. Let us adopt the notation of the previous Theorem. Then the sequence 
(/C"((7 V 0))„gpj^^ is monotone on [ao,am], bounded and dominated by h. Therefore we 
have the existence of a limit on [ao,am] given by 

Vx e [ao, a„] q{x) lim K."{g V 0)(x) = sup /C"(g V 0)(x). 

This limit is an element of Q and therefore can be canonically extended to the whole 
o/M. By the continuity of IC, q is a fixed point of IC. On [ao,a„i], the convergence in 
the last equation will be uniform. 

Proof. The only part of the Corollary that does not follow directly from the preceding 
Theorem 15. II is the uniformity of the convergence and that q will be harmonic on each 
of the intervals [0^,0^+1] for i < m. However, monotone convergence on compact sets 
preserves harmonicity and is always uniform (cf e g Meyer [2] - or, more directly. 
Port and Stone [251 Theorem 3.9] if P is the Brownian semigroup). □ 

Lemma 5.6. In the preceding Corollary's notation, q is the minimal nonnegative fixed 
point of IC. 

Proof. The proof partly copies the one for Lemma 11.41 Any nonnegative fixed point 
p of IC must be greater or equal g on [ao,am]. Therefore the monotonicity of K. on 
[ao,am], implies 

VneNo p = /C>>/C"(gVO) on [ao,am], 

yielding 

P > sup /C"(5 V 0) = g on [ao, a^]. 

nSNo 

□ 

Example 5.1 (Bermudan vanilla call on a dividend-paying asset in a special 
Black-Scholes model). Assume 

P ■■= {Pt)t>0 {l^t,t,aH * ■)t>o ' 

where 

cr > 0, fi:^r-—, 
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thus P can be perceived as the semigroup associated to the logarithmic price process 
under the risk-neutral measure in the one- dimensional Black-Scholes model). We will 
assume that (possibly after a linear change of the time scale) a = 1 and we assume 
that r has been cut to discount dividends. Define 

g := exp —K 

(the payoff on exercise of a one- dimensional call option with strike price K). The 
infinitesimal generator of the Markov semigroup P is 

Thus we obtain 

Lg = rcxp > 0, 

hence g is, P-subharmonic. We can find the P -harmonic functions for ^ ^ ( otherwise 
they are simply the affine linear functions) by observing that for all a e M, 

= Lexp(a-) = — (a^ + 2iiaj exp(a) 
<^ a e {0, -2fj,} . 

If fj, 7^ 0, the functions 1 : a; i— > 1 and exp(— 2/i-) : x i-^ ^-'^f^^ qj-^ i^^q linearly 
independent harmonic functions, thus by Corollarv \5.1\ we have already found a basis 
for the space of harmonic functions. If fj, = 0, the harmonic functions are exactly the 
affine linear functions. In order to obtain the setting of Theorem \5.1\ we will assume 
/J, < ^ such that the sum h o/exp(2/i-) and a sufficiently large positive constant is a 
harmonic function dominating g — exp— if . In order to satisfy the conditions on c we 
could simply take c = —K for instance. 

5.3 Reduite-based approximation of Bermudan op- 
tion prices 

Suppose P is a Markov semigroup on [d € N) and L is the infinitesimal generator 
of P. We will call a function / : M'* ^ M subharmonic if and only if 

Vt>0 Ptf>f 

holds pointwise. A function / : K'' ^ M will be called superharmonic if and only if 
— / is subharmonic, and / : M'' ^ M will be called harmomic if it is both super- and 
subharmonic. 

Let U denote the operator of upper-semicontinuous regularisation, that is, for all 
functions / : R'' ^ R, 

Uf ^ini{i> f : ^ : M'' ^ M subharmonic} 

(of course, this is a priori only defined as a function taking values in R U { oo}). 

Consider a harmonic function /i : R'* — > R and a closed (and therefore F^) set B 
and define the reduite operator Tl — TZh,B on the set of all subharmonic functions 
/ : R'' ^ R dominated by h via 

TZf := U (sup {e<h : ^ : R'' ^ R subharmonic, i < f on B}) . 

It is a well-known result from potential theory (cf e g the work of Paul- Andre Meyer 
[22I Theoreme T22]) that there will be a greatest subharmonic function dominated by 
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f on B and that this function will be equal to TZf. Moreover, we have that / = TZf 
on B except on a set of potential zero, in probabilistic/potential-theoretic jargon 

/ = 7^/ q.e. on B, 

where "q.e." is, as usual, short-hand for "quasi-everywhere". Now define 

Q ■■= {f <h : / : M"* ^ R subharmonic} . 

Then our definition of the reduite operator TZ implies TZf < h {as h is dominating the 
function whose upper-semicntinuous regularisation is, according to our definition, the 
reduite TZf of /) and our potential-theoretic characterisation of the reduite - as the 
greatest subharmonic function dominated by / on S - ensures the subharmonicity of 
TZf. Therefore, 

TZ-.Q^Q. 

We also have that U is monotone (in the sense that for all /o < fi, Ufa < Ufi) so 
that TZ must be monotone as well (from the C-monotonicity of sup and the definition 
of TZ). 
Hence 

Lemma 5.7. Adopting the notation of the preceding paragaph, TZ : Q ^ Q and when- 
ever /o < /i, TZfo < TZf I. 

Let g : M'' ^ R be a subharmonic function such that g < h and let r > 0. The 
next step is going to be the consideration of the following family of operators: 

for t > 0. If f < h, Ptf < Pfh — h for all t > 0, since the operators Pt are positive and 
linear, and h was assumed to be harmonic. Thus, since g < h and r > 0, one must have 
4>tf <h for all / < ft. and i > 0. Moreover, the operators Pt preserve subharmonicity 
and the maximum of two subharmonic functions is subharmonic again, therefore <j)tf 
must be subharmonic for all subharmonic /. Finally, since Pt is monotone, 04 has to 
be monotone for alH > Summarising this, we obtain 

Lemma 5.8. Using the notation introduced previously, 4>t '■ Q ^ Q and whenever 
f(><fi,4>tf^<<t>tfifor allt>0. 

As a consequence, we derive from the two Lemmas 15 . 71 and 15 . 81 the following: 

Corollary 5.3. If we define K-t TZ o (jjt (adopting the notation of the previous 
paragraph), we have K-t Q ^ Q and whenever /o < /i, ICtfa < /Ct/i- 

Corollary 5.4. The map f i-^ ICtf V is a sound iterative Bermudan option pricing 
algorithm for the payoff function g V (in the sense of Deftnition \1.6\) . 

This already suffices to prove the following 

Theorem 5.2. Lett>0. Then for all n e Nq, 

/Ct"+i(5V0)>/C4"(5V0). (5.1) 

Furthermore, 

q := sup /Ct"(5 VO) 

neNo 

(which a priori is only defined as a function with range m R U {+oo}^ is an element 
of Q and indeed is the least nonnegative fixed point of K-t . 

Proof. 1. Relation H5.1II follows from the fact that Kt is a sound algorithm and 
R emark 1 1.41 
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2. Since K-t maps Q to itself, the whole sequence (/Ct"(g V 0))^g[^^ is bounded by h. 
This entails q < h a.s well. Applying Beppo Levi's Theorem on swapping sup and 
/ -dfj, - for bounded monotonely increasing sequences of measurable nonnegative 
functions and an arbitrary measure /i - to the measures Pt{-,x), x G R"* and the 
sequence (A^t"(5 V 0))„gpj^ , we can exploit the subharmonicity of the functions 
/Cf"(g V 0), n e No, to deduce 

Vx G Ptq{x) = sup Pt {JCt"{g V 0)) (x) 

neNo 

> sup /Ct"(5V0)(a;) =g(a;), 

nSNo 

which is the subharmonocity of q. As we have already seen, q < h, so q G Q. 

3. K we employ Beppo Levi's Theorem again, we can show that ICt and sup^gp^j^ 
commute for bounded monotonely increasing sequences of functions. Thereby 

ICtq^ sup /Ct/Ct"(.gVO) = sup/Ct"(gVO) = q. 

4. That g is the least nonnegative fixed point is seen as in the proof of Lemma, lOl 
Any nonnegative fixed point p of K-t must be greater or equal 5 V 0. Therefore 
by the monotonicity of sup and ICt , 

sup /Ct> > sup /Cf"(5 V 0) = q. 

□ 

Example 5.2 (Bermudan call option with equidistant exercise times in i • Nq 
on the weighted arithmetic average of a basket in a special Black-Scholes 

model). Let (3i,...,l3d £ [0,1] fee a convex combination and for simplicity, assume 
that the assets in the basket are independent and each follow the Black-Scholes model 
with one and the same volatility cri = ■ ■ ■ — Ud —'■ <J , and let r > be the interest 
rate of the bond. We may assume that, possibly after a linear change of the time-scale, 

(7 = 1. Then {Pt)t>o = ( '^t(j._i]'' 1 1* ' I semigroup of this Markov (even 

Levy) basket. Then one has 

(cf e g Revuz and Yor's exposition ]EEj), and for 

d 

g : xi-^'^/3i exp (x^) -~ K 

i=l 

we obtain 

L9 = i:(^+(r-^))c.pii-).) 



i=l ^ ^ ' ' 

which is pointwise nonnegative if and only if 

max,,g{i,...,rf}A _ 1 EtiA' 
2 2d 2d ' 

Hence, if r is sufficiently large, g is subharmonic and we can apply the theory developed 
earlier in this Chapter, in particular Theorem \5.'A 



54 



Chapter 6 



Soundness and convergence rate 
of perpetual Bermudan option 
pricing via cubature 

When Nicolas Victoir studied "asymmetric cubature formulae with few points" (22] 
for symmetric measures such as the Gaussian measure, the idea of (non-perpetual) 
Bermudan option pricing via cubature in the log-price space was born. In the following, 
we will discuss the soundness and convergence rate of this approach when used to price 
perpetual Bermudan options. 

Consider a convex combination (ai, . . . , am) G [0, 1]'* (that is, X^fcLi Q^fc = 1) and 
xi, . . . ,Xm € IR''- Then there is a canonical weighted arithmetic average operator A 
associated with a, x given by 

m 

V/eM« Af = Y,akf{--Xk). 

k=l 

Now suppose c € (0, 1), g,h : R''^ ^ Ag > g, Ah ^ h and QW g < h. Define an 
operator V on the cone of nonnegative measurable functions by 

V: f ^{c-Af)\J g. 

Since A is positive and Hnear, thus monotone (in the sense that for all /o < /i, 
Afti < Afi)^ it follows that V must be monotone as well. Furthermore, whenever 
Af > /, we have that AVf > Vf, as the linearity and positivity of A combined with 
our assumption on g imply 

AVf>cA^f\/g>{c-Af)\/g = Vf. 

Finally, due to our assumptions on h and g, we have for all nonnegative f < h, 

Vf < cAh < h. 

Summarising this, we are entitled to state 

Lemma 6.1. Adopting the previous paragraph's notation and setting 

Q-={f<h : Af>h}, 

we have that 

V is monotone (i e order-preserving), and AD — V is nonnegative. 
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Corollary 6.1. The map f i-^ Vf V is a sound iterative Bermudan option pricing 
algorithm for the payoff junction (? V (in the sense of Deftnition \1.6p . 

This is sufficient to prove 

Theorem 6.1. For all n e No, 

I?"+i(.gVO)>I?"(.gVO) =:(?„. (6.1) 

Furthermore, 

q lim ^"(.g V 0) = sup V"{g V 0) G Q 
and q is the smallest nonnegative fixed point ofV. 

Proof. 1. Relation H5.1|l follows from the soundness of I? is a sound algorithm and 
Remark 11.41 

2. Since T> maps Q itself, the whole sequence {V^ig V 0))^gp,^ is bounded by h. This 
entails q < h as well. Using the linearity of sup and our previous observation 
that AD — I? > (Lemma inni) we can show 

WyeR^ Aq{y) = sup A (P"(g V 0)) (y) 

neNo 

> sup I?"(.gVO)(y) 

which means Aq > g. As we have already seen, q < h, so q Q. 

3. Again, due to the linearity of sup and the special shape of V that is based on 
a weighted arithmetic average operator, V and sup„gpj|^ commute for bounded 
monotonely increasing sequences of functions. Thereby 

Vq = sup X>X>"(g V 0) = supX>"(5 V 0) = g. 

4. Just as in the proof of Lemma lTT^ we see that q is the minimal nonnegative fixed 
point. For, any nonnegative fixed point p of V must be greater or equal 5 V 0. 
Thus, by the monotonicity of sup and V, 

sup !?> > sup V{g V 0) = (?. 

□ 

Lemma 6.2. Using the previous Theorem's notation, we have for all x € R''' and 
n e No, if qn+i{x) = g{x), then q„{x) = g{x). 

Proof. By the monotonicity of the sequence {qn)n&io (^Theorem 16. Ill , we have 

9{x) < qoix) < qn{x) < g„+i(a;). 

□ 

Theorem 6.2. For all n G N, 

WQu+I — <?n|lcO(R<*,K) — ^ ' ll^n " 9"- 1 1 1 (R<* ,R) ' 
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Proof. The preceding Lemma Fft . 21 yields 

Win+l — 9n|lcO(R'i,R) = \\ln+l — 1n\\ CO ({c- Aq,^> g] M) 

= \\C-Aqn- ((C • Aqr,-l) V g)|lcO({c.Ag„>ff},R) 

via the definition of <7i+i as [cAqi) V g for i = n and i = n + 1. But the last equality 
implies 

Ikri+l — 9n|lcO(R<i,R) ^ II ^ ' " C ' ^9n- 1 1| C0({c-ylg„ >g} ,R) 

< ||c • A(7„ — C • A(7„_i||pO(i{d R) ■ 

Since A is linear as well as an -contraction (and therefore a C°-contraction, too), 
we finally obtain 

||<Zn+l ~ 'Zn|lcO(R'i,R) — 11^ (^n ~ 9"- 1 ) II C" (R"* ,R) — ^ 11*?" ~ 9"- 1 II (R'' ,R) ' 

□ 

Example 6.1 (Bermudan put option with equidistant exercise times in i • Nq 
on the weighted arithmetic average of a basket in a discrete Markov model 
with a discount factor c — e^''* for r > 0). Let Pi,..., (3d G [0,1] fee a convex 
combination and assume that A is such that 

ra 

Vi€{l,...,4 ^a,e-(^-)' =1, (6.2) 

k=\ 

then the functions 

d 

g : X 1-^ K - ^^f3i exp (xi) 

i=l 

and h :— K (where K > 0) satisfy the equations Ah — h and Ag — g, respectively. 
Moreover, by definition g < h. Then we know that the (perpetual) Bermudan option 
pricing algorithm that iteratively applies V to the payoff function gWO on the \og-price 
space, will increase monotonely and will have a limit which is the smallest nonnegative 
fixed point of V. Moreover, the convergence is linear and the contraction rate can be 
bounded by c. 

The condition \f). ij)) can be achieved by a change of the time scale (which ultimately 
leads to different cubature points for the distribution of the asset price) 

One might also be interested in determining the convergence rate for the approx- 
imation of non-perpetual American option prices based on non-perpetual Bermudan 
option pricing via cubature. After proving a series of Lemmas we will end up with a 
Theorem that asserts linear convergence and also provides bounds for the convergence 
factor. 

From now on, c and A will no longer be fixed but their role will be played by 
e"""* and Pt (for t G sNo where s > shall be fixed) respectively, where r > and 
{Ps-m)mmo describes a Markov chain on W'- (By the Chapman-Komogorov equation 
this is tantamount to \ls,t > PsPt = Ps+t)- 
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Chapter 7 



Some convergence estimates for 
non-perpetual American option 
pricing based on cubature 

For this Chapter, let us consider an arbitrary but fixed translation-invariant finite- 
state Markov chain P := [Ptj^^i with state space (for d € N) where / = /iNq 
for some real number /i > 0, as well as a real number T > h > (the time horizon, 
or maturity), a real number r > 0, a continuous function / : R'' ^ [0,-|-oo) that is 
monotone in each coordinate, a nonnegative real number K >0 and let us set 

9:=K- f 

as well as defining a family of maps Bt : L° (M'', [0, -hcxo)) L° (M'', [0, -|-cxd)) , i > 0, 

by 

Vt > Bf.f^ max {e-'-*Pt/, g} = [e-'^'Ptf) V g. 
(Note that Btf will always be nonnegative for / > - hence, for all / > 0, Btf > gVO.) 
Furthermore, we shall denote by \xi \ . . . ,x^^\ \ the set of (distinct) states at time 

t after starting the process at time in e M** and by \oii \ ■ ■ ■ ,ct^*\ | C (0, 1] the 

weights for each of these states, thereby imposing on the sets I a^i\ . . . , a^*\ \ for 

t e I, in addition to it being a subset of (0, 1], the condition that they be a convex 
combination, viz. 

t 

Vie/ ^af^ = l. 

i=l 

Summarising this, we write 

t 

i=l 

For the whole of this section, the Lebesgue measure on M shall be denoted by A, and 
A"^ will be shorthand for the measure-theoretic power A®"^. 

The operators max and min when applied to subsets of will be understood 
to be taken componentwise. Analogously, we will interpret the relations < and > 
componentwise on R"^ x R"^. 

For convenience, we allow all i^'-norms (including the C° norm) of measurable 
functions to take values in the interval [0, -|-oo], thereby extending the domain for 
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each of the L^-norm to the vector lattice of measurable functions. Furthermore, 
any functions occurring in this Chapter will be assumed to be measurable. Thus, eg 
the relation /o > /i should be read as shorthand for /o G (A*^) fl {• > /i} for all 
functions /o,/i; analogously for the relation /o < /i- 

Finally, we will use the operation V in such a way that it is applied prior to +, but 
only after Pg and multiplication with other functions or constants have taken place: 

C ■ Psh V /3 • /i + /2 = max{(7 • PJo, h ■ fi] + f2 

In this Chapter we are aiming to understand the convergence behaviour of the 
sequence ^(-By.2-n)°^ {g V 0)^ . We will start by noting that this sequence is mono- 
tonely increasing: 

Lemma 7.1. The sequence ( (-B^.o-n)"^ /) is monotonely increasing 

\ /neNon{T-2-6/} 

for all functions / : M'' — > K. Furthermore, if there exists a function g > g ^ 
such that g is e"*"' P. -harmonic (ie e~'^^Phg = g) and f < g, then for all n € No, 

(St.2-0°'" / < .9- 

Proof Consider / : R'' ^ M and n € No such that T • 2-("+i) € / = /iNo. Then 

(Pt-2-("+i))°^ / = (^(-Bt.2— '+0°^) / 

and by the monotonicity of the operators Pg for s € I, 

(B^.2-„+0°' = e-'-^-^-<""''p^.2-(n+.) (e-'-^-2-"*^^'p^.2-(n+.)(-) V5) Vg 

> e-^-2-^"^^'p^.2-(n+., (e-^-^-'-^^'p^.^-cn+.iO) V5 

= e-'-^-2""p^.2-„(-)V5 = BT.2-", 

where the last line is a consequence of the Chapman-Kolmogorov equation. This com- 
pletes the proof for the monotonicity of the sequence ^(-87^.2-")°^ ^ y 

Now suppose there exists such a function g as in the statement of the Lemma. 
Then e^'^'^Psg = g for all ,s G / and therefore Bsg = g for all s G I. Also, the map Bs 
is monotone in the sense that go < gi always implies Egga < B^gi (because it is the 
composition of two monotone maps: e~'"*Ps and • V5) for all s & I. Thus we see that 
for all f < g, 

(PT.2-)°'"/<(ST.2-r'''9 = 5- 

□ 

Lemma 7.2. For all measurable functions /i > /o > 5 V 0, as well as for all t € I 

and p e {1, 00} one has 

ll-Bt/i - Pt/o||£,p(;^a) < e"''* ll/i - /o||^p(;^<i[{e-rtPj/i>g}n.]) 

< e"''* ll/i - /o||iP(A<i) 

(with the usual convention that x < +00 for a// a; S M U {±00} ). 

Proof. The map Bt is monotone. Thus we have 

{Btfi=g} = {Btfo < Btfi = g} 

= {5 V < Btfo < -Bt/i =g} = {Btfo =g}n {Btfi = g} 
C {Btfi-Btfo = 0} 
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for /i > /o > .9 V 0. Since Btfi > g, this implies 

Q<Bth~Bth = X{B,f^>a}{e-^'Pthy g-e-'-'Ptky g) 
= X{B,f,>g} {e-^'Pth - e-^'Pth V g) 
< X{e-p,/,>s} {e-^'Ptfi - e--*PJo) 

= e"''*X{e-'-*Pi/i>g}-Pt (/l - /o) 

which yields the assertion as Pt is an L^'(A'')-contraction (for p = oo this is immediate 
and for p = 1 it follows from the translation-invariance of both Pt and the Lebesgue 
measure). □ 



Lemma 7.3. Suppose 



max Xj, < U 
fee{i,...,m} 



componentwise, implying x'f' < componentwise for all s ^ I = /iNo anii i G 
|l, . . . , TO'^ } . T/ien /or a/Z s £ I , g is nonnegative on {Ps{g V 0) = Psg}- 

Proof. Recalling our notational convention that < as relation on x and max 
when applied to subsets of M'' are to be interpreted componentwise, we may write 

it) ^ Jt) 



Vie {l,..., mi} R'^ 9 > max ^ a;f ^ > x 



£ ' 



ie| l,...,m'H' I 

due to the componentwise monotonicity of g, yields for all s S / the inclusion 
{PsigWO)^Psg} - {Vze {l,...,m^} g (. - xj^') > o} 

C \gl- max 1 > 

C {5>0}. 

□ 

Lemma 7.4. Suppose there is a 70 > 1 (without loss of generality, 70 G [l,e'^)) such 
that 

Ptf > 70*/ 

for all i e (0, r] n/ (where I = /iNo wii/i /i > whence it is sufficient that this estimate 
holds for t = h). In addition, assume that g > on the subset {Pt{g V 0) = Ptg} ofM.'^ 
for all t £ (0,r] n / (this assumption being, due to Lemma l7.'A satisfied in particular 
i/maxig{i_..._^}xf ^ < 0;. Then for all t e (0,T] n /, 

{g>e-'-*Ft(gVO)} D {Pt{g \J Q) = Ptg} 

Proof. Let t G (0, T] n /. Due to our assumption of 5 > on {^'4(5 V 0) = Ptg}, one 
has 

{5 > e-'-'Ptig V 0)} n {Pt{g V 0) = Ptg} 
= {g>e-''Ptg]r\{Pt{gyQ)^Ptg} 
= {g > 0} n {(id - e-^'Pt) .g > 0} n {Pt{g V 0) = Ptg} 
= {5 > 0} n { (1 - e-'-*) if > (id - e-^'Pt) /} 

n{Pt{g\JO)^Ptg}. (7.1) 
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On the other hand 

{(1 - e--') if > (id - e-'-'Pt) /} = {{l- e-^') K > / - e"''* Pj 

>7o'/J 

3 {(l-e-'-*)if>(l-e-'-V)/}, 

that is 

{(l-e-')A-> (id-e-'P,)/) 

where we have exploited 70 < e*". Now 70 S [1, e' ) gives 

1 — e ''*7o 

since K >0. Combining this estimate with the previous inclusion, one obtains 

{(1 - e--^*) if > (id - e-'-'Pt) /} D {if > /} 

and hence 

{ (1 - e-'-') K > (id - e-^'Pt) /} r^{g>{)} = {g> 0}. 
This result, combined with the first equation Ij7.1|l in this Proof, yields 

{5 > e-'-'Ptig V 0)} n {Pt[g V 0) = Ptg} 
= {g>Q}n{Pt{gyQ) = Ptg}. 

However, one of our assumptions reads 

{Pt(.gV0)=Pt5}C{g>0} 

whence we conclude 

{5 > e-'-*Pi(g V 0)} n {Pt(5 V 0) = Ptg} 
= {Pt{gyQ)^Ptg}. 



□ 



Remark 7.1. The assumption of the existence of a < such that Ptf > 70*/ for 
allt^I is natural: If {Xt)tei "was a Markov process evolving according to {Pt)tei, the 
(stronger) condition 

ytel Pj^e+^'f 

simply means that the process f{X.) is, after discounting, a martingale. Now, if X. 
was a Markov model for a vector of logarithmic asset prices (a Markov basket in 
our terminology) and f would assign to each vector the arithmetic average of the 
exponentials of its components, this is by definition true if P. governs the process X. 
under a risk-neutral measure. Furthermore, the said assumption 

ytel PtJ^e+-'f 

trivially implies 

371 >o vtei PtI<ii'I 

and therefore provides us with some vindication for assuming the last assertion in some 
of the subsequent Lemmas of this Chapter. 
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Lemma 7.5. Suppose there is a 71 > such that 

for all t £ (0, T] n / (for which in our case of I ^ HNq with h > it is sufficient that 
this estimate holds for t = h), and let us assume without loss of generality that this 71 
be > . Then, setting 

71* - 1 

R:^K- sup — , 

te{o,T]ni t 

we have found an R < +00 such that for all s e (0, T) n (2 • /) C / and measurable 
/>5V0, 



<R- 



Proof. Let s G (0, T] n (2 • /) C / and consider a measurable f > gV 0. Then by our 
assumption of Ptf < 71*/ for 71 > 0, we firstly have (inserting | for t) 

g~e-^ip.g = K- / - e"''* {K - P.f) 
= if (l-e-'^t)+e-'^iPj/-/ 
< iffl-e-"*) + (e-''47i^ -1) •/ 



and therefore (using f > gV and 71 > as well as the monotonicity of P^ 



< X| 
< Xl 



g>e ^2P^ 



(g-e-^'^P.f) 



{9 



[a>e-'-iPsf} 

< X{g>o} -{Kil- e-^^) + (e-'-47i4 - l) • /) 

< X{f<K} -{Kil- e-'-i) + {e-^h^i - l) ■ K) 

< K ■ {(l^e-'i) + {e-'-i-fii 
= K ■ e-''^ ■ hi^ - 1) 



(7.2) 



Now, 



sup 



71 



te(o,T]n/ t 



< +00 



since i 1-^ 71* — 1 is right-differentiable in zero with derivative In 71. Therefore 



71-1 



R — K- sup 

te(o,T]n/ * 



< +00. 



Via estimate H7.2II . we arrive at 
But 



(5-e-fp./)<i?.-. 



(7.3) 



XI 



V] 



{g - e-^iPsf) = {e-^iP.f) V 5 - e'^ip^f, 



and - in combination with the linearity of P| and the Chapman-Kolmogorov equation 
- this implies 
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hence by equation IjT.HII 



3 ' 



> e-'-5Pj((e-'-tp./)V5)-e-"P,/V5 (7.5) 

(where in H7.5II we have exploited the fact that P| is an L°° -contraction). 

Now, again by the Chapman-Kolmogorov equation and the monotonicity of Pt for 
any i G /, we have 

and therefore (due to the estimate Btfo > 3 V which holds for arbitrary < G / and 
/o>0) 

{{B,rf-9} - {Bsf<{B,rf = 9} 

= {BJ = g}n{{B^y'f^g} 
C {{B.y'f-BJ = 0} 



But 



thus the last inclusion yields 



(B^y f>g\/0>g, 



< {B.y^f-BJ 



, o2 



= e-i (P. (e-tPj/ V,,) - (e-P.fVg)) 

where the last line has used the estimate (j7.5|l derived previously. □ 

Later on, in Lemma [7.61 we will see that it is impossible to obtain estimates for 
(B^)°^ f — Bsf that are both uniform in / > o V and of higher than linear 

order in s. 

We can draw from the proof of Lemma 1731 the following Corollary: 
Corollary 7.1. Suppose there is a 71 > such that 

for all t e (0,T] n / (for which in case I ~ /iNq with h > it is sufficient that this 
estimate holds for t = h). Then for all measurable f > gV 

< {B^y^f-BJ 

= ((^~^*^f /) V 9 e-^iPsf) x| ^^^^ 



(^{.>e-tp,(.vo)}-(.9-^-^^^f(5V0))) 



{B^y (5V0)-P,(5V0). 
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We will continue to assume / = hNo, h < T and define 
Definition 7.1. 

\ftGl E' ={Pt{g\/0)> Ptg}. 
Remark 7.2. Equivalent expressions for E' are: 

yt&I E* = {Pt{gVO)> Ptg} 

= L{Pt{gVO)=Ptg} 

= {Pt(5A0)<0} 

= {3iG {l,...,m*} .g(--xf^)<0| 

= u {^</(---'0}- 

i6{l,...,m7r } 

These formulae for E' imply, by the monotonicity of f, that E^ is north-east connected 
(that is E!^ C E* + a for all a < 0) for all t G I. Furthermore, if one had for all t € I 

and i G |l, . . . ,m^| an index k = k{i) G |l, . . . ,m^| such that xf^ + ar^^j < 

componentwise (for instance if the set ^Xi^\ . . . ,Xm^^ could be written as the sum of 

a reflection symmetric subset o/M'' and a componentwise nonpositive vector), then the 
north-east connectedness of the E* 's entails for all t € I and i, 

= U ^ {^</<(--f--S))} 

je{l,...,m^:} 

^ u ^ {^</<(--4:^-4f)} 

= U {^</<(-4"0} 

£e{i,...,mTr} 
= E'^' 

where for the penultimate line we have used the Chapman-Kolmogorov equation, of 
course. Therefore 

Xe* (• + xf^^ = X^t.^W < XB^t- 
Also, if there exists an io € {1, . . . ,m} such that 



< 

■ J 



yje{l,...,d} {x\'J)^ 

has - due to the monotonicity of f in each coordinate - first of all f ^- — > / 
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and thence for all n €N the inclusion 

E"'' = y [k < f - xi''''^)} 



ie{l,...,m"} 

u {^</( 

iiT..,in€{l,...,m} 
ii,...,i„_iG{l,...,m} 

U {^</(--^' 



ii,...,i„_ie{l, •••>"»} 
^(n-l)/» 



{ A' </(.-x<;" -<>,)} 



Tftis means 



and for all T € [h, +00], 



i?'* t as s t 00 in / 

Pi = E^. 

se(o,T]n/ 



The reason for E'^ not being the whole space is that the measure B 1— > PjXs on 
the Borel cr-algebra of M has compact support. 

If one interprets g as a logarithmic payoff function (eg g = K — cxp, = 1 in case 
of a vanilla one-dimensional put) and P as a Markov chain that models the stochastic 
evolution of the logarithmic prices of assets in a given portfolio, then the set P*, for 
t £ I consists of all those vectors of logarithmic start prices where the probability of 
exercising the option at time t is strictly positive. 

Lemma 7.6. Suppose there is a 70 > 1 (without loss of generality, 70 S {l,e^]) such 
that 

Ptf>Vf 

for all t € (0, r]n/ (where I = HNq with h > whence it is sufficient that this estimate 
holds for t = h). Assume furthermore that maxjgji a;|'''' < 0, implying that g > 
on the set {Ptig V 0) = Ptg}. Then for all ei > there is an Eq < T independent 
of h < T such that for all t G 2 • ((0,eo) CM) and A D LE^ (with positive Lebesgue 
measure), 



sup 

/>svo 



{Birf-BJ 



> 



I,~(A) 



L~(A) 



> 



as well as 



> 



mm a 
ie{i, .••,'"} 



(h) 



e ''2i4:(ln7o-£i)- -, 



sup 
/>svo 



B, 



f-Bsf 



> 



L\A) 



Li(A) 



T 

min a^'^n^ [(Pi (5 V 0) = P* gll ■ e'^'i K (In^io - eA 
ie{i,...,m} J LL 2 2 JJ 



(the left hand side, following the usual convention, being +00 ifX^ |^|Pt (g' V 0) = Pi^l 
+00, m > 1 and e\ < In 70/ 
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Proof. Let us first remark that, due to Corollary 17. II we have 



sup 

/>9V0 



{Bsy{gVO)-Bs{gVO) 



as well as 



sup 

/>svo 



(B^Vf-BJ 



f t 1 • (g-e-''iPt(g\/0)) 



{B.y{gVO)-B,{gVO) 



L~(A) 



J •(<?-e-'-tp.(.gVO)) 



Li(A) 



for all io G |li • ■ • , ?ti2t:| (using the translation invariance of /jj^j •A''). 

Next let us note that by our assumption of a;-^'' < componentwise for all io G 
{l, . . . , TO^^ combined with the north-east connectedness of E'^ (which entails south- 
west connectedness of C£"*), we have 

h " 

Therefore we may conclude that for all t G 2 • / and jq G |lj • ■ ■ , Jriarr |^ 



(7.6) 



> 



> 



> 



mm a^^'x 



ie|l,...,m2*fi I 



lg>e '^2Pt(gvO) 



mm a 
;e{i,...,m} 



{h) 



mm a 

ie{l,...,m} 



(h) 



( > 0) 



2h '° 



5-e-'-5F.(gV0))xc^i 



211 



(7.7) 
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on CE" 

on C^;" 



Now, off one has due to Lemma 17.41 (which may be appUed thanks to our 
assumption max^ xf^"^ < 0) the following situation: 

X{5>e-"P,(3vo)} • (.9 - e^^'^'Psig V 0)) 
= X{P,(3VO)=P,g}n{g>o} ' (s - eT'^" Ps{g V 0)) 
= X{p,(3vo)=p,3}n{g>o} • (g - e""Fs5) 
= X{P,(3V0)=P,3}n{g>0} ■ {K - f- e-'^'K + e-'^'Psf) on CE" 
> X{P.(3V0)=p,3}n{<,>0} • (^-/-e-'^^if + e-"7o'/) on C^;^ (7.8) 

However, one can also perform the calculation 

X{g>o} (K - f^e-^'K + jo'e-^'D 
= X{K-/>o}((i^-/)(l-7o^e''-1+ife— (70^-1)) (7.9) 
> A'e-'-^ (70^ - 1) (7.10) 

(where we have used the assumption 70 < to get from 117. 9|l to Ij7.1()ll l. Combining 
estimates Ij7.1()ll and 117. 811 . we arrive at 

X{s>e-"P,(gvo)} (.9 - e^'^'Psig V 0)) 

> ife-" (70" - 1) on CE' 

> Ke-"^' {\nja - ei) ■ s onZE' 

for every s < Eq for some £0 > dependent on ei > and finally (using estimate Ij7.7ll . 

min,jg{i ,„} a^''^ < 1 and T >t) 



XA ■ e 



'f t 1 • (g-e-''^Pt(g\/0)] 



> 



mm a 

ie{l,...,m} 



(h) 



■ X 



2h 



> 



mm a 
;s{i,...,m} 



(h) 



for alH G 2 • (/n (0,eo)) and iq G |i, . . . ,toA|. 

This yields - due to the translation-invariance of the Lebesgue measure (which gave 
us estimate H7.6|l l - the first line of the Lemma's norm estimate. It also implies 
the L°° norm estimate of the Lemma since for all s G / (in particular for s = |), 



ZE' 



= A'*{P,(gVO)=P,g} 



s . (h) 

— mm x) 



> ^ > 



(a consequence of the monotonicity of g) , and therefore 



mm a 

ie{l,...,Tn} 



{h) 



X 



(fc) • Ke ''2 (In 70 - £1) 



mm a 

ie{l,...,m} 



(h) 



Ke '■2 (In 70 - ei) • -. 



□ 
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Remark 7.3. Assume f is not strictly less than K, say f{z) > K for some z e . 
We can use the property of f being monotonely increasing in each component to see, 
via E,em.ark l7.S\ that 

Ce'^ - {Ph{gVO)=Phg} 

= C{3ze{l,...,m} 5(.-xf))<0} 

m 

= n{/(-!'")s''-} 

m 
i=l 

3 f]lyje{l,...,d} (•-xf)).<z,} 

= f]hje{l,...,d} (•), < + ') . [ 

[Vje{l,...,d} {■)^<z,+ min (xf)) 

'Oo, Zj + mm x) 
ie{i m} V 



where the set in the last line has infinite Lebesgue measure. 
Thus, X'^ [ZE'^] — +00 whenever f < K fails to hold. 

Keeping Corollary 17. II in mind, our next step shall consist in proving 

Lemma 7.7. Let Tel. Suppose there is a 71 > such that 

Pt.f < 71*/ 

for all t E (0,T] n / (where I = KHq with h > and therefore it is sufficient that this 
estimate holds for t — h). Let us define 

D ■■= X(o,en M . , inf / + X[e-,+oo) (71) sup />0. 

Dte(o,T]niE [Jte{0.T]niE* 

Then there is a constant Cq G M given by 

f 1 - e^''' \ / 71*6^'"'' - 1 A 

Cq :— K \ sup ^ + ^'^P ln7i + r 

\se{a,T]ni ^ J Yse(o,T]n/ ^ J 

such that for all s G (0, T] n / and measurable A, 

< ||x{3>e-"n(gvo)} (g-e-'-'PsigVO))]]^,^^^,^^^ 

< X'' [{e'-'g > Ps{g V 0) > Psg} n A] ■ ((In 71 - r) ^ + ri^ + Co) • s. 
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Proof. For all s £ (0,T] n /, the following estimates hold on E'^: 
< X{g>e-"P.(gvo)}(5-e-"Ps(5VO)) 

= X{e"(,>p,(gvo)>p,g} (.9 " e''''' Ps{g V 0)) on E'' 

< X{e-ff>p,(svo)>p,g} {g - e^^'^'Psg) 

< X{e"g>P,(gVO)>P,g} " e^'^^Psg) 

= X{e"ff>P,(gVO)>P,g} {K - J - K + e-""" Psf) 

< X{e.-=,>P.(.vo)>P.,} (i^ (1 - e-'-^) + (7i^e-- - 1) /) 

< X{e-=s>P.(,vo)>P.,} (i^ (1 - ) + (7i^e-- - l) d) 

< X{e"ff>p,(gvo)>p,g} (tKs + (In 71 - r) D ■ s + C ■ 

for some real constant C > that can be bounded by 

/ 1 - e^''* \ ~ / ^i^e"''''-! A 

C<K \ sup r]+D-\ sup (ln7i-r) ^Cq. 

\se(o,T]n/ s J \se{a,T]ni « J 

This gives a uniform pointwise estimate for the nonnegative function X{e'-=g> Ps(gvo)} (5 ~ ^~^'^Psg) 
on -E* from which the Lemma's estimate can be derived immediately. □ 

Corollary 7.2. Let T E I . Assume there exists an iQ G {1, . . . , m} such that 



< 0. 



Then one has 



D = X{o.e-) (71) inf / + Xfe-.+oo) (71) sup / 

and for all s and measurable A, 

||X{g>e-"P.(gvo)} " e'^' Ps{g V 0)) H^ij^^nA) 

< A"^ [{e^'g >Ps{gVO)> P^g} D A] ■ ((In 71 -r)D + rK + Co) ■ s 

Proof. The assumption about io implies that i?* t as t t 00 by Remark f7.2ll . hence 
UtG(o T]ni P* ~ P^ which suffices to prove the Corollary. □ 

Lemma 7.8. // there is an ii G {1, . . . ,m} such that x^l*^ > componentwise (which 
entails {g > 0} C nse(o T]ni {Psig V 0) > 0} by the monotonicity of g in each compo- 
nent), one will have the following upper bound for the measure of the set occuring in 
the preceding Lem,m,n i7.7\i 

X'' [{e'-'g> Ps{gWO)> PsgjnA] 
< [{P,{g V 0) > Psg >0}nA]+X'' [{P,(g V 0) > 0} n {P,g < 0} n A] 

for all measurable A CM."^ . 

Proof. We shall establish an upper bound for the set {e^'^g > Ps{gV 0) > Pg}. Since 
by our assumption 

{g>0}c {F,(.9V0) >0} 
for all s e (0, T] n /, we may, once again exploiting Ps{g V 0) > 0, derive 

{e"5 > Ps{g V 0) > P,.g} n {P,g < 0} 
C {.g>0}n{P,(,gVO) >P,.g}n{P,(7<0} 
c {P,{g V 0) > 0} n {Psg < 0} n {Ps{g V 0) > Psg} 
= {Ps{g y 0)>0}n {Psg <0} 
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This implies 

{e-'g> Ps{9V0)> Ps9} 
C {P,{g V 0) > Ps9 > 0} U {{Ps{g V 0) > 0} n {Ps9 < 0}) . 

□ 

This and Lemma, ITTtI rea.dilv yield, via Corollary [Ol the following 

Lemma 7.9. Suppose T € I and max^x^''^ < componentwise. Assume furthermore 
that there exists a real number 71 > such that 

Phi < 11'' I ■ 

Then for all s £ (0,r]n(2-/), 



Li £;'-nnr=i A+x 



(i) 



< 



, o2 



B.) f-BJ 



it) 



{e'-^/2.g>P,/2(.gV0) >P,/2.g}nA 



e ' 2 



((In 



ji-r)D + rK + Co) ■ -, 



wherein 



and 



D = X(o.e'-) (71) inf / + Xle-.+oo) (71) sup / > 



/ 1 - e-''" ^ 
Cq = K sup r 



sup in 71 + r 

^se(o,T]n/ s ^ 



Proof. Consider t & I. Via our assumption of maxigji ... „j} x\ <Q componentwise, 
one has 



(t) 



t 



(h) ^ rx 

maxxL = — maxx,- < 

" hi' 



componentwise. Since the set is north-east connected, this yields E* + xf^ D E* 
for all k G |l, . . . ,m~|, which in turn - via xe<^ ~ ^e*+x'*^ — 

/c e |l, . . . , I - gives 



t_ 

m h 



fe=i 

XE^Ptf 



for all / > 0. This yields, replacing / by fxA, 

Pt iXE^nAf) = Pt (xb* • XaI) > XiJ* • Pt iXAf) 

> Xe* ■ (■ - 4*0) ^'•^ 



> XEtX 



n.(A+4") 
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for all / > 0. Therefore - using in addition the translation- invariance of Pt and A"* 
(which makes Pt a map that preserves the (A'^)-norm of nonnegative measurable 
functions) - we deduce that for all measurable / > 0, 



I^JII 



< \\Pt{fXE*nA)hH 



\\fXE'nA\\L^{Rd) < ll/llLi(B'nyl) 



From this, using Corolla,rv l7.1l we derive 

,o2 



< 



e-'-^n (x 



Li B^nnJ A+x 



(i) 



(^{,>e-^fp.(,vo)}-(5-^"^^^f(ffV0))) 



L^\E2nnJA+xl 



X) 



g>e "^P, 



.(9VO)} ' 



(<7-e-'-4p.(5V0)) 



'j(^EinAj 



This is enough to prove the Lemma once one takes advantage of Lemma 17.71 and 
Corollary E2 □ 

Remark 7.4. Let the translation-invariant Markov semigroup P be derived from a 
cuhature formula for the Gaussian measure with points {yi, . . . , ?/„i} in such a way that 



a geometric Brownian motion with logarithmic drift /i = ( r — 



ke{i,...,d} 



(r > 



and a € E+ being the interest rate of the price process and the volatility vector, 
respectively) shall be approximated, that is to say 

Vi e {!,..., m}Vj e{l,...,d} (xf))^ =/ifc/i + afc/i^ (y,)fe- 

Then the assumption that all the x^''' be componentwise nonpositive for i £ {1, . . . , m} 
reads 



^ ' /fce{i, 



,d} 



< 



and therefore simply means that fih is componentwise at least as small or even smaller 
than — ( maxi (j/i)^ • (Tfe/i2 ] which, needless to say, equals 

^miui (j/i)^, ■ akh^^ ^ ^ in case of an axis- symmetric cubature formula for the 
Gaussian measure. This assumption is tantamount to 



Vfce{l,...,d} ak^-2h 2 max (y,)fc • cTfe - 2- > 0, 



that is 



\ 



\ 



h 2 max, - h ^ ■ (max^ (yOfe)^ + 2^, 



maxi (yi)^ + ^J /i^i • (max^ + 
^ h-^ ^maxj {yi)^ - y^^niax~(y^)jN^^ , ^ 
f maxi {y.i)f^ + ^J ■ (max, (y^fc)^ + 2r 



for allk ^ {1, . . . ,d} , entailing that P models a basket of logarithmic asset prices whose 
volatilities are bounded below by the positive number 

fmaxi (y,)j, + ^- (max^ (yj)^)^ + 2r 
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Now, emphasising again that our investigations are only concerned with discrete 
translation-invariant Markov chains [Pt)^^i (Markov chains which are derived from 
cubature formulae, for instance), we can use rather elementary inequalities to find 
upper bounds on the subsets of W'- occurring in the estimates of Lemma, ITTzl 

We will start with the simple, nevertheless practically important, example of a 
one-dimensional American vanilla put: 

Lemma 7.10. Suppose d — 1 and f = exp. Under these assumptions there exists a 
7i > such that Ptf ~ li' f , and furthermore, one has for all s G /, 



1 



{e'-'g> Ps{g WO) > PsgjC In K 
Proof. The real number 71 is given by the relation 



min x-'^-'jO 

h ie{l....,m} 



71 



E 



a) 'e 



that is 



71 



Next we observe that on the one hand by Remark [221 

{Ps{gWO)> Psg} = {3fce {l,...,mt} 5(.-4^))<o} 



and secondly 
thus 



mm g [ ■ 

fcej i,...,mT I 



K — max exp 

feej l,...,mt| 



K — exp 



InK <■ 



mm 
fce( l,...,mt 



x['^ 1 < 



4'^ I < 



mm X 



is) 



fcej i,...,m h I 



mm xj^ , +00 
feej i,...,mT I 

= \ \nK + ^ ■ min x.-'^^' , +00 

h ie{l,...,m} , 



> 0} = {X > exp} (-00, In i^) , 



{e"g > (g V 0) > P,.g} C \\nK+^-- min xf^ , In K 

h ze{i,...,m} 



In A' 



1 



h je{i,...,m} 



min X , ) • s 



□ 
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Applying the preceding Lemmas and using Corollary 17. II we conclude by stating 

Theorem 7.1. Suppose d — 1 and f — exp. Under these assumptions there is a 
7i > such that Ptf — 71*/ for all t £ I . Assume, moreover, that 



Vze {1,. 



i}Vje{l,...,d} 



< 0. 



Then there is a real constant D such that for alls e (0,r]n(2-/) and for all f > .gVO, 



We can compute D explicitly as 



< 



B, 



f-Bsf 



D 2 



D 



{\n-fi-r)b + rK + Cq 



mini X 



(h) 



Proof. One only has to apply Lemma FTfll for A 
Lemma Oni 



which one is entitled to by 

□ 



Now we shall proceed to establish convergence estimates for the sequence (i?T-2-"/)ngp 
in the L^{E^ n ^)-norm, for all measurable f > gV and measurable A C R''. 

Lemma 7.11. Suppose d = 1 and f = exp. Under these assumptions there is a 71 
such that Ptf = 7i*/- Assume, moreover, that 



yie{i,...,m}yj e{i,...,d} (xf')^ 



< 0. 



Under these assumptions there exists a real number D > (the same as in Theorem 
E2P such that for all keNo, se (0, T] n (2'=+i • /) and measurable f > gW 0, one has 



Li(Ai[-B''n-]) 



The proof is contrived inductively, the base step being Theorem l7.lL and the induc- 
tion step being the first part of Lemma f7.1^il However, the second and more general 
part of Lemma I7.1^il - which we will need later on in this Chapter when we study 
options on multiple assets - requires the following auxiliary result. 

Lemma 7.12. Let t £ I, A CM.'^ measurable, and assume 

max a;^'') < 

ie{l,...,m} 



(which due to I = HNq is equivalent to max^^j^ x*^"^ < for all s £ I). Th 

Clseio T]ni = ^'^ Remark\^^ and for all /i > /o > .9 V and p £ {I, +00}, 
\\BJi ^ BtfoW 



en 



(A+.«)n. 



< e" 



/i /oIIlp(A''[£;'' 



nAn-]) 



Proof. Consider a measurable set A C M'' and measurable functions /o,/i > V 0. 
Similarly to the proof of Lemma 17.21 we observe that due to the monotonicity of Bt 
and the fact that B*/ > .g V > g for all / > 0, 

{Btfi=g} = {Btfo < Btfi = g} 

= {g < Btfo < Btf, = g} 

= {Btfi^g}n{Btfo^g} 

C {Btfi-Btfo = 0}, 
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that is 

{Bth - Btfo ^ 0} C {Btf, ^g} = {Bth > 9} 

Combining this with the monotonicity of Pt as well as the fact that E* — xf^ C for 
all i (which in turn is a consequence of the north-east connectedness of - cf Remark 
17.21 - and the assumption that < for all i), we obtain 



< {Btfi - 5t/o)X£.nn,(A+4'>) 
= (St/i -BJo)x^.nn,(A+4'))n{B,/i><,} 



-rt 



-rt 



e 



2:»."'*..nn.(.M")(^-(-^'")-^»(-^-")) 

t 

t 



-rt 



Now, since 



i,...,m7r I 

for alH G |l, . . . , mi | and /i — /o > 0, this means 

t_ 

m h 



E "^'^ (^(£'^-4') nA - /o)J (• - ^'t^) 



Combining this pointwise estimate with the translation-invariance of the Lebesgue 
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measure yields 



< e 



/ {Bth-Btfo)x 

t 

t 



i=l 

t 



(/1-/0)) (•-xf))rfA'^ 



< e 



< 



[ e-^*XE''nA{fi-fo)dX'' 

JR'' 

= e"*^* ll/i - /ollii(^hn^) , 

where we have used the inclusion E"^ — .Tj*-' C E''^ which - owing to the north-east 
connectedness of the sets E"^ and our assumption maxjg{i_ „j} x'-^^ < - holds for 

arbitrary k e |l, . . . , | and s,t € I as well as the assumption /i — /o > 0. Similarly, 
the translation- invariance and the sub-Hnearity of the ess sup Rd-norm imply 

ll-^*-^i-^*-^oliL~(A'^[B''nn,(A+xW)n.]) 

(Bt/i-Bt/o)x^.nn,(A+x«) 



= ess sup Rd 



< ess sup 



-rtS^ (t) 



E 



(/1-/0)) {--4') 



(t) 



ess sup Rd 



t 



(A - /o) 



„(t) 



< 



(t) 

'ess sup 



(/i - /o) 



i=l 

i 

< e"^^ ^ af ^ess sup [xe^ha ih - fo)] 

i=l 

t_ 

m h 

= e"''* E 'ess sup Ef^nA (/i - /o) 
= e"''*ess sup e^ha (/i - fo) 



where again one has exploited the inclusion — a;^*^ C E^ that holds for any k € 
|l, . . . ,m^| and s,t G I. 

□ 

Lemma 7.13. Let T G I and p € {l,-|-oo}. Consider a real number D' > and a 
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measurable set C CM.'^. Suppose one has an estimate of the kind 



v/>5V0Vse (2-7)n(o,r) 



D' 



LPi\''[Cn-]) 2 



Assume, moreover, niaxjg|i a;*^''-' < (which by the Chapman- Kolmogorov equa- 
tion is firstly equivalent to a;^*^ < for all k G |l, . . . ,m"s^| and t € I and secondly 
also entails nte(o T]ni ~ E^). Then we get for all measurable f > g\/0 and for all 
fc e No, s > such that s e (0, T) n (2'=+^ • /), the estimate 



LP(A<i[Cn-]) 



< D'-s^. 2- 
Furthermore, if one assumes in addition 

e ^ : i e 



{l,...,m}|, 



then one has a related implication for j fl Q 



He 



{l,...,mi } 



A + x. 



(f) 



instead of 



LP (C) for all measurable Ad M"; // under these assumptions the assertion 
V/>ffVOVse (2-7)n(0,T) 



B.y'f-BJ 



D' 2 
< — ■ s 

Lp(E''nni(A+x<''')) ~ 2 



holds, then the estimate 



V/>5V0 



LP I E' 



'"nni(-4+a:<°')) 



holds for all fc € No and s>0 such that s e (0, T) n (2'=+i • /) . 

Proof. For both parts of the Lemma, we will conduct an induction in fc G No , the initial 
(or base) step being tautological each time. We have for all s € (0, T) fl (2*^+^ • /) and 
f > gV the estimate 

(B,.2-(.+i)r(''^')/-(B«.2-o°'V 
= (B,.2_(.+i)r(^')o(B,.2-(.+i,)°('')/ 



(s,.2-(.-l)) 



>(2'=-i) 



Bs 



0(2'=) 



.2-('=-i) 



3(2'=-i) 



o(B«.2-0 



0(2'=-!) 



/ 



which plays a crucial part in both the first and the second part of the Lemma. For, 
we can first of all note that the induction hypothesis in the situation of the first part 
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of the Lemma reads 



yf>g\/ ovt e (2*= • /) n (0, T) 



LP{C) 



< D' -t^ ■ 2-''. 



(7.11) 



And if one now applies this induction hypothesis Ij7.11ll for t = f (recaUing that by 
assumption s e 2'''+^ •/, thus | € 2'^ •/) to the previous two equations and uses Lemma 
.121 then one gets by the triangle inequality for the (C)-norm, 



< 



.2-k 



+ 



< 



{Bg.2-{k + l) 

0(2") 



{Bs.2-(i' + ^))° ' f - {Bs-2-*')° f 

0(2 



LP iO 
(2'=-i) 



Bs 



(2*=) 



,(2'=) 
.(2'=-i)^ 



B 



.2-('=-i) 



LP{C) 
f 



LP(C) 



B. 



y''yiB^.2-^)<'-)f 
0(2"-!) 



LP{C) 



.2-('=-i) 



< D' ■ — ■2~'' + D' ■ — ■2-'' = D' ■ ■ 2-^^+1). 



LP(C) 



In order to be entitled to apply Lemma r7.12l in this situation we have successively used 
the fact that 

yt e /v^ > .g V Btt>gy 0. 



This completes the induction step for the first part of the Lemma. 

Turning to the proof of the second assertion in the Lemma (where maxigji x^^^ < 
is assumed), we remark that 



Vt e /Vfc < 
Ait) fl 

ie| l,...,mT I 

- n 



In < — In 



In 2 



A + x 



fcG| l,...,m7r 



n 



ii,...,ir^k'^<, l,...,m7r - 



A 



(t.2-('=+i)) 



(7.12) 
.(7.13) 



In particular, if G ^xl'^^ : i e {1, . . . , ^(s) is decreasing in s: 

ys,t e I {s <t^ A{s) D A{t)) . 
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Similarly to proof of the first part of the present Lemma, we deduce 



.(2'= + i) 



LPiEt'nAis)) 



< 



< 



0(2") 



+ 



(^f-2-fc) - (-B|.2-('=-i) 



LP{E'^nA(s)) 



{B^.,-,.^.,r^ ' o IBs 



.2-k 



- B 



°(Ss.2-r^ >f 
(2-^) 



.2-('=-i) 



LP(£;''nA(s)) 
/ 

Lp(_E''nA(2'=+if )) 



< 



^s..2-<'=-)) o(Ss.2-r' V 



LP{E''nA(s)) 



Lp(£;''n/i(s)) 



from the triangle inequality. But by a successive application of Lemma [7.1 21 combined 
with the properties II7.13|I of A{-), we have for all /i > /o > V 0, 



< 



< 



< 



(B..2-(^+i))°^'')o(/i-/o) 

(i?,.2_(.+ „)°(^'")o(/i-/o) 
(i?,.2_(..„r(^')o(/i-/o) 

(B..2-('=+i)r^''"')°(/i-/o) 



Lp(£;''nA(s)) 
LP I _E''nn 



e l,.-.,»n2h 



A(f)+x 



LP E'^nf] 



2fc+l 
-1 



< ll/i - /ollLp(£;'>nA(f)) • 
In light of the inclusion A{s) D A (|), we also have 

{Bg.2-{k + l))°^^ ^ f - {Bs.2-k)°'^ f 



LP(E''nA{s)) 



< 



02*= 



I O {B,.2-k) f 



+ 



B^. 

Bs,2-(k-l) 

0(2") 



LP(Ei^nA(s)) 



.2-k 



-[Bs2-,k-^A ']J 



Lp{EfnA{^)) 
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Combining the previous two sets of estimates leads to 

(B,.2-(.+i,r('')o(/i-/o) 
/ 



LP(E''nA{s)) 



< 



(B4.2-.) 



02" 



(B. 



s-2-'' ) 



-Bs 2-(fc-l) 



.2-('=-i) 



3(2'=-i) 



Lv[E<^nA{^)) 
f 



Lp{E'^nA{i)) 



< D' ■ — ■2-'' + D' ■ — ■2-'' = D' ■ ■ 2-('=+i\ 
- 4 4 



where in the last line we have taken advantage of the induction hypothesis 

Vfc e NoV/ > 5 V OVi > 



V 



te{2'' -1)0 (0,T) 
< D' ■ 2-^ 



(2*=*) 



for the special case t — 



□ 



ih) 



< 



The assumption of € ^x['^'' '■ «€{!,..., C M"* while maxjgji ... „j} a;^ 
componentwise corresponds to the volatility attaining a certain critical value: 

Remark 7.5. Consider a cubature formula for the one- dimensional Gaussian measure 
with cubature points {zi, . . . , z^} which will then give rise to a new Markov chain via 



V* e {1, 



. , m 



} 



(if simply {zi,Z2} = {il}; then this was a discrete model for a logarithmic asset 
price evolution that converges weakly to the Black-Scholes model with volatility gq and 
discount rate r > v)hen h [O). In this setting, the set of pairs (r, (Tq) G such that 







has at most two elements, ie it is a Lebesgue null set. However, in practice, we will 
not have the exact values of the volatility ctq ( and if the maturity is sufficiently large, 
one will not even have an exact value for the interest rate r), but we will only know 
that (To e ((T ~ £, (7 + e) for some £ > 0. So, given r > the set of volatility parameters 
that both fit the model and allow for the previous Lemma to be applied will equal 



(To € (a — £, (7 + e) 



o-o 



h + (To/i2 • maxZj 







If e > and the equation characterising this set has a solution cr' e {a — e,a + s), this 
set will at least have positive Lebesgue measure, so that there is some hope that our 
condition of max^ xf^"^ = (which we had to impose in the second part of the previous 
Lem.m,n i7.1Sp can be satisfied in practice at least occasionally. 

With the first half of Lemma fZ.l.SI we have completed the proof of Lemma !?. Ill We 
shall now apply this result to finally get to a convergence bound for {Bx.2-p (g V 0))„ 
- which can be conceived of as a sequence of non-perpetual Bermudan option prices 
when successively halving the exercise mesh size. 
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Lemma 7.14. Let p G [1, +00]. Consider a real constant D > as well as a measur- 
able set C and a set E of nonnegative measurable functions, and suppose one has an 
estimate of the kind 



Vfc e NoV/ e EVs e (2 



fc+i 



/)n(o,T) 



{B. 



LP(A''[Cn-]) 



<i?.s2.2-('=+i). 



Then for all N > M eN, s e (0, T) n (2^ • /) and f £ E, the estimate 



o(2«) 



LP(C) 



-M 



-M 



1 - 2"(^"*^"i) 



holds. 

Proof. With M, TV, s, / as in the statement of the Theorem, we obtain by the triangle 
inequality 



(S,.2-«) 



oi2" 



.0(2-) 



LP(C) 



< 



c=Af 



s.2-(fc + i) j ^ W - [Bs.2->') ^ ' 



k=M 
N-1 



LP{C) 



LP{C) 



k=M 



D.-.2 



< D 2" 

2 



A/ 



• 2. 



1-2-1 

This suffices to prove the Theorem. 

Thus, if we combine this last Lemma FZ. 141 with Lemma FZ. Ill we arrive at 



□ 



Theorem 7.2. Suppose, as before, d — 1 and f = exp. Under these assumptions there 
is a 71 such that Ptf — 71*/; '"^'^ let us suppose this 71 G (0,e'']. Assume, moreover, 
that 

Vie{l,...,m}Vje{l,...,d} (xf')^<0. 

Under these assumptions there exists a real number D > such that for all N > M G 
N, s € (0, T] n (2^ • /) and monotonely decreasing f > gV 0, one has 



< Z) • • 2 



< D-s^-2-^'. 



Analogously, we may proceed to prove convergence of higher order in s for / ~ 
J2j=i '^j 6xp {{■)]), where wi, . . . , is a convex combination (the weights for a weighted 
average of the components / assets in a d-dimensional basket) , as well as for the choices 
/ = min^g^i £;}exp((-)j) and / = maxj-gj^^ ,i}exp((-)j). However, this time, we 
shall employ different norms: [X'^ [i?'' n ^ n •]) for a compact subset A C M'' such 

that A'' nse(o.T]n/ ET\ A e (0, +00). 

The first part of this endeavour will be to prove certain generalisations of Lemmas 
O and 1232 
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Lemma 7.15. If f = Wj exp the 



whe 



71 := 



max 'S^ a^l^\ 



V 

Proof. We have for all s € / the estimate 



>o 



Psf = 



m h 
i=l 

7=1 i=l ^ ■'^ 



< 



j=i i=i 



max > 
.•e{i,...,d}^ 



' 3 Wj exp I (• 



>o 



in particular this estimate holds for s = h. But this is to say 

hence we have proven the estimate in the Lemma for s = h. This readily sufSces to 
prove the Lemma's assertion in its full generality, as {Ps)g^j is a Markov semigroup 
and by applying the Chapman-Komogorov equation (and the monotonicity of Ph) 
inductively, 

Vn e N PnJ =Ph---Ph!<li''--- li f = 7i''"/" 



□ 



Lemma 7.16. // / = min^g^i^ i^-j. exp then 

Psf<ll'f 

where 



( 



71 := 



E 

i=l 



a. 



is) 



max 
?e{i,...,d} 

V 

>o 
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Proof. We have for all s G / the estimate 



Psf = E"i^V-(--^rO 



i=l 



E(s) 
al _ mm ^ exp 



i=l 

m'h 



je{i,...,d} 



< 



•e^""^' )*exp((-)^) 



i=l 



'Ifee-fl d> • "'"'je-[i,...,d}«p(( )3) =«p(( )fe). "|(^) 

\ Vn > mln^g{i^ ,jj exp ((Oj) < exp ((•)„) J 



< 



E 

i=l 



max 
^e{i,.--,rf} 



mm exp 

3e{i,...,d} 



(<•).) 



>o 



in particular the estimate holds for s = h again. But this means 

hence we arrive at the estimate of the Lemma for s = h. This readily suffices to prove 
the Lemma's assertion in its full strength, as (Ps)j,g/ is a Markov semigroup and by 
applying the Chapman-Komogorov equation inductively, 



Lemma 7.17. If this time f = max^g^i^ rfj exp then 

Pel < I 

where 



( 



71 



Oif' max 



d} 



V 

Proof. We have for all s e / the estimate 



>0 



i=l 

= y^a-*^ max exp ( ( • 

^ ^•6{l,...,rf} 



< 



E 

i=l 



a,- max 



max exp ((Oj) 



i6{i,...,d} 



a; ' m 



max 

d} 



' f max exp 

j 



(a) 



□ 



>o 
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in particular the estimate holds for s — h. But this is - as it was in the proofs of the 
two preceding Lemmas - to say 

A/<7i'/" 

hence we have proven the estimate in the Lemma for s = h. This readily sufBces to 
prove the Lemma's assertion, as {Ps)s^j is a Markov semigroup and by applying the 
Chapman-Komogorov equation inductively, 

Vn e N PnJ =Ph--- Ph / < 71'' • • • 71" / = 7i''7" 



□ 

Thus at least for certain choices of / - viz. weighted arithmetic average of the 
exponential components, minimum of the exponential components and maximum of 
the exponential components - we can apply Lemma FZ. 71 

Therefore we shall next turn our attention to deriving upper bounds for the mea- 
sures of the sets in the estimates of Lemma, l7.7l for the said examples of / = ^xp {{■)]), 
f = min^gji (;}exp((-)j) and / = maxjg{i ,^}exp((-)j). We continue to use the 
notation / = hNo and 

m'h 

i=i 

where {Ps)s(zi = (Pnh)neNo = I ^'^ "J j '^^ ^'^^ Markov chain generated by Pu- 

\ n / neNo 

Lemma 7.18. If f = niinjgji exp then for all s G I , 



InK -\ max fa;,-''') ,+00 

{l,...,m} V / 



s 

h ie{ 



as well as 



|vi e {l,...,m*} .9(--a;|''^) >0} 

max \ X; ] 
yi- . ' V h^e{l,...,m}\ ' Jj^ 
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Proof. Let s G /. Then 

{Vie{l,...,m*} 5(--arf))<0} 

= < max g{■—x^^^\<Q\ 

1 ie{l,...,mt| ^ ^ J 

= Ik min min exp ( ( • — xf^] I < > 

I .G{i,...,m*}ie{i.-..,d} hj J 

= I K — min min exp ( ( • — x^^M | < > 

= <-?!'— min exp I (•)■ — max (x^"-' | 1 < > 
I ie{i,...,d} \^'^ ie{i,...,mi] ^ ^ij ~ J 

= IhiK < min I (■),— max (a;-"'') > 
I - ie{i,...,<i} ie{i,...,mt} ^ J 
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and also 



{Vie{l,...,m} g(--x\'^)>0} 
K — max min exp ( ( • — x^^^] ] > > 

^e{l,...,mi}i^{^'■■■'d} \\ JjJ J 



c 



K — min min cxp 

»e{i,...,m*} 

K — min min exp 
j^{h-^d} ie{i,...,mi} 



x^/'> 1 I > 



> 



K— min exp | (•),- — max [xl^^] \ 
\^ ' ie{i,...,mi} ^ ^i) 

\aK > min (•),■ — max (a;-^M > 
= M i lni^+ max (xf^) > (^j 



> 



d 



l,...,m I 

Milni^+J max (xfA > (•) 

^ 1 /l ie{l,...,m} \ ' Jj 



Corollary 7.3. If f = min^g^i^ exp then for all s G I, 

{Ps{g V0)>0} 



= IJ Mx---xMx( -oo, InK +^ max (x'^/'A 
I ' » ' V ^ ie{l,...,m} V * J j 



X M X • • • X 



^■=1 \ i-1 

Proof. Let s G /. We simply remark that 

{P,(ffVO)>0} 
= {3iG {l,...,m*} ,9 (• -.xl'^) > O} 

= C{vi e {l,...,m^} 5 (• - x,^'^) < o} 
Lemma 7.19. If f = max^gji cxp ((Oj), ^'^e'^ for all s € /, 

{vie{i,...,mt} <o} 

lnK+^ min fa;,^'*^') , +00 | x K 
/i ie{i, .••,"»} ^ h J ^ 



C y R X ■ ■ ■ X R X 
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as well as 



-ooAnK + — mill 

h ie{i,...,m} V * /j 



Proo/. Let s e /. Then 

{Vze {l,...,m} g(.-xf))<o} 

max g (• — x[^^j < 
ie{i,...,mt} ^ ^ 



c 



ii' — min max exp f ( • — a;,-"'' ) ) < 

K — max max exp f (• — x,-*'' ) ) < 
,e{i,...,mt}ie{i,...,d} h) 

K — max max exp ( ( • — x-*'' ) | < 

K — max exp ( (■)■ — min fxl''''' | 
\^ ' ie{i,...,mt} ^ 

IniiT < max (•),• — min (x-"') > 
je{i,...,d}' ie{i,...,mi} ^ h\ 

\\\\nK+ min (x^"') < 
y=i [ ie{i,...,mt} ^ 

y.{'°^-+L.,=„,(»s"),^()4' 
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and also 



{Vi€ {l,...,m} g(- 



min 5 ( • — x^^^ ) > 
iej i,...,m's' I 

if— max max exp ^^• 



K — max max exp 
je{i,-,<i},g|i,...,™t} 



((■ 



x^/A I >o 



x.^'M 1 > 



jell,. ..4} 



iei l,...,m h 



K — max exp I (•) ■ — min [x\^^ ) I > 



K > exp max 
- ^ ' je{i,...,d} 1 



mm I 

ie{l,.../mt} ^ 'J 



In if > max (■)■ 
- je{i....M} 1 



mm I xl I 

ie| l,...,m"E I ^ 



= n^^^+ min (xi^)) >(.), 

j=l [ te[l,...,-mi] ^ ' ^ 



n Inif + J min (xf^\ > {■), 

11] /l ie{l,...,m} \Jj 

3 = 1 



Corollary 7.4. If f = maxjg{i_ exp then for all s G I , 

{P,(5V0) > P,g} 

= \ \ \ Rx ■■■ xRx (hiK + ^ min ( x^,^'^] , 

^ . ' V h^e{l,...,m}\ ' Jj' 

J = l \ i-1 

Proof. Let s G i. We simply remark that 

{F,(<?V0)>P,<7} 
= {3ze {l,...,mt} <o} 

= C{vi e {l,...,m*} 5 (• - xl'^) > o} 




□ 



□ 



These estimates lead to the following Corollary that will enable us - under the 
assumption of 

max a;, =0 

ie{l,...,m} 

(in order to be entitled to apply eg Lemma r7.13|l - to prove an L^-convergence estimate 
(on a particular subset of R'^) for (i?T-2-"/)„gN ^'^^ measurable / > .9 V 0. 
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Corollary 7.5. Suppose f — max^gji ,i].exp((-)j) and consider any compact set 
B CR"^. Then 



1 / 

ft, ^ V ie{l,...,ni} 
.7 = 1 



vo 



for all s e /. 

Proof. Let s G /. Since 

{e''^5 > P,(g V 0) > P.g} C {P,(5 V 0) > P.g} n {g > 0} 
by the monotonicity of P, , we only have to observe that 



{ff >0} 



ii' > max exp 
je{i,...,d} 



= n{^'> ((•).)} 

d 

= (g)(-cx),lni^) 

to arrive - after taking advantage of the preceding Corollarv 17.41 - at 

{e"g>P,(.gVO) >P,.g} 

/ (— oo, hiX) X • • • X (— oo, Ini^) \ 



d 

^ u 

3 = 1 



(h) 



X l\nK + f min,,g{i^...^„} | x. 

X (— oo, hiif ) X • • • X (— oo, \nK) 



,\nK 



d-3 



d 

u 



S<o X • • • X K< 



^ \n- ie{i,...,r?x} V / j 



• s X M<o X • • • X K<o 

d-j 



+ (lni^),li. 



However, by our assumption that B be compact, there is some P > such that 
P- (In/O^^i C [-R,Bf. Thus 

{e''^g>P,(5V0)>P,.g}nP 



c y(i-«-''r'>=(l,,r„,(-S''')/'')^ 



s X 



and from this inclusion we may deduce the estimate given in the Lemma. 



□ 



The inequality we have just derived impHes that the 



E' 



n pl^ {^B + xf^^ n • -volume of the set occurring in Lemma [7.71 is of order s 

for any compact B and for / = max^gji (;}exp((-)j). Hence again by Lemma [7.71 
(which is applicable because of Lemma l7.17|l we obtain that the difference Bs/2f — Bgf 

is of order (this time, however in the (^E'^ n f]g(B + x^^^YVnorm) . This esti- 
mate on the norm of Ps/2/ — B^f leads, via Lemmas 17.141 and 17. 131 to the result that 
the analogon of the difference in Theorem 17.21 is of order • 2~^^, too: 
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Theorem 7.3. Suppose f = laaxj^^i^ ^^ exp and consider a compact set B. 

Assume that 



max x] = 0. 



ie{l,...,m} 

Under these assumptions there exists a real number D > such that for all N > M G 
N, s e (0, T) n (2^ • /) and f>gVO, one has 



(i3,.2-«)°('")/-(i3,.2-M) 



>{2") 



(l _ 2-(JV-M-i)^ 



Remark 7.6. This -convergence result has some (however, because of our assump- 
tion maxigji ,„} xf^"^ = 0, fairly limited) practical interest, as in practice quite fre- 
quently the exact start price of the (multiple) asset on which an option is issued, is 
unknown. Instead, one will have the logarithmic start price vector x gW'' a short time 
S > before the actual option contract becomes valid. Now, asuming that P^^ has a 

continuous density this function will be bounded on fl ( -B + xp J by 
some constant 



C :-- 



sup 



< +00. 



B''nn^(.B+x^''') 

One will therefore have for all f > g V 0, s e {0,T] n I and N > M gN, 



(s,.2-«)°('")/-(i?..2-My 



,(2-) 



< 



(B,.2-«r('")/-(s,.2-M) 



.0(2") 



< C- 



.o(2«) 



/ - (-Bs-2-m) 



0(2") 



Li(A<i[i=;''nn,(B+x^''')n-]) 



< C-D-s^-2- 



-M 
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Part IV 

Numerical analysis of 
cubature-based American 
pricing 
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Chapter 8 



Motivating Bermudan pricing 
based on cubature 

8.1 The general setting 

This Chapter has been designed to elaborate the idea of Bermudan option pricing via 
cubature and to put it into a context of other Bermudan option prcing algorithms. 
Given its informal character, it can also be perceived as another introductory chapter. 

Consider a basket of d assets. A d-dimensional Bermudan option is an option that 
can be exercised at a discrete set of exercise times, yielding payoff g{x) if a; e M'' 
is the vector of logarithmic prices (of stocks in the d-dimensional basket) at that 
time, for some g : ^ M. which will be called the payoff function. In the case 
of a one-dimensional put option with strike price K for example, one would have 
g — {K — exp)+ — [K — exp) V 0. In case of a call on a stock index, 5 : R*^ ^ R would 
be the positive part of the difference between a weighted sum of exponential functions 
of the coordinate entries and the strike price. Unless specified otherwise, we will from 
now on assume the exercise times to be equidistant with an exercise mesh size h> Q. 
Adopting the notation of Chapter this is to say J = /iNq. 

We regard such a Bermudan option as a binary tree of European options. This 
means that at each exercise time one has to decide whether it is more rewarding to 
keep the option or to exercise it - in other words, whether the payoff at that exercise 
time is less than the value of the (European) option to exercise at the next exercise 
time. A recursive algorithm is thus obtained. Note however that this binary tree of 
European options has continuum many nodes at each level, one for each price vector 
at the subsequent exercise time. 

Let us now describe this recursion in detail. Suppose the option is non-perpetual, 
i.e. it has a maturity time T £ (0,oo), and assume furthermore that T = hN for 
some TV e N. Then the Bermudan option price for a start price vector G M'' 
will be Vo(a:o) where the Vfc, /c € {0, . . . , N} are computed according to the following 
backward recursion (where for simplicity we assume the logarithmic discount rate to 
be a constant r > and (It)* denotes the M+'^-valued process of vectors comprised of 
the asset prices in the basket): 

Vn = g, 

Vfce {1,...,7V} VxeR'' Vk-i{x)^ina^{g{x),E^e-''''Vk{Yh.k)}. (8.1) 

Recursion formulae of this kind for the pricing of Bermudan options are fairly 
standard and can be found for example in textbooks such as Hull's jl^^j or Wilmott, 
Howison and Dewynne's jHH]. To use these recursion formulae practically, one needs a 
way to summarise or approximate the state Vk in a way that permits the equivalent 
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summarisation or approximation for Vk-i (this for all positive integer k). One method 
of achieving this in a one-dimensional setting is the application of Fourier-Hermite 
expansions to the functions Vk, as studied in the paper by Chiarella, el-Hassan and 
Kucera (Hj- Our goal is to develop their approach; in particular we will employ cuba- 
ture formulae for symmetric measures. These methods of approximating integrals by 
weighted (finite) averages can be computationally efficient, and with increasing dimen- 
sion may be superior to other approaches. Victoir introduced a vital improvement 
by constructing sequences that scale well. This route to high-dimensional Bermudan 
and American option pricing was proposed for the first time by my supervisor |2(1] . 



8.2 Application to the Black-Scholes model 

Let us in this section work within the situation of the multi-dimensional Black-Scholes 
model, that is to say that the logarithmic price processes of the d assets in the basket 
are independent Brownian motions with drift. Let us assume the volatilities of the 
assets to be constants at for i e {1, . . . ,d}. 
We set 

for i e {1, . . . , c?} and denote by i^a^a^ for a € M'', cr > the d-dimensional Gaussian 
probability measure of variance ct^ centered at a. We will assume that - possibly 
after an appropriate change of the time scale and the discount rate r (by a linear 
transformation from the left) - we have ai = 1/Vd for alH S {1, . . . , d}. /i will change 
according to its definition. 

Then Ito's Lemma implies that the logarithmic price process in the z-th coordinate 
is - with repect to the risk-neutral measure - just a Brownian motion with drift fit 
and volatility l/Vd, thus the process of logarithmic prices of the assets in the basket 
is a Brownian motion with drift /i and volatility 1. Therefore our recursion fomula 
118. l|l becomes 

Vfc e N n [0, T] e R'^ 



max 



-rh 



Vk{^ + ^,h + h'^^■)dlyo,l,g{0 



Now, if the points S,j £ with respective weights aj > 0, j G {1, . . . , m}, determine 
cubature formulae for the standard Gaussian measure vo.i, we can approximate the 
previous recursion by the following formula: 

Vfc e N n [0, T]V77 e M.'^ 

8.3 Exploiting combinatorial aspects of Gaussian cu- 
bature 

Thanks to the work of Nicolas Victoir (which has later been extended by Christian 
Litterer), there are "cubature formulae with few points" [32^ for the integration of 
polynomials with respect to the standard Gaussian measure up to a certain degree. 
Although "asymmetric" (22], their shape is quite regular and uniform. Since the re- 
cursion following the previous recursion formula amounts to the evaluation of payoff 
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functions at (modified) sums of these cubature points and we therefore desire recom- 
bination of these sums, this will turn out to be a computationally palpable advantage. 

The commutativity of (R"*, +) and the equidistance of the exercise times already 
enable us to perform a geometric argument based on the regular and uniform shape 
of the cubature points, which results in 

Theorem 8.1. Let d = 3fc — 2 for some fc g N. The recursion according to \H.^) . 
using the cubature formula for the integration of degree 5 polynomials with respect to 
a standard Gaussian measure from Victoir's example ]3'iA 5.1.1], is polynomial in j^. 

Proof. The cubature points of the cubature formulae referred to in the Theorem form 
a finite subset of \/3{0,±l}'^. Sums of length (provided this fraction is an integer) 
of the cubature points are therefore always elements of (Z'^ n {|-| < -^j) (and this 
set has only (2-^)'' elements), and the points used in the recusion formula stated above 
are comprised of a subset of /i^/^ • v^S (Z'' n { H < i}) + • {O, . . . , i} +Co- □ 

However, this is not the only recombination that can be accomplished in the case 
where d — 3k ~ 2 : 

Remark 8.1. Let us look at the tree obtained from starting at some point £,q G M'' 
and then at each node letting exactly \{xo} U v^S • tJaXi] branches leave (where {xo}Li 
• G3X1 in Victoir's notation is the set of cubature points he uses in the example 
/.Viil 5.1.1] we are referring to), exactly one branch for each element of the set ^0 + 
{{xoyu^-GsXi). 

If we intend to find and eliminate the branches of the tree that are computed "waste- 
fully", it is reasonable to divide the sums (of length j^) of the cubature points by and 
consider them coordinate-wise modulo 2. Then one is looking at elements of the vector 
space (Z/2Z)'^. For the sake of simplicity, let d = 7, that is k — 3 in the notation of 
the previous Theorem \8.1[ The coordinate-wise projection of the -multiple of our 

set of cubature points {xq} U \/3 ■ G3X1, where 

Xi := . . .,xij} 

and 

xo:=0, 03 = ({±1}',*) 

(in the notation of JSSl 5.1.1]), into the vector space (Z/2Z)^ now contains only eight 
points (instead of 57 as before). 

Thus, using basic linear algebra in a 7 -dimensional 2, /2'L-vector space, we are eas- 
ily able to classify the non-trivial zero representations from elements of the projected 
cubature points. 

Perceiving Xi as a 7 -element subset o/(Z/2Z)^, we see that {x)xeXi is an invertible 
{Z/2Zy^'^ -matrix. Therefore we cannot expect any recombination from representations 
of zero by nontrivial linear combinations of elements of Xi C (Z/2Z)^. Moreover, the 
fact that A := {x)xeXi is invertible, shows that xq — can only be written trivially as 
a sum of elements of Xi . Hence we have shown that we exploit symmetries optimally 
if we use: (i) the commutativity of (M'',-|-); (ii) the obvious symmetries due to the 
construction of the cubature formulae by means of the action of a reflection group on 
certain points; (Hi) the fact that addition of xq does nothing at all. 
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Chapter 9 



Numerical results 



In this Chapter we shall present some numerical results. We have decided to choose a 7- 
dimensional example, since (1) most previous research has stopped short of numerically 
tackling American options on baskets with more than 5 assets, (2) it is the smallest 
dimension d > 5 in which Victoir's cubature formulae for the normal Gaussian measure 
of dimensions rf = 3fc — 2 (where k e No) Example 5.1.1] hold. 

We shall assume that the basket X as a logarithmic price process follows the Black- 
Scholes model for independent assets with discount rate r > and volatilities > 0, 
i £ {1, . . . , d}, that is 



(where B is the d-dimenional Wiener process). Given a payoff function, a strike price, 
a maturity T > 0, and logarithmic start price vector x S K"* we shall vary the exercise 
mesh size h > (say h G {ho, . . . , /i„}) and compute approximate Bermudan prices 
C/™°(T)(x). Then we will extrapolate the function h ^ U™°{T){x) to /i = by 
assuming assume that h i-^ [/''^o (T) (x) is a polynomial of degree n in h" for a given 
a > that finally shall be varied as well. 

Unfortunately, it is difiicult to find data on American option prices for dimension 
d = 7. However, one can of course use our algorithm sub-optimally for d = 5 through 
letting the payoff function only depend on the first five coordinates. Then a comparison 
with the numerical value computed by the SOS algorithm (as stated in Rogers [27J 
sadly yields a 3.64 % difference after 9.87 seconds of computations on a 1.4 GHz 
Personal Computer (whereas 50S needed 14 seconds on a 600 MHz PC). 

More extensive numerical experiments (on computers of better performance) may 
find, however, that a cubature-based algorithm is superior to a Monte-Carlo routine 
when higher dimensions than d = 7 are considered. On a different note, recall that in 
practice for the vast majoriy of derivative options, pricing algorithms are only used as 
part of hedging programs - and with hedging, the accuracy of the prices computed is 
of lesser importance than the processor time the algorithms actually requires. 

We conclude this Chapter by stating some Bermudan and American option prices 
computed through our cubature-based algorithm. 

For a min-put on a basket of seven independent assets with discount rate r — 0.06, 
maturity at time 0.5 and strike price K — 100, one will get the following numerical 
results. (Here, extrapolation I is the extrapolation of ?7™°(T)(a;) from h G {t' T' i"} 
to ft. = with scaling exponent a — 1.0, and extrapolation II is the corresponding 
extrapolation with scaling exponent a — ^. The amount of time elapsed during each 
computation is given in seconds.) 

If all volatilities txi, . . . , (Xd are equal to 0.4, then 
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start prices 


Bermudan 
= 3) 


rjxtrapolation 1 


iijxtrapolation 11 


± 11116 


80, . . . , 80 


42.7981 


44.4427 


45.7474 


12.35 


90, . . . , 90 


36.0074 


37.8231 


39.2082 


11.52 


100 100 


29.2172 


31.2058 


32.6757 


11.81 


110,..., 110 


23.2196 


24.2340 


25.3323 


12.07 


120, . . . , 120 


17.3103 


17.2211 


17.8944 


11.78 


80,90,90,100, 
110,110,120 


34.1646 


34.1625 


34.2569 


12.89 


case the volatilities ai,. . . ,ad are all equal to 0.6, then 




Start prices 


Bermudan 
= 3) 


Extrapolation I 


Extrapolation II 


± imtj 


80, . . . , 80 


53.0963 


55.1235 


56.2195 


11.94 


90, . . . , 90 


47.5948 


49.8475 


51.0137 


12.49 


100, . . . , 100 


42.0951 


44.5798 


45.8316 


12.18 


110, ...,110 


36.9143 


39.1641 


40.3082 


11.90 


120,..., 120 


32.2508 


33.5128 


34.2440 


11.60 


80,90,90,100, 
110,110,120 


45.9639 


45.8861 


44.4914 


11.56 


nally, if ai — 0.3, 


0-2 = 0.4, (73 = 


0.5, (74 = 0.6, (75 = 


0.7, erg = 0.8, 0-7 


= 0.9, one 



has the following figures: 



Start prices 


Bermudan 

(f-3) 


Extrapolation I 


Extrapolation II 


Time 


80, . . . , 80 


56.6740 


57.1258 


55.2092 


11.50 


90, . . . , 90 


51.6227 


52.1136 


49.9160 


13.04 


100, . . . , 100 


46.5720 


47.1037 


44.6290 


11.51 


110, ...,110 


41.5507 


42.1357 


39.3828 


11.81 


120,..., 120 


36.7872 


37.2093 


33.9543 


11.52 


80,90,90,100, 
110,110,120 


45.1416 


44.8038 


42.5739 


11.51 



These data suggest that the optimal scaling exponent a for the extrapolation from 

Bermudan to American min-put prices will have to depend on both the volatility vector 

(<Ti, . . . , ad) and the vector of (logarithmic) start prices (xi, . . . , Xd). 

Our second example concerns itself with the pricing of Bermudan and American 

put options on the arithmetic average of a basket of independent assets. 

. , Gd are equal to 0.4, then 

Bermudan ^ , . ^ 

Extrapolation I 



If all volatilities ai , 
Start prices 



(^ = 3) 



Extrapolation II Time 



90, . . . , 90 10.0000 10.0000 

100, ...,100 3.15446 3.13173 

110,..., 110 0.715959 0.681342 

120, . . . , 120 0.113121 0.149574 

80,90,90,100, 
110,110,120 
In case the volatilities (7i 



3.17743 3.22211 

, (7d are all equal to 0.6, then 



10.0000 

2.44775 
0.793472 
0.214042 

2.70252 



10.52 
10.68 
10.24 
10.86 

10.44 
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start prices 


Bermudan 
= 3) 


rjxtrapolation 1 


rjxtrapolation 11 


± ime 


80, . . . , 80 


20.0000 


20.0000 


20.0000 


10.18 


90, . . . , 90 


10.8834 


10.8347 


10.4006 


10.62 


100, . . . , 100 


5.34809 


5.16546 


4.54056 


10.30 


110,..., 110 


2.22539 


2.35572 


2.18241 


11.09 


120, . . . , 120 


0.881495 


0.758719 


0.791531 


10.85 


80,90,90,100, 
110,110,120 


5.38701 


5.21894 


4.58038 


10.67 


AnH if (Ti — 




5 ith — 6 iTc — 


7 rra — R it^ 


= 0.9, 


obtain the following 


results: 








Start prices 


Bermudan 
= 3) 


rLxtrapolation 1 


rjxtrapolation 11 


± ime 


80, . . . , 80 


20.0000 


20.0000 


20.0000 


11.32 


90, . . . , 90 


11.3241 


11.4565 


11.4368 


10.55 


100, . . . , 100 


5.69733 


5.67916 


5.47655 


10.98 


110,..., 110 


2.44459 


2.52824 


2.35541 


10.52 


120,..., 120 


0.937346 


0.915058 


1.10870 


11.04 


80,90,90,100, 
110,110,120 


6.24348 


6.26401 


6.26811 


10.33 



The first line of each of these sets of figures of course simply means that immediate 
exercise is optimal if the start price of each asset is at 80 or below (and, in case cr, = 0.4 
for alH e {1, . . . , 7}, even if each asset start price is at 90). 

Again it is apparent from these numerical data that the scaling exponent needed for 
the extrapolation from Bermudan to American prices put-on-the-average option prices 
has to be varied with the vector of asset start prices and possibly the volatilities of the 
underlying assets (otherwise the American price computed by extrapolation would be 
at times very significantly below the approximate price of a Bermudan option on the 
same basket and with the same payoff function) . 
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Part V 



High-dimensional approximate 

A-hedging 
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Chapter 10 

Hedging options on multiple 
assets - a suggestion for further 
research 



10.1 Theoretical suggestions 

Up to this point, the subject of our investigation has been the pricing of high-dmensional 
American and Bermudan options. In practice, there is at least as much (if not even 
significantly more) interest in the heding of such options as in finding out their price 

- the latter task often being simply left to the markets. To this extent, any pric- 
ing algorithm gains much of its practical interest merely from being employable as a 
subroutine of a hedging algorithm. 

The canonical way of hedging - that is replicating a portfolio, ideally without risk 

- that does not need to introduce utility functions for portfoHos which sometimes may 
not be that easy to justify themselves is A-hedging. Unfortunately, however, there 
is no straightforward multi-dimensional generalisation of A-hedging in the discrete 
binomial model (in the sense of eg Hull or Wilmott, Howison, Dewynne [SI)- For, 
if the price processes of all of the assets in a portfolio of d different types of shares 
each follow the binomial model, then at each time step n where the vector of current 
asset prices equals Xn S K'' there are 2'* possible states of the market that may be 
encountered at the next time step (given by a set of the form 



as each asset i G {1, . . . ,c?} is assumed to move either by a factor > or by a 
factor Pi > where without loss of generality one may assume a,; Pi for all i) 
compared with only d+1 elements in the portfolio (including the bond). Assuming 
translation-invariance of the Markov chain, we introduce the notation 



Then the volatility ai of the z-th asset is defined to be the square root of the variance 
of the one-dimensional random walk with steps ai , Pi and transition probabilities 







and 1 — Pi respectively, on the set a^No + PiNg. This is to say. 
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We shall define for each such vector y in the set 

Y -.^{zeW : Vz e {1, . . . , d} e ft}} 

the overall absolute correlation by 

(\ ~ [Oi^ ■ Pi + ■ {'^ - Pi))\ 

piv) -iV fpr— ■ 

One can now think of various approximate A-hedging algorithms - previsible trans- 
action policies that whilst being unable to eliminate the A altogether, reduce it sig- 
nificantly. At least three classes of such algorithms come to one's mind: 

1. At each time step A-hedging of proper subsets of the portfoHo, possibly changing 
the subset with time. 

2. Removing 2'^ — [d+l) of the elements of Y (thereby making the market model 
a "d -t- 1-nomial" one) via a correlation analysis (cf Section FlO. 211 . 

3. The use of cubature formulae to achieve this elimination. 

A natural method of comparing these hedging algorithms will be to look at the 
£P-norms of the resulting sequences of A's. The algortihms ofQ] and El are straight- 
forward modifications of standard A-hedging algorithms for d-component portfolios in 
market models that only allow for d+1 possible states of the market at the respective 
subsequent time step. 

We will therefore dedicate the rest of this short Chapter to the algorithm suggested 
in El 



10.2 A A-hedging algorithm based on a correlation 
analysis 

Following suggestion[2l of the preceding paragraph, we will now propose an algorithm 
that constructs a new set Y' with cardinality d -I- 1 from Y (in the notation of the 
previous Section). It will then be possible to apply A-hedging to the Markov chain 
market model that is given by 

In order to describe the algorithm that produces Y' from y, two cases according 
to the size of d have to be distinguished. 
Case r. d < 3. 
In this situation, 

|r \ r'l = 2'' - (d + 1) < d + 1 = |r'| 

and one will determine the set of those 2'* — (d + 1) vectors in the set y C M'' which 
will be removed, that is the elements of y \ F' (as opposed to finding the elements of 
Y' themselves). 

The set Y \ Y' will comprise exactly the d + 1 elements yofY with the smallest 
overall absolute correlation p{y). 

Note that a priori there can he y z Y such that p{y) = p{z). Therefore, in 
order to get the procedure of constructing Y \Y' from Y well-defined, it is necessary 
to first define a well-ordering ^ on y and to define a Hnear order <p on Y by 

yy,zeY y <p z -.-^ p{y) < p{z). 
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Then, the product {<p x ^) will be a well-ordering and we will define Y \ Y' to be 
the set of the (<p x ^)-smallest 2'' — (d + 1) elements of Y. 

Case IT. d > 3. 

In that case 

|y\y'| = 2''-(d + i)>rf + i = |y'|. 

For this reason it is faster to single out the elements of Y' directly, rather than deter- 
mining the elements of its complement Y \Y' first. 

Using the well-ordering defined above, one will thus choose the (<p x ^)-greatest 
d+1 elements of the finite set Y (the elements of Y with the largest overall absolute 
correlation p, that is). 
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Part VI 

Appendix 



Appendix A 

Re-formulation of the perpetual 
Bermudan pricing problem in L 
and 



A.l Non-applicability of the L^(M^) Spectral Theorem 

Consider a c?-dimensional Levy basket X with associated family of risk-neutral prob- 
ability measures P' and discount rate r > 0. 
Fixing h> Q and defining 

P ^L^(E^\G), A A" := I - e-'^''Px„-x, * • = (^o - g-^'^Pxo-xJ * 

we can rewrite the result of Lemma [231 as follows: 

PA{vlj-gi)=-PAgi (A.l) 

where we assume that g has a square-integrable extension from G to the whole of 
M^; given this assumption, the gi G L'^{W^) of the previous identity can be any such 
extension. 

We will suppress the superscript of A for the rest of this paragraph. 
Also, without loss of generality, we will assume in this Chapter that the components 
of the basket X when following the Black-Scholes model all have volatility 1. 

Lemma A.l. Let X he a Levy basket with associated family of risk-neutral probability 
measures P' and discount rate r > 0. Then A and PA \ ^^(Rrf \ g) are invertible. 

Furthermore, the norm of A is bounded by (l + e^'''')^, if X — fi ■ +B (thus 

/i=(r — where B is a standard Brownian motion. Moreover, A is a contraction 

^f^l^O. 

Proof. Suppose 7^ u e L'^{R'^\G) and u is bounded. Then ess sup \u\ ^ and we may 
choose a set -ff C G of positive Lebesgue measure such that e~''''ess sup |m| < \u{x)\ 
for all X G H (this is possible because r,h > and therefore e"'''' < 1), we deduce 

yxe H ler'^'Txa-x^ * u{x)\ < e^'^'^ess sup |u| 

< u{x), 

which means that 

\/xe H PAu{x) = u{x) - e-^'^Fx^-x^ * u{x) ^ 0, 
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hence PAu ^ (for H has positive Lebesgue measure). So 

kerPA \ L2(Rrf \ G) = kerPA n L2(R'^ \ G) = {0} 

and we are done for the invertibiUty of PA \ L?{^'^ \ G). Similarly, one can prove the 
invertibility of A. Finally, A is seen to be a contraction by appHcation of the Fourier 
transform: The Fourier transform is an isometry (by Plancherel's Theorem), thus 



-f\ 



/) 



/ 

\-\''h/2 



1 - e" 



-rh-\-ih fix—- 



L2(]Rd) 



fix) 



dx 



Now, the factor in front of 



f{x) in the last Hue can be bounded by (l + e , and 
it is strictly less than one for /i = 0. Using Plancherel's Theorem again, this yields the 
result. □ 



Now, this is sufficient to apply a Wiener-Hopf factorisation (for a general treatment 
of this kind of factorisations, one may consult eg Speck j2Hl, our application uses in 
particular |2Sj 1.1, Theorem 1]) and state 

Theorem A.l. Let G Q and let X he a Levy basket with associated family of 
risk-neutral probability measures P' and discount rate r > 0. Then Vq, the expected 
payoff of a perpetual Bermudan option for G with exercise mesh size h > and payoff 
g, is - using the above notation - given by 



Vg = ,91 - {PA \ i'(M'^ \ G)) ' PAg, = .gi - A^^PAZ^PAg, 
where A = A^A^ is a Wiener-Hopf factorisation of A. 
We observe 



Lemma A. 2. The Hilbert space operator A : L'^{W^, C) L'^ 

Proof. We define p := er^^—^j^j^ (where A'' is the c?-dimensi 
and via the Fubini Theorem one has for every f,gG L^(R'^) 



{Af,g) = {5o-p){x-y)f{y)dyg{x)dx 



fiy) i^o ~p){x- y)g{x)dxdy 



= (/,(<5-po(-I)*5)), 

that is 

A* = (,5o-po(-I))*. 
But since the convolution is associative and commutative, this implies 

A* A = (<5o-po(-I))*(5o-p)*- 
= {So - p) * {So - P o {-'f)) * ■ 
= AA*. 



□ 
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However, it will not be possible to find a basic system of eigenvectors and eigen- 
values for this operator, since 

Lemma A. 3. The operator A fails to be compact. 

Proof. Any normalised basis provides a counterexample for the compactness assertion. 

□ 

Therefore, the equation ljA.l|l cannot easily be applied to compute the expected 
option payoff by means of a spectral analysis. Thus, our examination of the Hilbert 
space approach in the second part of this Chapter has led to a negative outcome. 

However, one can also conceive of the operators A'^ as operators on the Banach 
space (R'') : 

A. 2 The operator equation: analyticity in the ex- 
ercise mesh size 

From now on, h will no longer be fixed and we will therefore write A'^ instead of A. 

If we now assume gi to be an integrable extension of G to the complement of CG 
as an element of Quite similarly to lA.ll we can prove 

Theorem A. 2. Let G C and let X be a Levy basket with associated family of 
risk-neutral probability measures P' and discount rate r > 0. Then Vq, the expected 
payoff of a perpetual Bermudan option for G with exercise mesh size h > and payoff 
function g, is - using the above notation - given by 

V^ = gi- {PA \ L^R" \ G)y' PAgi = 51 - A'^PAZ^PAg^, 
where A = A^A^ is a Wiener-Hopf factorisation of A. 

It suffices to observe that A'' is - due to the norm estimate for the convolution 
of two integrable functions (as the product of the norms of the convolved functions) - 
also a bounded operator on L^ (R'^) . 

We shall now identify g and gi. 

Theorem A. 3. With the notation previously introduced, we define E to be the semi- 
group 

where 

g^,,, = (27ri)-'^/2e-lM*-lV(2t) 

is the distribution of the logarithmic price vector at time t. Suppose x ^ G and, with 
the notation from the previous chapters, X is a (normalised) Brownian motion with 
(possibly zero) drift (TBlack-Scholes model/ Then h 1-^ Vq{x) is real analytic in h on 
(0, +00) as function with range in the Banach space L^{R'^). 

Proof. It is obvious that i5 is a semigroup. According to ^Ol Theorem 1.48], the set 

L\R'^) : t^Etf E L\R'') entire} 

is dense in L^{R'^). Hence it is possible to approximate every g by a sequence {gk}k C £ 
in L^{R'^). Since 

yt>0 ||i?t|lLi(E.)<e-'-*<l, 

we obtain Etgk — > Etg for fc — > cxd uniformly in t on R+, where Etgk is entire for every 
t > and k G N. Thus, 1 1-^ Etg, and thereby 1 1-^ Pg — PEtg, is an analytic function 
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on (0, +00) taking values in the Banach space L^{R'^). Now observe that for arbitrary 
open U CC M+ (the symbol "CC" indicating that U is contained in a compact subset 
of M>o) the following equations hold: 

yteU = g-{l-PEt\ PL\R'')y' (Pg-PEtg) 

= 9 (fyPEt)'^ (Pg - PEtg) 

00 oc 

= g-Y,{PE^)'^Pg + Y.^PE^f+'g 

00 

= Y.iPEt)" (g ~ Pg) , (A.2) 

k=0 

since the sums converge uniformly int onU CC M+, yielding the analyticity oft 1-^ Vq 
as a function whose range lies in the Banach space L^(R'^). □ 

Lemma A. 4. Let u > 0, n G N. Then the equation 

d" / I \" 

holds (where *y denotes the transpose of a vector y). In particular, i/^A/ — /iV/ — Xf 
for some A, /, 

—Euf = {-r + \Tgt^,t*f- 

Proof. According to Davies 'JH, Proof of Theorem 2.39], we have 

111 

—E^^{ZE^/.^)\ (A.3) 

where Z denotes the infinitesimal generator of the semigroup E. Now, define C to be 
the convolution operator semigroup {gtfj.,t * ■)t>o (normalised) Brownian motion with 
drift /X (as before denoting by g^^cr^ the Lebesgue density of the Gaussian distribution 
centered around z of variance cr^ for all z e and a > 0). It is well-known (cf e g 
|2()[ p. 352]) that the infinitesimal generator of this semigroup C is 

L := ^A + VV. 

By our requirements on /, Lf — on U. Furthermore, L and C commute: 

> CtL = LCt. 

Thus, 

yt>0 ZEt = -.Et^^ie""' -Gt) 
d dt^ ' 

_ —rt/~i I —rt ^ /~i 

= -re Of + e — 6t 
dt 

= e~-*Ct{-r + L), 

which due to equation IIA.3p already suffices for the proof of the Lemma in the general 
case. And if / is an eigenfunction of L for the eigenvalue A, one has (— r + L)"/ = 

i-r + xrf. 

□ 
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Theorem A. 4. The Taylor series for the expected payoff of a perpetual Bermudan 
option as a function of the exercise mesh with respect to a fixed exercise region G is 
for all s > 0; 

'^■=0 h + --- + Irn = n ^'=^ ' ' 

ai,...,Z™)eNo™ 



where, in order to avoid confusion with pointwise exponentiation, A°'' denotes A'^ for 
any operator A. 

Proof. We know about the real analyticity of 1 1-^ Et on and even, thanks to the 
previous Lemma, the expUcit Taylor series. Thereby we also have the Taylor series for 
t h-^ PEt . So we can use equation IIA.211 and see by means of a binomial expansion 



A:=0 k=Q \ e=0 



oo oo 



EE E U^-^e-Pi-r + LY-C. 

Gi,...,Z™)eNo™ 



E(*-^rE 



{h,...,l^) gNo 



ll + ---+lrn=n 



□ 



This Taylor series fails to provide any straightforward possibility for the computa- 
tion of Vq- Instead we state the following immediate Corollary of equation IIA.211 : 

Corollary A.l. With the notation as in the previous Theorem, 

m— 1 



106 



Appendix B 



An algebraic perpetual 
Bermudan pricing method and 
its natural scaling 

In this Chapter, we will, in a more algebraically flavoured way, present an approach 
that approximates the perpetual Bermudan option price as the fixed point of some map 
on the space of polynomials that is defined by means of not only the max operator, 
but also interpolation with respect to a given, fixed, set of interpolation points, as well 
as convolution with one and the same Gaussian (not necessarily normalised) measure. 
This set of interpolation points could, for example, be a set of cubature points for the 
distribution of the time h increment of the logarithmic price process (if this process is 
assumed to be Gaussian with stationary increments). 

Let us, for this purpose, adopt Victoir's notation and denote the space of all 
polynomials of degree m and degree at most m by Rm[X] and ]R<m[-'^], respectively, 
for all m S N. We will write polynomials in the form p[X] and denote by p{-) the 
associated polynomial function from R to M. 

Now, let TO be any positive integer. We introduce the interpolation map 

X : X M™+i ^ M<™[X] 

that assigns to each pair (/, x) of a real- valued function / and a vector x of interpolation 
points the well-defined (see e g Stoer and Bulirsch (5^1) polynomial p{X) G Rm[X] 
satisfying 

Vje{0, •••,m} p{xj) = f{xj). 

Now let r > be a positive real number (interpreted to be the discout rate), 
/i := ^ — r e M (which is the drift of the logarithmic price process in a Black-Scholes 
model with normalised time scale - ensuring the volatility a to be equal to one), and 
(7 : M ^ M a continuous function (interpreted to be the payoff function of an option 
defined on the space M of logarithmic underlying asset prices) . Consider a vector- valued 
function y : R>o ]R'"+^ such that the range of the function y consists exclusively of 
vectors with mutually distinct entries. The elements of the range oi y can in this case 
serve as sets of interpolation points (these points also called support abscissas). Hence 
using the notation of previous paragraphs, we may define another map 

for all /i > by 

Hh : P[X] ^ I ((p(-) * u-,,h.h ■ V 5, m) ■ 
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Note that 

This can be shown using the hnearity of the convolution and the fact that for all 
k E N, the convolution v-f_ih,h * {'Y of the measure V-ij,h,h with the function x x^ 
is again a polyonmial function of degree k. Writing down an explicit formula for the 
map V-fj,h,h * • : ]Rm[^] — > we see that this function is continuous with respect 

to the Euchdean topology on the (to + l)-dimensional real vector space ]R<m[X]. 

This ushers in the proof of the following Lemma which is one of the first observations 
leading to the fixed point equation mentioned at the beginning of this Chapter. 

Lemma B.l. Consider any h > and arbitrary p[X] € 'Hh is continuous 

with respect to the Euclidean topology on the (to + I) -dimensional real vector space 
M<„[X]. Also, if and only if T-lhp[X] = p[X], that is, p[X] is a fixed point of Hh, 
there will exist a polynomial q[X] G R„[X] such that p[X] = lim„^oo {T~i-h)" q[X] (in 
the Euclidean toplogy of the (to + 1)- dimensional vector space Mm[X]). 

Proof. The continuity of Hh is a consequence of the continuity of the map V-^h,h * • '• 
—^ Mm[X]. For the second part of the Lemma observe that provided the 
existence of such a q[X] as in the statement of the Lemma, we can deduce the fixed 
point equation of the Lemma's statement from the continuity of the map Hh- For the 
converse implication, simply take q[X] =p[X]. □ 

Notational convention B.l. For any vector q = ' (qa,- ■ ■ ,qm) € ffi™^^, q[X] will 
be understood to be the polynomial q[X] = X^^g Ik ' X^ € M<f„[X], and for all q{X\ € 
]R<to[X], q{-) shall be understood to denote the polynomial function 

rn 

q{-) : X ^^qk ■ x'', 

k=0 



Lemma B.2. For all h > and each p[X] G R^l-'^], the function {x,h) i-^ p(-) * 
V-nh,h{x) is a polynomial function in both x and h. Its degree in x is m, its leading 
coefficient in x being the leading coefficient of p[X]. Furthermore, the function h 

p{-) * v-fih,h{x) is o{h) + for all a; € K. Moreover, if x = o {h^^, then the first 
two leading terms of {■Y' *V-nh,h{x) as a function of h are x'' = o (h^^ O'l^d a term of 
order o ^ respectively, if k is odd - and if k is even, the first two leading terms 

of h<-^ {■)'' * I'-fxhthix) are x^ = o {h^^ and some term of order o {h^^ , respectively. 

Proof. First, let us once again remark that y-nh,h*- is, as a convolution operator, hnear 
and that we therefore may restrict our attention to the functions V-f^hM * {x i-^ x'^) 
for k <m. Consider any xq S K. Using the transformation y ^^^^^^ we find that 

J^-ij.h,h* {x ^ x^){xo) = j (^xq + nh - zV2h}j e'"^ dz 

= x'' + polynomial terms in .xo°, • • • , x^^^^ , h, 
where the last line follows from expanding the binomial and using the identity 

e No / dz = Q 

JTS. 



(an immediate consequence of the "oddness" of the integrand) , which entails that all 
odd terms in /i^/^ will be cancelled out (as they have to be odd terms in z as well). □ 
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Remark B.l. The assumption of x = o {h^^ will hold in particular for any non-zero 
element of a set of cubature points for the measure v^^h.h derived from a cubature 
formula for i^o.i • 

Definition B.l. For all h > 0, a <E Sm+i (the symmetric group 0/ {0, ... , m} ), and 

n < m, we define the (m + 1) x (m + 1) -matrix 



I 



\ 



+ya(n+l){hf 



+ya{rn){hr J 



and the following vector: 



( 9{ya{o){h)) \ 

9 {ya(n){h)) 




V {y, n, h, a) . 



V / 

We also define p{y,n^h, a) to be the solution q of 

^nW ■ 9 = v{y,n,h,(T) , 

provided there exists a unique solution to this equation. 

From the matrix formulation of the interpolation problem (again cf Stoer and 
Bulirsch (22]), the following Lemma is immediate: 

Lemma B.3. The polynomial q[X] = J2k=Q li^-^'' ^ l^m[-'^] is a fixed point of Hh if 
and only if there is an n < m and a a E Sm+i (the symmetric group of {0, . . . ,m\) 
such that 

KXh) • *(qo, ■ • ■ ,9m) ^v{y,n,h,a) 



ana, m a. 



well 



Vje {0,...,n} 

9{ya(j){h)) 

Vj e {n, + 1, . . . , m} 

9 {y.U){h)) V e-'-'' {q{h){-) * i^.^hm) {y.u){h)) 



-rh 



{q{h){-)*v^t,h,h) {yc(j){h)) . 



Lemma B.4. Suppose there is ana E (0, 1) such that for all j £ {0, . . . , m}, h ^ yj{h) 
is a non-constant polynomial in h" , except for possibly one ja G {0,...,to} where 
yj„{h) = for all h > 0. Then for all sufficiently small h, detA^Jh) ^ 0. The upper 
bound in M>o U {+cxd} on all those h that satisfy the previous inequality detyl^(/i) 7^ 
for all n < m and a G Sm+i shall be denoted by Hq (y). 
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Remark B.2. According to B.emark \B.1[ the assumption of h yj{h) being non- 
constannt polynomial in h" for all j G {0, . . . ,m} (apart from possibly one zero coor- 
dinate) for some a holds in particular for any set of cubature points for the measures 
v-nh,h that is derived from a cubature formula for the normalised Gaussian measure (in 
that case a ~ ^). Similar assertions hold if one replaces v-^h.h by rah where {mt)t>o 
is the convolution semigroup associated to some other symmetric stable process. 

Proof sketch for Lemma \B.4\ The function h i-^ detA'^{h) is, by our assumptions on 
the functions yj on the one hand polynomial in /i", as one can see exactly as in the 
proof of Lemma FB. 21 On the other hand, one can show, using the polyonmiality in h 
and the assumpion that all the entries of y{h) are mutually distinct for all h > 0, that 
the function h i-^ detA'^{h) is non-constant. Hence, h i-^ detA'^{h) is a non-constant 
analytic function in /i", and note that hi-^ h"' is a bijection on the unit interval (0, 1). 
Therefore, h det Af^{h) cannot be constantly zero on the open unit interval (0, 1), 
but it can also only have finitely many critical points on that interval. Hence there 
must be an /i > such that either detA^{h) < for all h e (O, h) or detA^{h) < 
for all he {0,h). □ 

On the other hand, we have the following 

Lemma B.5. Suppose the real-valued function h i— > yj{h) is a non- constant polynomial 
in h i—> h°' for all j G {0, . . . ,m} for some a e (0, 1), except for possibly one jo G 
{0, . . . ,to} where yj„{h) — for all h > 0. Furthermore, take g to be analytic. Then 
for all n < m and all a G Sm+i there exists a vector TZ{n, a) G {<, >j.'"+i of relations 
such that for all sufficiently small h > 0, 

Vi G {0, . . . , m} 

(g - e-~''^P {y, n, h, a) (•) * i^^f^hji) (yz(h)) 7^(cr, n)j 0. 

The upper bound in M>o U {+00} on all those hi such that for all h < ho (y) A hi the 
relations in the previous line hold for all n < m and a G Sm-i-i shall be denoted by 
hi (y). (Here ho (y) is the strictly positive constant of Lemma \B.4\ ) Thus, hi (y) < 
ho (y)- 

Proof. We have already defined p[y, n, h, u) to be the solution q of 
^nih) ■ * ((7a(o), ■ • • ,9<7(m)) = v{y,n,h,a) . 

From Cramer's rule, Lemma lB.2l and our assumptions on y (coordinatewise polynomial 

in (•)") as well as g (analyticity) , we derive that h^ [g — e~''''p (y, n, h, a) (•) * I'-^h.h) {yj{h)) 

is analytic in (•)" and therefore has only finitely many critical points on (0, 1) for all 

j G {0,...,r7i} and arbitrary choice oi y,n,a. Thus, when approaching zero, these 

functions must eventually stay on either side of nought. Put more formally, there 

must be for all a G Sm+i and n < m a vector TZ{n, a) G {<, >}™+i such that for all 

sufficiently small h, and for all i G {0, . . . , m}, 

{g " e"''''p (y, n, h, a) (•) * V-^h,h) {yi{h)) TZ{(t, n)i 0. 

□ 

Corollary B.l. Let the assumptions of the previous Lemma hold. If there is a fixed 
point of Tihi for an hi < hi{y) (the strictly positive constant of Lemma [_B.,5|) . all 
T-Lh with positive h < hi must have fixed points as well. There exist a permutation 
do G Sm-i-i as well as a natural number no < m such that the coefficient vectors 
q = *{qo, ■ ■ ■ ,qm) to all these fixed points q[X] — ^^^olkX^ are solutions q to the 
linear equation 

^nl^h) ■ q^ v{y,no,h,ao) 
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and satisfy 

Vje{0,...,no} 

9 {yMj){h)) > e-''^ {q(h){-) * v^^^hm) {yaoU)ih)) 

as well as 

Vj e {no + 1, . . . ,m} 

9 (y<ToO)(/i)) < e^'''' * l^-^,h.h) {yao{j){h)) ■ 

Proof. By virtue of Lemma IB. 81 a fixed point is a polynomial q[X] = X^I-Lo Ik^'' 
whose coordinate vector q = ''^{qo, . . . ,qm) solves the linear equation 

^nW ■ q = v{y,n,h,a) 

and satisfies, moreover, inequalities of the form 

{g~e-^''q{-)*,y^^h,h) {y^{h)) 7^, 

for some vector of relations TZ £ {<, >}™. Now apply the previous Lemma FB. 51 

□ 

Theorem B.l. Suppose there is a vector G such that yj{h) = —^h + £^j ■ h^/'^ 

for all j e {0, . . . , m} (where the mutual distinctness of the entries of y{h) for all 
h > entails that S,o, - ■ ■ are mutually distinct as well), and assume furthermore 
that g is analytic and satisfies g{0) ^ 0. Suppose, moreover, that there exists a fixed 
point of 'Hh2 for some /i2 G (0, /ii (y)) (where hi (y) is the strictly positive constant 
from Lemma TB.!^) . Then by Corollarv \B.l\ the maps Hh do have a fixed point for every 
h < /i2, and let n,a be the natural number < m and the permutation whose existence 
is stated in Corollarii \B. 1[ respectively. Also write p{h)[X] :— p{y,n,h,a) [X] for all 
h G (0, /i2). Then there is a polynomial p'^[X] such that 

p(ft)[x]=/[x] + o(;i^) 

componentwise. 

Remark B.3. Such a ^ £ M™+^ exists in particular whenever {ya, . . . , ?/„i} is a set of 
cubature points for the measure v^^h,h derived from a cubature formula for t'o,!. 



o(/i^) 
(l-e''''') •o(/i^) 

(l-e-''^) •o(/i^) J 



Proof of Theorem. \B.l\ Observe that by Lemma lR2l 
Am 



( 



1 



o\h 



oh 



/ 1 o(h^ 



1 
o ( 



o I /l 2 
O ( /l3 



^ o o (^h2 



o{h'T) 



oh 
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Next, we will use Cramer's rule to determine if there is a limit for the solution of 



Al{h)-p{h) =v{y,n,h,a) 



as h tends to zero and if so, what the convergence rate will be. For this purpose, 
we have to consider the determinant of the matrix A'^^^{h) which is defined to be the 
matrix coinciding with A^(ft) in the columns 0, . . . ,i — l,i + 1, . . .m and having the 
vector 

/ 9{yo{h)) \ 



givnih)) 





V J 

as its i-th column. Then, since g is right-differentiable in and by assumption g{0) ^ 0, 



o[hi) 




/ o(/.S) 



O /l 2 




V 



o h 2 



[h^^^ o(l) o(/i^) 



O ( ft. 2 

o (h^ 



o { h 2 



o(l) o(/i^ 
o(/i^^ 



o ( hr^ 



o{h'i) 
o {h^^ 



o 



I 



From this, we can conclude that for all z e {0, . . . , to} there is a constant G M such 
that 

det<^,(fe) / 



det Al^ih) 



Example B.l (m = 2, the quadratic case, when = 0). H^e /ia?;e 
V/u G M Va: G M 



□ 



implying 



Va;G 



{■f * y-^h,h{^) =x^ + 2fj,hx + h + h'^ 



{■f*'^o,h{x) =x'^ + h + h'^ 
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and 

Wh > m^l ^i-h^ ■ 
\a-h-^ J 

Suppose the parameters from Corollary \B. 1\ are in our example — 1 and (To — id. 
Furthermore, in our case, for all h > 

( 1 iohh ^o^h 

Af{h) = 1 ^^h 

Therefore for all h > 0, 

det A\^{h) = ^ihi {{1 ~ e''-'') ^2'' ~ e-^'\l + h)) 

Ho^i'hi (1 _ e--'^) + (1 - e---") ^o'^/i^ 
-Co/i*((l-e-'-") ^,^^e-^\l + h)) 

= (1 - e-"-") hi + eoa' + ^0^6 - eo6' - - ^o'Ci) 

TVoie that 

66' + 66' + 6'6 - 66' - 6'6 - 6'6 = dot b (e) , 

where B is the matrix of the interpolation problem with support abscissas 6;6i6- 
Due to the unique solvability of the interpolation problem (see again e g Stoer and 
Bulirsch /&9|/). this determinant det B never vanishes unless the support abscissas 
67616 fO''^l to be mutually distinct. 

Remark B.4. Similarly one can prove that the function h 1-^ P^[X] is differentiable 
in if one assumes the support abscissas to be polynomial in h rather than h^. 
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